Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,799,370 Δ1 bp 9.0% intergenic (‑100/+23) pheM ← / ← rplT phenylalanyl‑tRNA synthetase operon leader peptide/50S ribosomal subunit protein L20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,799,3700A. 9.0% 597.5 / 50.5 178intergenic (‑100/+23)pheM/rplTphenylalanyl‑tRNA synthetase operon leader peptide/50S ribosomal subunit protein L20
Reads supporting (aligned to +/- strand):  ref base A (77/85);  new base . (10/6);  total (87/91)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

Only 100 of 179 total aligned reads displayed.

TTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAAT  >  NC_000913/1799225‑1799516
                                                                                                                                                 |                                                                                                                                                   
 tCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAg                                                                                                                                                 >  3:1285141/1‑149 (MQ=255)
       gCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCct                                                                                                                                          <  2:187318/149‑1 (MQ=255)
         tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCctct                                                                                                                                        <  3:221464/149‑1 (MQ=255)
         tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCctct                                                                                                                                        <  3:709995/149‑1 (MQ=255)
         tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCctct                                                                                                                                        <  1:741405/149‑1 (MQ=255)
           aaGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCttt                                                                                                                                      <  2:1321534/149‑1 (MQ=255)
           aaGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCCTCttt                                                                                                                                      <  2:87725/149‑1 (MQ=255)
            aGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCCTCtttt                                                                                                                                     <  2:308189/149‑2 (MQ=255)
              tAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCaa                                                                                                                                   >  1:242048/1‑149 (MQ=255)
               aaaaaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAAc                                                                                                                                  >  3:208297/1‑149 (MQ=255)
                aaaaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGCTAGTCTCCctctagatcgg                                                                                                                                 >  2:341743/1‑142 (MQ=255)
                     agaagCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCtt                                                                                                                            >  3:1126738/1‑149 (MQ=255)
                      caagCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTACTGTAGAGGCTACTGCGATAAATGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGTAGACTAGCTCCCTCTTTCAACTGGCTTa                                                                                                                           <  2:1057225/148‑1 (MQ=255)
                           ggAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCa                                                                                                                      <  3:1166907/149‑1 (MQ=255)
                            gAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCATAAGCCag                                                                                                                     <  1:991337/149‑1 (MQ=255)
                             aaaaTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCATAAGCCagt                                                                                                                    <  4:1143683/149‑2 (MQ=255)
                                  agcagcATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAgctgc                                                                                                               <  2:242048/149‑1 (MQ=255)
                                     agcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCttt                                                                                                            >  4:1440441/1‑149 (MQ=255)
                                      gcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTc                                                                                                           <  1:1281070/149‑1 (MQ=255)
                                      gcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTc                                                                                                           >  2:1124146/1‑149 (MQ=255)
                                       cattcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCg                                                                                                          <  1:297227/149‑1 (MQ=255)
                                        attcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGc                                                                                                         <  4:1148281/149‑1 (MQ=255)
                                           cattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                      <  1:1363891/149‑1 (MQ=255)
                                                cttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAGGGCTATGGCGATAAAAGTAAATGTTTTGATGGCGCTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                 >  3:1219557/1‑149 (MQ=255)
                                                cttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                 <  1:1376867/149‑1 (MQ=255)
                                                 ttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAAc                                                                                                <  1:653957/149‑1 (MQ=255)
                                                 ttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAAc                                                                                                <  4:660544/149‑1 (MQ=255)
                                                 ttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCGCCCTCTTTCAAc                                                                                                                                  <  1:328592/115‑1 (MQ=255)
                                                       tACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGc                                                                                          <  2:211058/149‑1 (MQ=255)
                                                        aCCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCt                                                                                                        >  2:949053/1‑134 (MQ=255)
                                                          cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGt                                                                                       <  3:1450920/149‑1 (MQ=255)
                                                          cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGt                                                                                       <  1:202601/149‑1 (MQ=255)
                                                          cTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGt                                                                                       <  3:803213/149‑1 (MQ=255)
                                                           ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTg                                                                                      >  3:1398271/1‑149 (MQ=255)
                                                                gTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAAc                                                                                                >  2:1180454/1‑134 (MQ=255)
                                                                gTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAAc                                                                                                <  1:1180454/134‑1 (MQ=255)
                                                                gTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGc                                                                                 <  4:733628/149‑1 (MQ=255)
                                                                 tACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCt                                                                                >  2:306200/1‑149 (MQ=255)
                                                                 tACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGAGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCt                                                                                <  2:521318/149‑1 (MQ=255)
                                                                  aCTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTa                                                                               <  4:1126738/149‑1 (MQ=255)
                                                                   ctctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTAc                                                                              >  3:293872/1‑149 (MQ=255)
                                                                    tctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACt                                                                             >  2:1103119/1‑149 (MQ=255)
                                                                    tctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGACAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACt                                                                             >  4:428791/1‑149 (MQ=255)
                                                                     ctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACtt                                                                            >  3:1416029/1‑149 (MQ=255)
                                                                      ttCAGAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACttt                                                                           >  3:964937/1‑149 (MQ=255)
                                                                       tCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCt                                                                                <  2:833590/143‑1 (MQ=255)
                                                                        cAAAAGACCTTTATTGAAAAGGCTACGGCGATTAAAGTCAATGTTTTGATGGCGTTGAAACGATAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                      <  4:12828/120‑1 (MQ=255)
                                                                        cAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                         >  1:934646/1‑149 (MQ=255)
                                                                         aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                        <  1:1103119/149‑1 (MQ=255)
                                                                         aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                        <  4:525318/149‑1 (MQ=255)
                                                                          aaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCg                                                                       >  1:628871/1‑149 (MQ=255)
                                                                          aaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCg                                                                       >  4:580398/1‑149 (MQ=255)
                                                                                cTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACt                                                                 >  4:501178/1‑149 (MQ=255)
                                                                                cTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACt                                                                 <  1:740846/149‑1 (MQ=255)
                                                                                cTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACt                                                                 >  2:465336/1‑149 (MQ=255)
                                                                                 tttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaataaa                                                                                             >  3:637550/1‑116 (MQ=255)
                                                                                 tttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTg                                                                >  1:372441/1‑149 (MQ=255)
                                                                                  ttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTTCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                               >  1:461544/1‑149 (MQ=255)
                                                                                  ttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                               <  3:531917/149‑1 (MQ=255)
                                                                                   tATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCg                                                              <  2:991960/149‑1 (MQ=255)
                                                                                    aTTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGa                                                             >  2:584298/1‑149 (MQ=255)
                                                                                     ttGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGat                                                            >  2:958264/1‑149 (MQ=255)
                                                                                        aaaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATc                                                         >  2:711380/1‑149 (MQ=255)
                                                                                        aaaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACCAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCg                                                                                         <  2:1208230/117‑1 (MQ=255)
                                                                                        aaaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACAAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCg                                                                                         >  1:1208230/1‑117 (MQ=255)
                                                                                         aaaGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCa                                                        <  2:158716/149‑1 (MQ=255)
                                                                                             gCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCa                                                    <  3:139215/149‑1 (MQ=255)
                                                                                                   gCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCt                                              >  1:972818/1‑149 (MQ=255)
                                                                                                     gATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCtca                                            >  2:1274401/1‑147 (MQ=255)
                                                                                                     gATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTa                                            <  4:209109/149‑1 (MQ=255)
                                                                                                     gATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTa                                            <  4:613263/149‑1 (MQ=255)
                                                                                                      aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                        >  1:1052422/1‑120 (MQ=255)
                                                                                                      aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTAc                                           >  3:224949/1‑149 (MQ=255)
                                                                                                          aaGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTTGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCTACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGCATCTTACGGTc                                       <  1:1207201/149‑1 (MQ=255)
                                                                                                             tCATTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                      <  3:537001/83‑1 (MQ=255)
                                                                                                             tCATTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                      >  4:537001/1‑83 (MQ=255)
                                                                                                               aaTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGAttt                                  <  3:414153/149‑1 (MQ=255)
                                                                                                                aTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTc                                 >  3:874149/1‑149 (MQ=255)
                                                                                                                aTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTc                                 >  4:1326498/1‑149 (MQ=255)
                                                                                                                  gTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCaa                               <  4:1093964/149‑1 (MQ=255)
                                                                                                                   ttttGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCGACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAAc                              >  2:421451/1‑149 (MQ=255)
                                                                                                                              ttGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCttttt                   >  2:1478328/1‑149 (MQ=255)
                                                                                                                                  aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                               <  3:1270316/101‑1 (MQ=255)
                                                                                                                                  aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                               >  4:1270316/1‑101 (MQ=255)
                                                                                                                                  aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg               <  4:293872/149‑1 (MQ=255)
                                                                                                                                  aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg               <  4:964937/149‑1 (MQ=255)
                                                                                                                                  aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg               <  3:1097847/149‑1 (MQ=255)
                                                                                                                                  aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg               >  2:641607/1‑149 (MQ=255)
                                                                                                                                  aaCGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg               <  1:1124146/149‑1 (MQ=255)
                                                                                                                                      aaaaGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCa           <  1:1218350/149‑1 (MQ=255)
                                                                                                                                       aaaGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAt          <  3:356116/149‑1 (MQ=255)
                                                                                                                                       aaaGAGGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAt          >  4:1445187/1‑149 (MQ=255)
                                                                                                                                        aaGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt         <  2:942534/149‑1 (MQ=255)
                                                                                                                                         agagGGAGACTAGCTCCCTCTTTCAACTGGCTTAAGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACCCTACTTTGTCGAATACTGCGATATCAGCCAGGAACTTACGGTCGTTTTCAACAGAGGCTTTTTTCAGGCCATtg        <  2:138880/149‑1 (MQ=255)
                                                                                                                                         agagGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg        >  2:813970/1‑149 (MQ=255)
                                                                                                                                         agagGGAGCTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg        >  1:1305280/1‑149 (MQ=255)
                                                                                                                                           agGTAGACTAGCGCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACGGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtgat      <  3:1025656/149‑1 (MQ=255)
                                                                                                                                           agGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtgat      <  4:329169/149‑1 (MQ=255)
                                                                                                                                           aaggAGACTAGCTCCCTATTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCGTtgat      >  4:417842/3‑149 (MQ=255)
                                                                                                                                               agaCTAGCTCCCTCTTTCAACTGGCTTATGCCAGTGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAAt  <  1:958264/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAGCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAAT  >  NC_000913/1799225‑1799516

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