Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 1,485,365 C→T 19.2% D178D (GAC→GAT hrpA → ATP‑dependent helicase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077791,485,3650C→T19.2% 158.2 ‑1.2 ‑1.7287D178D (GAC→GAThrpAATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (30/25):  ref base (106/126):  total (136/151)
Fisher's exact test strand distribution p-value = 2.93e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.83e-02

Only 100 of 297 total aligned reads displayed.

TAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGA  >  NC_007779/1485265‑1485465
                                                                                                     |                                                                                                    
tAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGAc                                                                                                      <  1:6520986/101‑1 (MQ=255)
 aaGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACg                                                                                                     >  1:10484389/1‑101 (MQ=255)
   gTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGaa                                                                                                   <  1:10587101/101‑1 (MQ=255)
     gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgc                                                                                                 <  1:900435/101‑1 (MQ=255)
     gATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAgc                                                                                                 <  1:8876595/101‑1 (MQ=255)
        aaCACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGCGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCa                                                                                              >  1:1735524/1‑101 (MQ=255)
         acacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCAc                                                                                             <  1:5412860/101‑1 (MQ=255)
          cacGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACg                                                                                            <  1:9788019/101‑1 (MQ=255)
           acGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGa                                                                                           >  1:8738791/1‑101 (MQ=255)
             gATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                                                                         >  1:850216/1‑101 (MQ=255)
             gATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                                                                         >  1:2929491/1‑101 (MQ=255)
             gATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                                                                         <  1:1908685/101‑1 (MQ=255)
             gATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAAc                                                                                         <  1:14106979/101‑1 (MQ=255)
             ctgcattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAAc                                                                                         >  1:10192952/11‑101 (MQ=255)
              tgcattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACg                                                                                        >  1:4116630/10‑101 (MQ=255)
              aTGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACg                                                                                        >  1:13928412/1‑101 (MQ=255)
                ggTCAAGCGGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCa                                                                                      <  1:11308690/101‑1 (MQ=255)
                cattaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCa                                                                                      <  1:11621582/94‑1 (MQ=255)
                cattaaGCTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGCCACCATCATCATCGATGAGGCGCACGAACGCa                                                                                      <  1:8748868/97‑1 (MQ=255)
                 attaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAg                                                                                     <  1:13327135/95‑1 (MQ=255)
                  ttaaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGc                                                                                    >  1:1628036/6‑101 (MQ=255)
                  tCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGc                                                                                    <  1:924345/101‑1 (MQ=255)
                   taaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGcc                                                                                   <  1:9898031/97‑1 (MQ=255)
                   taaacTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGcc                                                                                   >  1:11159956/5‑101 (MQ=255)
                    aaGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCt                                                                                  >  1:14525449/1‑101 (MQ=255)
                     aacTGATGACCGACGGTATCCTGCTGGCCGGAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCtc                                                                                 >  1:11981933/3‑100 (MQ=255)
                     aGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTg                                                                                 <  1:13668354/101‑1 (MQ=255)
                       cTGATGACCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaa                                                                               >  1:9707292/1‑101 (MQ=255)
                        tgatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAat                                                                              <  1:5646730/101‑1 (MQ=255)
                        tgatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaac                                                                              <  1:10810858/101‑2 (MQ=255)
                        tgatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaac                                                                              <  1:6212919/101‑2 (MQ=255)
                        tgatgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCaac                                                                              >  1:15081843/1‑100 (MQ=255)
                         gatgaCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAata                                                                             <  1:4212323/101‑1 (MQ=255)
                           tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                                                           <  1:9534727/101‑1 (MQ=255)
                           tgaCCGACGGTATCCTGCTGGCCGAAATTGGCCGCGACCGACTGTTGCGTGAGTACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATc                                                                           <  1:2057405/101‑1 (MQ=255)
                             aCCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGa                                                                         >  1:4524842/1‑101 (MQ=255)
                              ccGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAt                                                                        <  1:5847932/101‑1 (MQ=255)
                                gACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAttt                                                                      <  1:11388385/101‑1 (MQ=255)
                                 aCGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtttt                                                                     <  1:2995911/101‑1 (MQ=255)
                                 aCGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtttt                                                                     >  1:10359572/1‑101 (MQ=255)
                                   ggTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGAtttttt                                                                   >  1:14686481/1‑101 (MQ=255)
                                    gTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTg                                                                  <  1:8149494/101‑1 (MQ=255)
                                      aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                                                                >  1:4869352/1‑101 (MQ=255)
                                      aTCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCt                                                                >  1:12474847/1‑101 (MQ=255)
                                       tCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTc                                                               <  1:4576531/101‑1 (MQ=255)
                                        ccTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCg                                                              <  1:7396235/101‑1 (MQ=255)
                                         ctgctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCgg                                                             >  1:12118113/1‑101 (MQ=255)
                                           gctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCt                                                           >  1:6492332/1‑101 (MQ=255)
                                            ctgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTa                                                          >  1:6315258/1‑101 (MQ=255)
                                             tgGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAt                                                         >  1:9299394/1‑101 (MQ=255)
                                               gCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTAttt                                                       >  1:7694378/1‑101 (MQ=255)
                                                cGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTg                                                      >  1:13165848/1‑101 (MQ=255)
                                                 ggAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGa                                                     >  1:8625423/1‑101 (MQ=255)
                                                  gagaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGaa                                                    <  1:12165496/101‑1 (MQ=255)
                                                  gagaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGaa                                                    >  1:3444531/1‑101 (MQ=255)
                                                    gaTCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAag                                                  <  1:13967165/101‑1 (MQ=255)
                                                      tCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag                                                >  1:1300767/1‑101 (MQ=255)
                                                      tCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAagag                                                <  1:5469780/101‑1 (MQ=255)
                                                       ccAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGt                                               <  1:5563301/101‑1 (MQ=255)
                                                       ccAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGt                                               >  1:12944291/1‑101 (MQ=255)
                                                       ccAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGt                                               <  1:10688109/101‑1 (MQ=255)
                                                       ccAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGt                                               >  1:14465312/1‑101 (MQ=255)
                                                        cAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGtt                                              <  1:262700/101‑1 (MQ=255)
                                                         aGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtg                                             >  1:3357257/1‑101 (MQ=255)
                                                         aGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtg                                             >  1:243344/1‑101 (MQ=255)
                                                          gCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgc                                            >  1:12979355/1‑101 (MQ=255)
                                                           cAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgct                                           <  1:7736796/101‑1 (MQ=255)
                                                            aaGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTtgctg                                          >  1:13515477/1‑101 (MQ=255)
                                                                ccGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCcgc                                      >  1:2317145/1‑101 (MQ=255)
                                                                  gCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcgg                                    >  1:12810760/1‑101 (MQ=255)
                                                                   ccTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggc                                   >  1:13192783/1‑101 (MQ=255)
                                                                   ccTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCAGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggc                                   <  1:11941338/101‑1 (MQ=255)
                                                                    ctgctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg                                  >  1:6086091/1‑101 (MQ=255)
                                                                    ctgctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg                                  >  1:2790409/1‑101 (MQ=255)
                                                                    ctgctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCgcggcg                                  >  1:8578834/1‑101 (MQ=255)
                                                                      gctgATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTc                                >  1:9499171/1‑101 (MQ=255)
                                                                         gATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTg                             >  1:3278864/1‑101 (MQ=255)
                                                                              cAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCt                         <  1:2537861/101‑1 (MQ=255)
                                                                                gTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTaa                       <  1:5198445/101‑1 (MQ=255)
                                                                                gTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTaa                       >  1:2824217/1‑101 (MQ=255)
                                                                                 tACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTaaa                      <  1:7179327/101‑1 (MQ=255)
                                                                                 tACGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTgaa                      >  1:10262897/1‑101 (MQ=255)
                                                                                  aCGACACCATCATCATCGATGAGGCGCACGAACGCAGCCTCAACATCGACTTCCTGCTGGGATATTTAAAGCATATTCTGCCGCAGCGCCCAGATCTgaag                     >  1:11052337/1‑100 (MQ=255)
                                                                                    gACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAAt                   >  1:11129884/1‑101 (MQ=255)
                                                                                    gACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATc                  <  1:11445825/101‑1 (MQ=255)
                                                                                     acacTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATc                  >  1:3532770/1‑101 (MQ=255)
                                                                                      cacTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCa                 >  1:6727268/1‑101 (MQ=255)
                                                                                      cacTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCa                 <  1:14544982/101‑1 (MQ=255)
                                                                                       acTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCAt                <  1:7132019/101‑1 (MQ=255)
                                                                                          aTCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTAt             <  1:10059177/101‑1 (MQ=255)
                                                                                              ttattGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACt         >  1:1167372/1‑101 (MQ=255)
                                                                                                attGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTc       <  1:3144038/101‑1 (MQ=255)
                                                                                                attGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTc       >  1:12673729/1‑101 (MQ=255)
                                                                                                 ttGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTcc      <  1:3966340/101‑1 (MQ=255)
                                                                                                  tGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCcg     <  1:14909177/101‑1 (MQ=255)
                                                                                                  tGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCcg     <  1:12250263/101‑1 (MQ=255)
                                                                                                   gACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCcgc    <  1:7994976/101‑1 (MQ=255)
                                                                                                   gACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCcgc    <  1:5322772/101‑1 (MQ=255)
                                                                                                    aCGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTACcgcg   <  1:8069402/101‑1 (MQ=255)
                                                                                                     cGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGa  >  1:3697813/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
TAAGTGATAACACGATGGTCAAGCTGATGACCGACGGTATCCTGCTGGCGGAGATCCAGCAAGACCGCCTGCTGATGCAGTACGACACTATCATTATTGACGAAGCGCACGAACGCAGCCTGAATATCGATTTTTTGCTCGGCTATTTGAAAGAGTTGCTGCCGCGGCGTCCTGACCTAAAAATCATTATCACTTCCGCGA  >  NC_007779/1485265‑1485465

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.