Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065843,513,9180G.26.8% 101.3 / 39.2 42intergenic (‑94/+245)tldD/LY180_16725protease TldD/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (11/19);  new base . (4/7);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTT  >  CP006584/3513624‑3514200
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                          
acCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGttt                                                                                                                                                                                                                                                                                         <  1:271420/298‑1 (MQ=255)
        ccTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCt                                                                                                                                                                                                                                                                                 >  1:92783/1‑298 (MQ=255)
           gATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGAt                                                                                                                                                                                                                                                                              <  1:39047/298‑1 (MQ=255)
                catcgtttttgctactcgttagttactgcagtagaggattttttacggctgccggatgcggcgtgaacgccctatccagcctacggttatgttccggtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgtcggatgcggcgtaaacgccttatccgacctacggttatgttccgtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgccggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGAc                                                                                                                                                                                                                                                                         <  1:176087‑M1/36‑1 (MQ=255)
                      tAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGACGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCAGCTAGCAATTGTTCACTTACCAGGTTAAGCCTCATCGTCCTTGCTCCTCGTTAGTTGCTGCCATAGAGGAATTTTTACGGCTCCCGGCTGCGGCCCGCACGCCCCATCCACCCTCCGGTTTTGTTCCGCTTTGTGGGCCttaataacccgcccg                                                                                                                                                                                                                                                                   >  2:242618/1‑284 (MQ=255)
                       aaGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGc                                                                                                                                                                                                                                                                  <  1:352938/298‑1 (MQ=255)
                         cgtgcgattttAATAATGCAGTGTTCTACAACCCAGGTGTCGTGATAGCGCGACGGAAAAAGGAGATCGCCATAATCAAGGGGCCGTTCGGCCAGTTGACCGAGGTTCGTGAACAAGTCCTGATGTTTCAGGCCGTTCGCCACTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTGGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTGCCGGTTTGTAGGCCTGATAAGACGCACAGCGt                                                                                                                                                                                                                                                                <  2:62458/290‑1 (MQ=255)
                          gctactcgttagttactgcagtagaggattttttacggctgccggatgcggcgtgaacgccctatccagcctacggttatgttccggtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgtcggatgcggcgtaaacgccttatccgacctacggttatgttccgtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgccggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTc                                                                                                                                                                                                                                                               <  1:128175‑M1/46‑1 (MQ=255)
                                 tttAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAg                                                                                                                                                                                                                                                        >  1:358444/1‑298 (MQ=255)
                                    aataatGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCACTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTcagatcggaag                                                                                                                                                                                                                                                     >  1:62458/1‑288 (MQ=255)
                                      taatGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAAc                                                                                                                                                                                                                                                   <  1:117944/298‑1 (MQ=255)
                                          tgcagtagaggattttttacggctgccggatgcggcgtgaacgccctatccagcctacggttatgttccggtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgtcggatgcggcgtaaacgccttatccgacctacggttatgttccgtttgtaggcctgataagacgcacagcggcgcatcaggcaacggctgccggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAGAACGCGCAGCGTCGCATCCGGCAACGGCt                                                                                                                                                                                                                                               <  1:224825‑M1/62‑1 (MQ=255)
                                           cGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTg                                                                                                                                                                                                                                              <  1:261731/298‑1 (MQ=255)
                                            ggTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGt                                                                                                                                                                                                                                             <  1:364762/298‑1 (MQ=255)
                                                ttctAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGa                                                                                                                                                                                                                                         >  1:55526/1‑298 (MQ=255)
                                                tcctAAACCCCAGGATACGGGGTGGCAGACTGGAAAATAAGGATCGCATTAATCAAGGCGAGGTTCGGCCAGTTGACGGAGAAGAGCGAACAATGCCCGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACGGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGa                                                                                                                                                                                                                                         <  2:387440/296‑1 (MQ=37)
                                                                      tAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAAAGCCTTATCCg                                                                                                                                                                                                                   >  1:387154/1‑298 (MQ=255)
                                                                                ttaaaaTAGATGTCGCCAGAATCCAGGTCACGGTCGGCCAGTTGACCGAGGATTGCGAACAAGTCCTGATGTTTCCGGCCGTTCGCCGCTAGCCATTGTTCACTTACCAGGGTAAGACTCCTCGTTTTGGCTACTCGTTAGTTACTTCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAGCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGt                                                                                                                                                                                                          <  2:189969/296‑1 (MQ=37)
                                                                                   aaaTAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTACACGCCTTATCCGACCTACGGTTAt                                                                                                                                                                                                       >  1:179023/1‑298 (MQ=255)
                                                                                           accGCCATAATCAAGCGGCCGTTCGGCCAGTTGACCAGGGATCGCGAACAATTCCGGTTTTTTCAGGCCGTTCGCGGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGAGAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTCATGTTCCGtt                                                                                                                                                                                               <  2:142606/296‑1 (MQ=255)
                                                                                            tCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTCATGTTCCGttt                                                                                                                                                                                              >  1:142606/1‑298 (MQ=255)
                                                                                               ccATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTa                                                                                                                                                                                           <  1:374827/298‑1 (MQ=255)
                                                                                                 aTAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGCGAACCCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGCTCACACGCCCAGCGTCCCCTCCGGCAACGGCTGTCGGATGCGGCGGAAACGCCCTACCCCGCCTACGTCTATGTTCCGTTTGtctg                                                                                                                                                                                         >  2:403235/1‑295 (MQ=255)
                                                                                                                 gcggcacattcaCCGGGGATCAGGAACTAGGCCTGAGTGTTGGGGTCGTTGGCCGGTGGTAATGGTTGACTTCCCATGTTAGGACTATTCTTTGTGGCTACTGGTTGGTTACTGCAGGGGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTCTCAAGCCGACGGTGATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                                                                                                                         <  2:220939/287‑1 (MQ=18)
                                                                                                                           gACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                                                                                                                               <  1:220092/298‑1 (MQ=255)
                                                                                                                                              agattttttacggctgccggaggcggcgtgagcagccaaaccatcctacgtttatgttccggtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgtcggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgc                                                                                                                                           >  1:37941‑M1/137‑298 (MQ=255)
                                                                                                                                                      tacggctgccggatgcggcgtgaacgccctatccagcctacggttatgttccggtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgtcggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGaa                                                                                                                                   >  1:296567‑M1/129‑298 (MQ=255)
                                                                                                                                                       tttCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTACGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGAACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGGaaa                                                                                                                                   <  1:52677‑M2/298‑4 (MQ=255)
                                                                                                                                                          gctgccggatgcggcgtgaacgccctatccagcctacggttatgttccggtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgtcggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGcc                                                                                                                               <  1:37378‑M1/174‑1 (MQ=255)
                                                                                                                                                                                     aCTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTaaacgccttatccgacctacggttatgttccgt                                                                                                     <  1:73705‑M2/298‑34 (MQ=255)
                                                                                                                                                                                      cTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTaaacgccttatccgacctacggttatgttccgtt                                                                                                    >  1:182442‑M2/1‑264 (MQ=255)
                                                                                                                                                                                                gttatgttccggtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgtcggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCtt                                                                                         >  1:331920‑M1/87‑298 (MQ=255)
                                                                                                                                                                                                      ttccggtttgtaggcctgataagacgcacagcgtcgcatcaggcaacggctgtcggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAAGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGGTATGTTCTGGCTCGTTGGCCTTTGGCaa                                                                                   <  1:331998‑M1/218‑1 (MQ=255)
                                                                                                                                                                                                       cATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTaaacgccttatccgacctacggttatgttccgtttgtaggcctgataagac                                                                                   <  1:176087‑M2/298‑52 (MQ=255)
                                                                                                                                                                                                                 gCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTaaacgccttatccgacctacggttatgttccgtttgtaggcctgataagacgcacagcgtc                                                                         <  1:128175‑M2/298‑62 (MQ=255)
                                                                                                                                                                                                                  ggcctgataagacgcacagcgtcgcatcaggcaacggctgtcggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTAt                                                                       <  1:151295‑M1/230‑1 (MQ=255)
                                                                                                                                                                                                                                 tGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGGCGCATCAGGCAACGGCTGCCGGATGCGGCGTaaacgccttatccgacctacggttatgttccgtttgtaggcctgatagaacgcgcagcgtcgcatccggcaacggct                                                         <  1:224825‑M2/298‑78 (MQ=255)
                                                                                                                                                                                                                                                 caacggctgtcggatgcggcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCaaa                                        >  1:265853‑M1/38‑298 (MQ=255)
                                                                                                                                                                                                                                                                   gcgtaaacgccttatccgaCCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCttt                      <  1:98171‑M1/279‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      ccTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACt   >  1:304919/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                       cTACGGTTATGTTCCGTTTGTAGGCCTGATAAGCCGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTGTGTTGTGGCCCGTGGGCCTTTGGCAACGATTATCCTAGGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCaaagcgt  >  1:212998/1‑293 (MQ=37)
                                                                                                                                                                                                                                                                                       cTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACtt  <  1:403235/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                          
ACCAACGCCCTGATCGATGTTGTAAGAACCATCTTTAATAATGCGGTCTTCTAAAACCCAGGATTCGTGATAGCTCGACTGAAAATAGAGATCGCCATAATCAAGGCGACGTTCGGCCAGTTGACCGAGGATCGCGAACAAGTCCTGATGTTTCAGGCCGTTCGCCGCTAGCAATTGTTCACTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTCTGGCTCGTTGGCCTTTGGCAACGATTATCCTATGAGATTGGGGCAATTACGCGCCCTCGTCAAATCATTGCGCTTTTTCTTTACGCGGTTGGCGCAACACTT  >  CP006584/3513624‑3514200

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