Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913249,5380TG33.8% 106.7 / 11.7 78pseudogene (533/1713 nt)lfhApseudogene, flagellar system protein, promoterless fragment; flagellar biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (18/33);  new base G (23/3);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 8.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGTTG  >  NC_000913/249306‑249774
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
gAGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCa                                                                                                                                                                                                                                             <  1:92299/239‑1 (MQ=255)
   ggggTTACCCTTTCCTTTTACCCCGCGCGCCCCCGTTTTGTAACCGAGTTTAAACCTGAGAGGCTGGATCGGCGGAATGTTTTCCCAACTGCCGCGCTGCTCGTTGGTTTCGGGGATGGTAAGGTTAACCTTTTCGAGGCTTTTCGCTTCCTCCACAGCCGCCGGCTGTTCGTCGCCACCCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCAtt                                                                                                                                                                                                                                          <  2:1570807/237‑1 (MQ=255)
   cGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCAtt                                                                                                                                                                                                                                          <  1:2077699/238‑1 (MQ=255)
       ttACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCg                                                                                                                                                                                                                                      >  1:1793303/1‑238 (MQ=255)
        tACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGGATTTCCCgg                                                                                                                                                                                                                                     >  1:1421586/1‑239 (MQ=255)
               tgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCAtt                                                                                                                                                                                                                                          <  1:348736/226‑1 (MQ=255)
               tgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCAtt                                                                                                                                                                                                                                          >  2:348736/1‑226 (MQ=255)
               tgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCCGTaccacc                                                                                                                                                                                                                              >  1:1776774/1‑238 (MQ=255)
               tgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGGCTTCCCCGGaccacc                                                                                                                                                                                                                              >  1:1380944/1‑238 (MQ=255)
               tgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGGATTCCCGGGaccacc                                                                                                                                                                                                                              >  2:1112782/1‑238 (MQ=255)
               tgCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTaccccc                                                                                                                                                                                                                              >  1:419687/1‑238 (MQ=255)
                gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCg                                                                                                                                                                                                                                                >  2:1311599/1‑220 (MQ=255)
                gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCg                                                                                                                                                                                                                                                <  1:1311599/220‑1 (MQ=255)
                gCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGGaccaccc                                                                                                                                                                                                                             >  1:478248/1‑237 (MQ=255)
                    tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCCAGCAGGGCGCTGAACAGCAAAAAAGGCAGGGGGGGGATTTCCCCGTacacccccca                                                                                                                                                                                                                          >  2:27001/1‑230 (MQ=255)
                    tGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCAccgcccgc                                                                                                                                                                                                                        >  1:1486010/1‑239 (MQ=255)
                    cgcccaccgGCGCCACCTTTTGGTACCCGAGGCTTAATCTGATGGGCCCGATCGGCGAATTGTTTCCCCAACGGCCCTGCTCCTCGTTGTTTTGGGTGATGGGGGGGTTCGGCTTTTCGAGGCTTTTCTCTCCCCCCCCCGCCGCCAGCGGCTGTTTCCTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGc                                                                                                                                                                                                                        <  2:2032150/236‑1 (MQ=255)
                               ccccccaTTTTGTACCGGAGGTTTAACCTGTTAGGCTGGATCGGGGAATTGTTTTCCCACCTGCCCTGCTGCTCGTTGGTTGGGGTGAGGGTGGGGGTCAGCGTTTCGGGGTTTTTCTCTCCCGCCCCCGCCGCCGGCGGCTTTTTGCTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACAt                                                                                                                                                                                                               <  2:2210766/232‑1 (MQ=255)
                                   ggtgTTGGTAACGGGGTCATAAACTTATGGGCGCGATCAGGGGAGTTGTTCCCAAGCTGACCTGCTGCTCGTGGTTTGGGGTGGTGGGGCGGGTAAGCTTTTCGAGGCTTTCCTCTCCCCCCCCCGCCGCCGGGGGTTTTTTCCTCACCCCCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATccc                                                                                                                                                                                                         <  2:2315927/235‑1 (MQ=255)
                                       ttttTACCCGAGGCTTAAACTGTTGGGCTCGATCAGCGTAATGGTTTCCCAACTGCCCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCAtt                                                                                                                                                                                                                                          <  2:19203/202‑1 (MQ=255)
                                                 aGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCACTCCCCGGACCACCGCCCGCACGAACATAAACCCGGGGGCGGTGGagg                                                                                                                                                                                           >  1:1714588/1‑237 (MQ=255)
                                                 aGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGCACATAATCCCGGTGGCgggg                                                                                                                                                                                               >  1:2315927/1‑232 (MQ=255)
                                                 aGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGGacca                                                                                                                                                                                                                                >  1:19203/1‑202 (MQ=255)
                                                             aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCCCCCGCCCGCACATAATCCCCGTGGCGGGGTAGAGCCCCCCCGGGCt                                                                                                                                                                              >  1:374355/1‑239 (MQ=255)
                                                             aTAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCACCACGAACATAATCCCGGGGGCGGGGTAGAGCACCGACGGGCt                                                                                                                                                                              >  1:1218204/1‑239 (MQ=255)
                                                                           cGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCCATCCCGGTACCCCCCGCAGCACGAACATAAACCCCGGGGGGGGGTAAGGCACCGCCGGGCTTGccctgcgctgcg                                                                                                                                                                >  2:167067/1‑229 (MQ=255)
                                                                                     tCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCggga                                                                                                                                                                                                                                             >  1:507130/1‑152 (MQ=255)
                                                                                     tCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGGAATCCCGGGACCCCCCCCCGCCCGCAAAAAAATCCCGGTGGGGGGGAGAGCGCCCGCCGGGTTGCCGGCACCCGGTGTAAGCACg                                                                                                                                                      >  2:597629/1‑238 (MQ=255)
                                                                                              aCCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCCCGAACATAATCCCGGTGGCGGGGTAGAGCACCGACGGGCTTGCCCGCAGCTGGTTACGCACGTCATGGGCGat                                                                                                                                           >  1:1643227/1‑239 (MQ=255)
                                                                                              aCCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCggt                                                                                                                                           >  1:1804299/1‑237 (MQ=255)
                                                                                              aCCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGACGGGCATTCCCTGTACCCCCCCCAGCACCAACATAAACCCCGCGGCGGTGTAGAGAACCGACGGGCTTGCCAGCAGCCTGGTACGCAAGGCTTAGGCGag                                                                                                                                           >  2:429666/1‑238 (MQ=255)
                                                                                                   gtgcctggtGGTTTGCGGTGAGGGGGCGGCTCGGGTTTTGAGGGTTTTCCCCTCCCCCCCCCGCCGCCGGCGCCTTTTTGCTCACCCGCCAGCCGGTAAACCCAGGCAGGGCGCTGAACAGAAAAAACGGCGGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCg                                                                                                                                       <  2:1776774/231‑1 (MQ=255)
                                                                                                      ctCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAAAAGCAAAAACGGCAGGGGCGGCATTCCCCGTACCACCGGCCACACGAACAAAACACCGGGGGCGGGGGAGGACCACGACGGGCGTGGCCGCCGATAGTTTCCCACGTCCTCGGCCGTATCGCCGt                                                                                                                                   >  2:2378622/1‑239 (MQ=255)
                                                                                                                     tGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCCCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCGCCGCCAGCACGAACATAATCCCgg                                                                                                                                                                                                       <  1:773914/159‑1 (MQ=255)
                                                                                                                      gATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGGACCCCCGCCAGCACGAAAATAATCCCGGTGGCGGTGTAGAGCACCCACGGGCTTTCCAGCAGCGTGGTTCGCACGGCCTGGGCGGTATTCCCCGTGTCGCTTAccccgg                                                                                                                    >  2:827402/1‑234 (MQ=255)
                                                                                                                      gATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCCCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCgg                                                                                                                     <  1:2176628/238‑1 (MQ=255)
                                                                                                                         ggggggggTCAGCTTTTCGAGGCTTTTCTCTCCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCcagcag                                                                                                                                                                     <  1:1683579/184‑1 (MQ=255)
                                                                                                                               ggTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGGATTCCCGGGACCACCGCCAGCACGAACATAAACCCgggggggggggt                                                                                                                                                                                             >  2:773914/1‑155 (MQ=255)
                                                                                                                                gTCAGCGTTTCGAGGCTTTTCTCTCCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTg                                                                                                                                                                            <  1:230776/175‑1 (MQ=255)
                                                                                                                                 tCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATc                                                                                                                                                                                                           >  1:1292478/1‑143 (MQ=255)
                                                                                                                                 tCAGCCTTTCCGGGCTTTTTCTCTTCCCCCTCCGCCGCCGGGGGCTGTTTGCTCACCCCCCAGCCGGTAAACCCAAGCAGGGCGCTGAAAAGAAAAACAGGCAGGTGAGGGTTACCCGGTACCACCGCCAGCACGAACATAATCCCGGTgg                                                                                                                                                                                                    <  2:279423/151‑1 (MQ=255)
                                                                                                                                   agagTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGGGGCAATACCGGGAACACCGCCAGAACGAACAAAAAACAGGTGGCGGTGAAGAGACCCGACGGGCGTGCCAGCAGCTGGGTCCGCACGACAGAGCGGAAATGCCCGGTGGCGTGGACGCGGGCGACGATAATcg                                                                                                       >  2:753275/4‑238 (MQ=255)
                                                                                                                                   aGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCCGTACCACCGCCCGCACGAACATAAACCCCGTGGGGGGGGAGAGCACCGACGGGCTTGCCCGCCGCTGGTTACGc                                                                                                                                                         >  2:1683579/1‑189 (MQ=255)
                                                                                                                                     cGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCg                                                                                                                                                                                                                                      <  2:2230744/113‑1 (MQ=255)
                                                                                                                                     cGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCg                                                                                                                                                                                                                                      >  1:2230744/1‑113 (MQ=255)
                                                                                                                                        tcgagggCTTTCCTCTCCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAACCCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATcgccgccg                                                                                                 <  1:1933366/234‑1 (MQ=255)
                                                                                                                                         tccgagttttttCTCCCCCCCCCCCGCCCCCGGCGGTTTTTGCCTCACCCGCCAGCCGGTAAAACCAAGCGGGGCGCTGAACAGCAAAAACGGCGGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATcgccgccg                                                                                                 <  2:419687/231‑1 (MQ=255)
                                                                                                                                          cGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGGTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCcgcg                                                                                               >  1:1085377/1‑239 (MQ=255)
                                                                                                                                          cGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCCTTCCCGGTACCACCCCCAGCACGAACATAATCCCGGGGGCGGGGTAGAGCACCGACGGGCTTGCCCGCAGCTg                                                                                                                                                                 >  2:230776/1‑175 (MQ=255)
                                                                                                                                          cGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCCGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGTGTagag                                                                                                                                                                                          >  1:279423/1‑151 (MQ=255)
                                                                                                                                                  tttgctagcccgcctccgccCCCGGGGGCTGTCTGCACACACGCGCGCCGGTAAACCCAAGCAGGGCGCGTAACAGAAAAAACGGCAGGGGCGGCATCCCCGGTACCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGgag                                                                                         <  2:880918/230‑1 (MQ=255)
                                                                                                                                                   ttCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCCGCACGAACATAATCCCGGTGGCGGGGTAGAGCACCGCCGGGCTTGCCAGCAGCTGGTTACGCACGGCATTGGCGATATCGCCGTTGTCGCTGACGCGGGGGAGGCTAAACGCCCCCCCGgggggggg                                                                                       >  2:478926/1‑235 (MQ=255)
                                                                                                                                                   ttCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAAAATAATCCCCGGGGCGGGGTAGAGCCCCGACCGGCTTGCCAGCAGCCGGTTACGCCCGTCATGGGCGGTATCGCCGTTTTGGCTGACCCGGGGGGCGGAAAACGCCGCCCCCGTggggg                                                                                        >  2:1037928/1‑235 (MQ=255)
                                                                                                                                                              cctccgccGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCCGGACCACCGCCAGCACGAACATAATCCCGGGGGCGGGGTAGAGCACCGCCGGGCTTGCCAGCAGCTGGTTACGCACGGCACGGGGGATATTCGCGTGGTTGCGTAAGCGGGGGCAGGACAACGCCGCCGCGGGGGTGGGCAGCcgggga                                                                             >  2:655865/1‑237 (MQ=255)
                                                                                                                                                                 ccgccgCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCCCCAGCACGAACATAATTCCGGTGGCGGGGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCa                                                                                                                                                        >  1:1690683/1‑162 (MQ=255)
                                                                                                                                                                                ttttGCTCACCCGCCAGCCGGTAAAACCAAGCGGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAgg                                                           <  1:1190307/238‑1 (MQ=255)
                                                                                                                                                                                     cTCACCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATcgccgccg                                                                                                 <  1:1478385/196‑1 (MQ=255)
                                                                                                                                                                                         tCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCa                                                                                                                                                                                                                       >  2:930310/1‑78 (MQ=255)
                                                                                                                                                                                         tCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCa                                                                                                                                                                                                                       <  1:930310/78‑1 (MQ=255)
                                                                                                                                                                                          ccGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGGGTAGAGCACCGACGGGCTTGCCCGCAGCTGGTTACGCACGTCCTGGGCGATATCGCCGTTGTCGCTGACGCGGGGGGCGATAATCGCCGCCGCGGGGGAGAGCAGGAGGGGAGGGGTCTGCGGCAACAGGccgccgccgg                                                 >  2:64955/1‑238 (MQ=255)
                                                                                                                                                                                            gCCAGGCGGTAAAACCAGGCAGGGCGCGAAACCACAAAACCGGCAGGTGCGGCATTCCCGGACCCACCGCAAGCACGAACATAATCCCGGTTGCGGTGTAGAGCAACGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATAGCGCCGTTGTCGCTGACGCGGGTGACGATAATCCCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATg                                               <  2:684311/239‑1 (MQ=255)
                                                                                                                                                                                            cccgcgcGGTAAAACCAAGCAGGGCGCGGAACAGCAAAAACGGCAGGTGCGGCATCCCCGGCCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATg                                               <  2:971410/233‑1 (MQ=255)
                                                                                                                                                                                               aGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCAGCACGAAAATAATCCCGGTGGCGGTGTAGAGCACCGCCGGGGTTGCCAGCAGCTGGTTACGCACGGCATGGGCGAGATCGCCGTTGTCGCTGGCGCGGGGGACGATAATCGCCGCCGCGGTGGggag                                                                                       >  2:1478385/1‑196 (MQ=255)
                                                                                                                                                                                                 ccGGTAAACCCAAGCAGGGCGCTGAACAGCAAAACCGGCGGGTGCGGCATCCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtca                                           <  1:429666/238‑1 (MQ=255)
                                                                                                                                                                                                  gggaaaacCCAGGCAGGCCGCTAACAAGCAAAACCGCCAGGGGCGCCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtcat                                          <  2:1380944/234‑1 (MQ=255)
                                                                                                                                                                                                  cGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGtca                                           <  1:1519907/237‑1 (MQ=255)
                                                                                                                                                                                                      aaaaCCAAGCAGGGCGCTGAACAGCAAAAACGGCGGGTGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGc                                      <  1:1037928/238‑1 (MQ=255)
                                                                                                                                                                                                        acccAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATcgccg                                                                                                    <  1:2313710/172‑1 (MQ=255)
                                                                                                                                                                                                        aaCCAAGCAGGGCGCTAAACAGCAAAACCGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGcaca                                   <  1:1488013/239‑1 (MQ=255)
                                                                                                                                                                                                                  ggCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGatat                                                                                                                                         <  2:1551928/128‑1 (MQ=255)
                                                                                                                                                                                                                  ggCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGatat                                                                                                                                         >  1:1551928/1‑128 (MQ=255)
                                                                                                                                                                                                                  ggCGCTGAACAGCAAAAACGGCAGGGGCGGCATTCCCGGTACCACCGCCCGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGGCATGGGCGATATCGCCGTTGTCGCTGACGCGGGGGACGAAAATAGCCGCCGCGggggg                                                                                          >  2:2313710/1‑173 (MQ=255)
                                                                                                                                                                                                                   gcgcTAAACAGCAAAAACGGCGGGGGCGGCATTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATActgctg                           <  1:167067/236‑1 (MQ=255)
                                                                                                                                                                                                                               aaaaaGGACAGGTGCGGCTTCCACGATCCAGCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGGGCAGCAGGGAAGGGATGTGCGCCGCCAGGCCGGCGCCGATGGGCATCAGCACATACGGCTGGAAGGCAGCAt                <  2:2076665/235‑1 (MQ=255)
                                                                                                                                                                                                                                 aaaCGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGGGGCGGGGTAGAGCACCGACGGGCTTGCCAGCCGCTGGTTACGCCCGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGGGACGATAATCGCCGCCGCGGGGGGGAGCAGCAGGGAAGGGGTCTGCGCCCCCAGGCCGTGGCCGGTGGTCATCAGCACATACTGCTGGAAGGGCGCATCgcgg           >  2:2393933/1‑235 (MQ=255)
                                                                                                                                                                                                                                        ggttGCGGCTTTCCCGGTCCCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGt    <  1:753275/235‑1 (MQ=255)
                                                                                                                                                                                                                                        gggTGCGGCATCCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGt    <  1:2393933/237‑1 (MQ=255)
                                                                                                                                                                                                                                        gggTGCGGCACTCCCGGTCCCCCCGCCAGCACGAACATAGTCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCGGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCCCCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGt    <  2:378281/237‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTGCGGCCTTCCCGGTCCCCCCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGtt   <  1:1112782/238‑1 (MQ=255)
                                                                                                                                                                                                                                           tGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGTTg  <  1:720700/237‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                 
GAGCGGGTTACCCTGTGCTTTGTCCACCAGCGCCACCAGTTTGTAACCGAGGCTTAAACTGATAGGCTCGATCAGCGGAATGGTTTCCCAACTGACCTGCTGCTCGCTGGTTTCGGTGATGGTGCGGGTCAGCGTTTCGAGGCTTTTCTCTTCCGCCTCCGCCGCCAGCGGCTGTTTGCTCATCCGCCAGCCGGTAAAACCAAGCAGGGCGCTGAACAGCAAAAACGGCAGGTGCGGCATTCCCGGTACCACCGCCAGCACGAACATAATCCCGGTGGCGGTGTAGAGCACCGACGGGCTTGCCAGCAGCTGGTTACGCACGTCATGGGCGATATCGCCGTTGTCGCTGACGCGGGTGACGATAATCGCCGCCGCGGTGGAGAGCAGCAGGGAAGGGATCTGCGCCACCAGGCCGTCGCCGATGGTCATCAGCACATACTGCTGGAAGGCAGCATCGGCGCTCAGGTTG  >  NC_000913/249306‑249774

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