Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,799,370 Δ1 bp 5.7% 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. 5.7% 673.7 / 30.4 193intergenic (‑100/+23)pheM/rplTphenylalanyl‑tRNA synthetase operon leader peptide/50S ribosomal subunit protein L20
Reads supporting (aligned to +/- strand):  ref base A (97/85);  new base . (6/5);  total (103/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.67e-01

Only 100 of 197 total aligned reads displayed.

GGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTT  >  NC_000913/1799222‑1799518
                                                                                                                                                    |                                                                                                                                                     
ggATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGgaga                                                                                                                                                       >  1:470255/1‑149 (MQ=255)
   ttCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTa                                                                                                                                                    <  3:355690/149‑1 (MQ=255)
      aGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑Ct                                                                                                                                                >  1:681735/1‑149 (MQ=255)
            tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCctct                                                                                                                                          <  1:1346034/149‑1 (MQ=255)
            tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCctct                                                                                                                                          >  1:1457747/1‑149 (MQ=255)
            tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCctct                                                                                                                                          >  2:10269/1‑149 (MQ=255)
            tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCctct                                                                                                                                          <  1:1243983/149‑1 (MQ=255)
            tAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCctct                                                                                                                                          <  3:625731/149‑1 (MQ=255)
             aaaGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCtt                                                                                                                                         <  2:470255/149‑1 (MQ=255)
              aaGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCttt                                                                                                                                        <  3:454852/149‑1 (MQ=255)
                gTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCa                                                                                                                                      >  2:85385/1‑149 (MQ=255)
                 tAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCaa                                                                                                                                     >  4:1428657/1‑149 (MQ=255)
                 tAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCaa                                                                                                                                     >  1:612900/1‑149 (MQ=255)
                 tAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCaa                                                                                                                                     <  1:257317/149‑1 (MQ=255)
                     aaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTgg                                                                                                                                 >  4:646581/1‑149 (MQ=255)
                     aaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTgg                                                                                                                                 <  1:10269/149‑1 (MQ=255)
                     aaaaGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTgg                                                                                                                                 <  4:471641/149‑1 (MQ=255)
                         gaagCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTa                                                                                                                             >  4:329572/1‑149 (MQ=255)
                          aagCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTAt                                                                                                                            <  1:695710/149‑1 (MQ=255)
                                 aaaTAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCagag                                                                                                                     <  3:1029819/149‑1 (MQ=255)
                                       cagcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCtt                                                                                                               <  1:221969/149‑1 (MQ=255)
                                       cagcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCtt                                                                                                               >  3:1536288/1‑149 (MQ=255)
                                        agcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAg                                                                                                                                                   <  4:796354/113‑1 (MQ=255)
                                        agcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAg                                                                                                                                                   >  3:796354/1‑113 (MQ=255)
                                         gcaTTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTc                                                                                                             <  4:732985/149‑1 (MQ=255)
                                           attcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGc                                                                                                           <  2:934780/149‑1 (MQ=255)
                                            ttcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAAc                                                                                                                                    >  4:363532/1‑123 (MQ=255)
                                            ttcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCt                                                                                                          >  1:503822/1‑149 (MQ=255)
                                            ttcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCtt                                                                                                                                         <  2:1186721/118‑1 (MQ=255)
                                            ttcattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAG‑CTAGTCTCCCTCtt                                                                                                                                         >  1:1186721/1‑118 (MQ=255)
                                              cattGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCttt                                                                                                        <  4:20697/149‑1 (MQ=255)
                                               attGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCtttt                                                                                                       >  1:886977/1‑149 (MQ=255)
                                                ttgcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCttt                                                                                                              >  3:253278/1‑141 (MQ=255)
                                                 tgcttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTc                                                                                                     <  2:294689/149‑1 (MQ=255)
                                                   cttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                   >  3:400781/1‑149 (MQ=255)
                                                   cttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCaa                                                                                                   <  4:267039/149‑1 (MQ=255)
                                                    ttgcGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCctct                                                                                                                                          >  2:1196631/1‑109 (MQ=255)
                                                         ttACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCagag                                                                                             >  1:1287429/1‑149 (MQ=255)
                                                          tACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGc                                                                                            <  2:544029/149‑1 (MQ=255)
                                                              ttttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTg                                                                                        >  1:1269406/1‑149 (MQ=255)
                                                                ttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGaa                                                                                      <  3:727298/149‑1 (MQ=255)
                                                                ttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGaa                                                                                      <  2:886977/149‑1 (MQ=255)
                                                                ttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGaa                                                                                      <  2:681735/149‑1 (MQ=255)
                                                                ttGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTCATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGaa                                                                                      <  4:908376/149‑1 (MQ=255)
                                                                   gTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGc                                                                                   >  3:981982/1‑149 (MQ=255)
                                                                   gTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGc                                                                                   >  2:1466772/1‑149 (MQ=255)
                                                                   gTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGc                                                                                   >  2:1452355/1‑149 (MQ=255)
                                                                    tACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAT‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCt                                                                                  >  3:704239/1‑149 (MQ=255)
                                                                     aCTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGTCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTa                                                                                 >  4:517252/1‑149 (MQ=255)
                                                                     aCTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTa                                                                                 >  4:1001838/1‑149 (MQ=255)
                                                                       tctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACt                                                                               >  4:121828/1‑149 (MQ=255)
                                                                       tctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACt                                                                               >  1:98426/1‑149 (MQ=255)
                                                                       tctTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACt                                                                               >  2:1407023/1‑149 (MQ=255)
                                                                         ttCAAAAGACCTTTATTGAGAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGaat                                                                                     <  4:1354611/141‑2 (MQ=255)
                                                                         ttCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACttt                                                                             >  2:781924/1‑149 (MQ=255)
                                                                         ttCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACAGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAAc                                                                                     >  3:1354611/1‑141 (MQ=255)
                                                                           cAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGt                                                                           <  4:400781/149‑1 (MQ=255)
                                                                            aaaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTc                                                                          <  3:517252/149‑1 (MQ=255)
                                                                             aaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCg                                                                         >  3:50173/1‑149 (MQ=255)
                                                                             aaaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCg                                                                         >  1:1048808/1‑149 (MQ=255)
                                                                              aaGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGa                                                                        >  1:555569/1‑149 (MQ=255)
                                                                                 aCCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATa                                                                     <  4:407241/149‑1 (MQ=255)
                                                                                   cTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCaga                                                                                                                      >  1:1197086/1‑98 (MQ=255)
                                                                                   cTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGc                                                                                                           <  1:1484437/109‑1 (MQ=255)
                                                                                   cTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACt                                                                   >  3:85442/1‑149 (MQ=255)
                                                                                     ttATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                                 <  2:1232653/149‑1 (MQ=255)
                                                                                       aTTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGa                                                               <  2:503822/149‑1 (MQ=255)
                                                                                        ttGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGat                                                              >  4:1458830/1‑149 (MQ=255)
                                                                                             aaGGCTACGGCGATAAAAGTCAATGATTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAg                                                         >  2:726060/1‑149 (MQ=255)
                                                                                                gCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCa                                                      <  2:612900/149‑1 (MQ=255)
                                                                                                 cTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTGCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAg                                                     <  3:1484052/149‑1 (MQ=255)
                                                                                                   aCGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGa                                                   <  4:705001/149‑1 (MQ=255)
                                                                                                    cGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGAt                                                  <  1:1407023/149‑1 (MQ=255)
                                                                                                         aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTAc                                             <  3:121828/149‑1 (MQ=255)
                                                                                                         aTAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTAc                                             <  3:371371/149‑1 (MQ=255)
                                                                                                            aaaGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGt                                          >  3:205411/1‑149 (MQ=255)
                                                                                                             aaGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTc                                         <  1:1023173/149‑1 (MQ=255)
                                                                                                               gTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGa                                       >  3:423701/1‑149 (MQ=255)
                                                                                                                tCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGc                                                                 <  4:466816/122‑1 (MQ=255)
                                                                                                                 cAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGAtt                                     <  3:557775/149‑1 (MQ=255)
                                                                                                                  aaTGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGAttt                                    >  4:631970/1‑149 (MQ=255)
                                                                                                                     gTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCaa                                 <  4:954455/149‑1 (MQ=255)
                                                                                                                      ttttGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAg                                                         >  4:962147/1‑124 (MQ=255)
                                                                                                                      ttttGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAg                                                         <  3:962147/124‑1 (MQ=255)
                                                                                                                      ttttGATGGCGTTGAAAAGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTgg                                                                                                                                 >  2:205513/1‑52 (MQ=255)
                                                                                                                      ttttGATGGCGTTGAAAAGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTgg                                                                                                                                 <  1:205513/52‑1 (MQ=255)
                                                                                                                         tGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACaga                             <  4:844477/149‑1 (MQ=255)
                                                                                                                          gATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACagag                            <  2:1483774/149‑1 (MQ=255)
                                                                                                                           aTGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAgg                           >  2:44748/1‑149 (MQ=255)
                                                                                                                              gCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCtt                        >  1:860581/1‑149 (MQ=255)
                                                                                                                                     aaCGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGCATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAg                 <  2:1367176/149‑1 (MQ=255)
                                                                                                                                      aCGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAgg                >  3:1000179/1‑149 (MQ=255)
                                                                                                                                      aCGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAgg                <  1:1466772/149‑1 (MQ=255)
                                                                                                                                       cGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGc               <  2:1287429/149‑1 (MQ=255)
                                                                                                                                       cGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGc               <  4:758840/149‑1 (MQ=255)
                                                                                                                                           aaGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt           <  4:1239178/149‑1 (MQ=255)
                                                                                                                                           aaGAGGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCAtt           >  3:839825/1‑149 (MQ=255)
                                                                                                                                            agagGGAG‑CTAGTCTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtg          >  4:1410977/1‑149 (MQ=255)
                                                                                                                                              agGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtgat        <  3:1428657/149‑1 (MQ=255)
                                                                                                                                               gggAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATtgatg       <  4:555084/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAGGCTACGGCGATAAAAGTCAATGTTTTGATGGCGTTGAAACGAAAAGAGGGAGACTAG‑CTCCCTCTTTCAACTGGCTTATGCCAGAGCTGCTTTCGCTTTTTCAACCAGAGCGGTGAACGCTACTTTGTCGAATACTGCGATATCAGCCAGGATCTTACGGTCGATTTCAACAGAGGCTTTTTTCAGGCCATTGATGAATTT  >  NC_000913/1799222‑1799518

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: ‑