Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,706,121:1 +T 5.2% coding (1004/1800 nt) lepA ← back‑translocating elongation factor EF4, GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,706,1211.T 5.2% 828.3 / 39.8 193coding (1004/1800 nt)lepAback‑translocating elongation factor EF4, GTPase
Reads supporting (aligned to +/- strand):  ref base . (90/93);  new base T (3/7);  total (93/100)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 195 total aligned reads displayed.

ACGTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAG  >  NC_000913/2705974‑2706268
                                                                                                                                                    |                                                                                                                                                   
aCGTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGcc                                                                                                                                                    >  4:637397/1‑149 (MQ=255)
 cGTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGccc                                                                                                                                                   >  1:1479025/1‑149 (MQ=255)
  gTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCa                                                                                                                                                  <  4:133556/149‑1 (MQ=255)
  gTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCCTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCa                                                                                                                                                  >  3:279473/1‑149 (MQ=255)
        aaCTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAgcgcgc                                                                                                                                            >  2:629706/1‑149 (MQ=255)
         aCTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCt                                                                                                                                           >  2:803987/1‑149 (MQ=255)
          cTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTg                                                                                                                                          <  3:468902/149‑1 (MQ=255)
           ttCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTgg                                                                                                                                         >  1:1364915/1‑149 (MQ=255)
           ttCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTgg                                                                                                                                         >  1:270547/1‑149 (MQ=255)
            tctcTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGa                                                                                                                                        <  3:213807/149‑1 (MQ=255)
              tctTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGCTCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGc                                                                                                                                      >  2:1362447/1‑149 (MQ=255)
                   aCGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTGGGCtttgg                                                                                                                                 <  2:1326354/149‑3 (MQ=255)
                   aCGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTcc                                                                                                                                 <  3:1569164/149‑1 (MQ=255)
                     gTGGTTTCAACTTCATACACTACGGTCGGCGCAGTTGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCgg                                                                                                                               >  4:105319/1‑149 (MQ=255)
                        gTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTGGGCtttggtttcc                                                                                                                            <  1:753563/149‑8 (MQ=255)
                         tttCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCa                                                                                                                           <  2:363588/149‑1 (MQ=255)
                           tCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATa                                                                                                                         <  4:930692/149‑1 (MQ=255)
                             aaCTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCgg                                                                                                                               <  3:1295700/141‑1 (MQ=255)
                             aaCTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGa                                                                                                                       >  3:1421575/1‑149 (MQ=255)
                             aaCTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGa                                                                                                                       >  1:877925/1‑149 (MQ=255)
                             aaCTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGa                                                                                                                       >  1:1435140/1‑149 (MQ=255)
                             aaCTTCATACACTACGGTCGGCGCAGTGGTGATCAGATACAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCgg                                                                                                                               >  4:1295700/1‑141 (MQ=255)
                               cTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAAc                                                                                                                     >  3:1505027/1‑149 (MQ=255)
                                ttCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACa                                                                                                                    >  3:1411547/1‑149 (MQ=255)
                                  cATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGt                                                                                                                  <  1:1288161/149‑1 (MQ=255)
                                  cATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGt                                                                                                                  <  1:361650/149‑1 (MQ=255)
                                   aTACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTg                                                                                                                 >  3:121426/1‑149 (MQ=255)
                                     acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAg                                                                                                               <  1:1523345/149‑1 (MQ=255)
                                     acacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAg                                                                                                               <  3:962908/149‑1 (MQ=255)
                                      cacTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAgg                                                                                                              <  4:389347/149‑1 (MQ=255)
                                       acTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTGGGCtttgg                                                                                                                                 <  3:868317/129‑3 (MQ=255)
                                       acTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGc                                                                                                             <  1:296547/149‑1 (MQ=255)
                                         tACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCAt                                                                                                           >  3:31473/1‑149 (MQ=255)
                                         tACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCAt                                                                                                           >  4:1225343/1‑149 (MQ=255)
                                               cGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCa                                                                                                     <  3:1375312/149‑1 (MQ=255)
                                                ggCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTGGGCtttggt                                                                                                                                >  4:663798/1‑118 (MQ=255)
                                                 gcgcAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAgg                                                                                                   >  4:1560844/1‑149 (MQ=255)
                                                     aGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGa                                                                                               <  4:502754/149‑1 (MQ=255)
                                                       tggtgATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCg                                                                                                                                >  2:1238441/1‑114 (MQ=255)
                                                       tggtgATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCg                                                                                                                                <  3:1320068/114‑1 (MQ=255)
                                                       tggtgATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCg                                                                                                                                <  1:1238441/114‑1 (MQ=255)
                                                         gtgATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGttt                                                                                           <  3:646105/149‑1 (MQ=255)
                                                           gATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTAc                                                                                         >  1:1211732/1‑149 (MQ=255)
                                                            aTCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTAcc                                                                                        >  2:1288687/1‑149 (MQ=255)
                                                            aTCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTAcc                                                                                        >  2:664061/1‑149 (MQ=255)
                                                             tCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTAccc                                                                                       <  1:1110603/149‑1 (MQ=255)
                                                               agatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAg                                                                                     >  3:123859/1‑149 (MQ=255)
                                                                gatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAgc                                                                                    >  1:679123/1‑149 (MQ=255)
                                                                 atccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAgc                                                                                    >  1:313279/1‑148 (MQ=255)
                                                                 atccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAgc                                                                                    <  2:313279/148‑1 (MQ=255)
                                                                   ccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGcgcg                                                                                 >  1:1457369/1‑149 (MQ=255)
                                                                     agatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTc                                                                               <  4:31473/149‑1 (MQ=255)
                                                                      gatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCa                                                                              >  1:429589/1‑149 (MQ=255)
                                                                         cGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACgg                                                                           <  2:720559/149‑1 (MQ=255)
                                                                             ttCACGTTCCAGACGGTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAgg                                                                       >  4:1503992/1‑149 (MQ=255)
                                                                                aCGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCtt                                                                    >  1:1021896/1‑149 (MQ=255)
                                                                                aCGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCtt                                                                    >  4:114842/1‑149 (MQ=255)
                                                                                 cGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTc                                                                   <  1:1335665/149‑1 (MQ=255)
                                                                                  gTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGGGTGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACTCAGCGCGTCAAGGAAGGCTTCa                                                                  <  2:198077/149‑1 (MQ=255)
                                                                                   ttCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCAt                                                                 >  2:1102432/1‑149 (MQ=255)
                                                                                       aGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCa                                                                                                                           <  4:1428157/87‑1 (MQ=255)
                                                                                        gACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATa                                                                <  3:1116153/145‑1 (MQ=255)
                                                                                        gACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATa                                                                >  4:1116153/1‑145 (MQ=255)
                                                                                        gACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAgtcg                                                            >  1:546383/1‑149 (MQ=255)
                                                                                         aCGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCTCCAGCGCGCTGGGCtttggtt                                                                                                                               <  2:142653/81‑5 (MQ=37)
                                                                                              ccTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGaa                                                      <  2:278324/149‑1 (MQ=255)
                                                                                              ccTGGATGATCTCCATATGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGaa                                                      <  3:1337890/149‑1 (MQ=255)
                                                                                                tGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGtcgtcg                                                         <  4:1346964/144‑1 (MQ=255)
                                                                                                tGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGtcgtcg                                                         >  3:1346964/1‑144 (MQ=255)
                                                                                                    tgatCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTAcc                                                <  1:863870/149‑1 (MQ=255)
                                                                                                     gatCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCg                                               <  2:429589/149‑1 (MQ=255)
                                                                                                      atCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGaa                                                                         >  4:1532353/1‑122 (MQ=255)
                                                                                                      atCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCgg                                              <  3:114842/149‑1 (MQ=255)
                                                                                                         tcCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGaa                                           <  2:1451949/149‑1 (MQ=255)
                                                                                                           cATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACa                                         >  3:365003/1‑149 (MQ=255)
                                                                                                           cATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACa                                         >  3:1461783/1‑149 (MQ=255)
                                                                                                                  agcagGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGc                                  <  2:877925/149‑1 (MQ=255)
                                                                                                                    cagGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGt                                >  3:139478/1‑149 (MQ=255)
                                                                                                                     agGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACa                                         >  1:989166/1‑139 (MQ=255)
                                                                                                                     agGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACa                                         <  2:989166/139‑1 (MQ=255)
                                                                                                                       gCCGAGGCAGCCGCAGCGGAAACCAAAGCCCAGCGCTCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGtata                             >  4:181217/1‑149 (MQ=255)
                                                                                                                       gCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGtata                             >  3:787198/1‑149 (MQ=255)
                                                                                                                       gCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGtata                             <  4:565474/149‑1 (MQ=255)
                                                                                                                          gAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCt                          <  4:859181/149‑1 (MQ=255)
                                                                                                                             gAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCa                                                                                                     <  3:1363537/71‑1 (MQ=255)
                                                                                                                                 ccGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGttt                   <  1:746862/149‑1 (MQ=255)
                                                                                                                                 ccGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGttt                   <  2:945450/149‑1 (MQ=255)
                                                                                                                                 ccGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGttt                   <  4:141414/149‑1 (MQ=255)
                                                                                                                                        ggAAACCAAAGCTCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTc            <  4:1313733/149‑1 (MQ=255)
                                                                                                                                        ggAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCg                                               >  2:59309/1‑114 (MQ=255)
                                                                                                                                        ggAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCg                                               <  1:59309/114‑1 (MQ=255)
                                                                                                                                         gAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCt           <  1:283064/149‑1 (MQ=255)
                                                                                                                                           aaCCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTa         >  2:1146580/1‑149 (MQ=255)
                                                                                                                                           aaCCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTa         >  2:524882/1‑149 (MQ=255)
                                                                                                                                             ccAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTaaa       <  3:220844/149‑1 (MQ=255)
                                                                                                                                             ccAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTaaa       <  2:1496548/149‑1 (MQ=255)
                                                                                                                                              cAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGGCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAg      >  2:100974/1‑149 (MQ=255)
                                                                                                                                               aaaGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGc     >  2:454616/1‑149 (MQ=255)
                                                                                                                                                  gCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAg  >  4:256757/1‑149 (MQ=255)
                                                                                                                                                  gCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAg  <  2:1211732/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGTAGATAACTTCTCTTGACGTGGTTTCAACTTCATACACTACGGTCGGCGCAGTGGTGATCAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTCAGGCTGAGTTTACCCAGCGCGTCACGGAAGGCTTCATAGTCGTCGGAACTTACCGGGAACAGACCGGCGTATACCTGCGGTTTGACTTTCTTAAAGCCAG  >  NC_000913/2705974‑2706268

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: