Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,2200GC28.9% 70.1 / 18.2 41intergenic (‑363/+132)lysO/aqpZL‑lysine/thialysine exporter/water channel AqpZ
Reads supporting (aligned to +/- strand):  ref base G (24/3);  new base C (0/11);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 3.02e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTT‑GCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAA‑TCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAAT  >  NC_000913/914921‑915516
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                        
gACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTT‑GCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAA‑TCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgccg                                                                                                                                                                                                                                                                                                          >  1:303094/1‑300 (MQ=255)
                                  gCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTT‑GCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAA‑TCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACANGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTANGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATg                                                                                                                                                                                                                                                                        >  2:588991/1‑300 (MQ=255)
                                              gTAGAGTATTCGTTTTCATTTAAAGATATTCTT‑GCGCTTTAATCACAAACTGCACCGATGTTGGTAGCGTCAAAA‑TCGACGAGGCGTTCCCTGAAGGCCGGGGCAGTCCACATGGATGTAGGCTNAGGGCGCGTTTTACAGGGATGTTACCNCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGGGGTGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                                                                                                                                                                                            >  2:501859/1‑300 (MQ=18)
                                                     gttttttttCATGGAAAGAT‑TTCTTTGCGCTTAAATAAAAAACTGAACCGTTTTGGGTGGCGTAAAAC‑TCGCCGAGTTTTTCCCGGAGGGCTGGGGCACCCCACATGGGTGGGGGTGGGGGGCGTTTTTGAAGGGGGTGTTACCTCACGCCCGACCAGGGAGCTGGAGGGGAGAAGT‑AAGGGACACCACGCAGTGGCGGTTTTGTTCGCAAGGGCCTGGGGGGGGGGGGGGTGGCGGCGGTTGGCCCCCCCCTACGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                                                                                                                                                                                     <  1:706577/295‑1 (MQ=255)
                                                     cttttttttCATTGAAAGATATTCTT‑GTGCTTTACTTACAAACTGCACCGATGGTGG‑GGACTCAAAAATCACTGATGCGTCCCCGGAGGGCTGGGACAGCCCACATGGATGTGGGCGGAGGGCGCTTTTTACGGGGATGTCACCTCGCGCCCGACCTGTTATCCGGAAGGGATAAGTCGAGGG‑CACCGCACAGCGGCGGTTTTTTCCGCCAGAGCCTGGGGGGGCGGGGGTTGGTGGCGATTGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAgg                                                                                                                                                                                                                                                     <  2:296347/295‑1 (MQ=255)
                                                               attgaaaGATATTTTT‑GTGCTTGAATAACAAACGGCACAGGTGTGGGGGGCGCAAAAA‑TCGCCGAGTCGTTCCCTGAAGGCTGGGACACCCCACATGGATGTGGGCGGGGGGCGCTTTTGACGGGGGTGTCACCTCGCCCCCGACCTGGTAGCTGGAAGGGAGAGTTCGAGGG‑CACCGTGCAGTGGTGGTTTTGTTCGCAAGAGCTGGGGGGGGGGGGGGGTGGTGGCGTTTGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaacaa                                                                                                                                                                                                                                           <  1:646572/296‑1 (MQ=255)
                                                                         gcttctttctttttaagtacaggcttcACCGATTTTGGTGTC‑TCAAAAATCGCTGGGGCGTTCCCGGAATTCCGGGGCAGCCCACGGGGATGTGGTTTGGGGTCGTGTTTGACAGGACTGTTACCTCGCTCCCAACCTGGTACCCGGAAGGGATGATTAGAGGG‑CACCGCGCAGTGTCGATTTTGTTCACAAGAGCCTGGGGGTGCGGGGGGCGGCGTCGCTTTTCTGCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATTAAAGGAGCGCAACAAGATGATGaa                                                                                                                                                                                                                                  <  2:191740/274‑1 (MQ=255)
                                                                                     tgaatgacaaacgtcACCGGTGTTGGTGTCGTCAAAA‑TCACCGAGTCGTTCCTTAAATGCCGGGGCAACCCACAGGGGTGGGGGTGGAGGGCGCGTTTTAAAGGGGTTTCACCTCTCCCCCGACCCGGTAGCCGGAAGGGAGAAGTCGGGGG‑CACCGCGCAGTGGCGATTTTGTCCACAAGAGCTCGGGGGTGCGGGGGGTGGCGGCGGTTGGCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaac                                                                                                                                                                                                                      <  1:214282/286‑1 (MQ=255)
                                                                                                        aTGTTGGTGGCGTCAAAA‑TCGCCGAGGCGTTCACTGGAGGCCGGGGCAGCCCACATGGATGTGGGCCGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCAAGACCCGGTAGCCGTCAGGGATAAGTTGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCACGGGGGTGCAGAGGGCGGCAGCGATTGGCCGCCCCCTACGCGCCCATCGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACCAGATGCTTAACTAACATATACAAAAATCTTAAAAAAAGGcc                                                                                                                                                                                                   >  2:236300/1‑300 (MQ=17)
                                                                                                         tGTTGGTGGCGTCAAAA‑TCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAAAAACATCTTAAAAAAAggact                                                                                                                                                                                                  >  1:24718/1‑297 (MQ=255)
                                                                                                                     taaaattcaccgagtctttCCCGGAAGGCCGGGGCACCCCACGTGGGTGTGGGTTGATGGTGCGTTTTAAGGGGGTGTTACTTCTCCCCCGACCCGGGATCTGGAGGGGGTAATTCGAGGA‑CACCCCCAATTGGCGTTTTTGTTCGCCAGATCTTGGGGGTGCGGGGGGTGGTGGCGGTTTTCCCCCCCTTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGcc                                                                                                                                                                                       <  1:268058/283‑1 (MQ=255)
                                                                                                                                                         gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcc                                                                                                                                                   >  1:321128/1‑299 (MQ=255)
                                                                                                                                                         gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagca                                                                                                                                                   >  1:401376/1‑300 (MQ=255)
                                                                                                                                                         gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTG‑CCGCCCCCTGCGCGCTCATTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag                                                                                                                                                  >  1:164514/1‑300 (MQ=255)
                                                                                                                                                          cccacggggatggggttggggggCGTGTTTGACGGGGATGTCCCCTCGCCCCCGCCCCGGGAGCTGGAAGGGAGAAGTCGAGGG‑CACCGCACAGTGTTGATTTTGTTCACAAGAGCTTGGGGGTGGGGGGGGTGGTGGCGGTTGGCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCcagcag                                                                                                                                                  <  1:602085/280‑1 (MQ=255)
                                                                                                                                                           ccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCNCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTCAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCAAGCTTTTCCAGCAgg                                                                                                                                                 >  1:493015/1‑300 (MQ=255)
                                                                                                                                                              caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATCTAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTCATCACGCTTTTCCCGCAGGGTc                                                                                                                                              >  1:117128/1‑300 (MQ=255)
                                                                                                                                                                  gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGTCTCAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTCAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTCATTAATCACGCTTTTCCAGCAGGGTCTGGt                                                                                                                                          >  1:295892/1‑300 (MQ=255)
                                                                                                                                                                  gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGt                                                                                                                                          >  1:237222/1‑300 (MQ=255)
                                                                                                                                                                  gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTACGCCAGTGGCAATATATTGCTTAGCTCATGAAAGGAGCGCAAAAAGATGATGAATCAACATCTAGCAACATCTTCAAAAAAGGCCTGACCTTACGCCAGGCCTTCTGCGTTAATTAATCAAGATTTTCCAGCAGGGTCTGGt                                                                                                                                          >  2:241983/1‑300 (MQ=21)
                                                                                                                                                                  gATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCAGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGTGCAACAAGATGATGAATCAAAATATAACAACATCTTAAAAAAAGGCCTGAAATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGt                                                                                                                                          >  1:294850/1‑300 (MQ=255)
                                                                                                                                                                  aaTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGTGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCATTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGt                                                                                                                                          >  2:257697/2‑300 (MQ=255)
                                                                                                                                                                       gggttggggtcgcgTTTGACGGGGGTGTTACCCCTCTCCCGACCCGTTAACCGGAGGGAACAAGCCGAGGA‑CACCTCCCATCGTCGCTTTTGTTCACCAGACCCTGGGGGTGCAGGGGTTGGCGGCGATTGTCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATc                                                                                                                                     <  1:99452/290‑1 (MQ=255)
                                                                                                                                                                       gggttGAGGGCGCGTTTTACAGGGATGTCACCCCGCGCCAGACCCGGTCGCCGTAAGGGAGAGGTCGAGGG‑CACCGCGCAGCGGCGGTTTTGTTCGCCAGAGCCTGGGGGGGCCGGGGGTGGTGGCGGTGGGCCTCCCCTTGCGCGCTNCTTGCGCCATTTGCAATATGTTGCTTNGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATc                                                                                                                                     <  2:490151/296‑1 (MQ=255)
                                                                                                                                                                                       ttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGANTCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATACTgc                                                                                                                     >  1:506410/1‑300 (MQ=255)
                                                                                                                                                                                               aTGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCAAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCAGCCGACa                                                                                                             >  1:733970/1‑300 (MQ=255)
                                                                                                                                                                                                                 cctgggagcgggaaGGGAGAAGTCGGGGA‑CACCACGCAGTGGCGGTTTTGTTCGCCAGATCCTGGGGGGGCAGGGGGTGGCGGCGGTTGACCCCCCCCTTCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCagaag                                                                                           <  1:708550/288‑1 (MQ=255)
                                                                                                                                                                                                                 cccgggaacctgaaGGGATAGGTCGAGGG‑CACCACCCAGTGTCGATTTTGTTCGACAAAGCCCGGGGGTGCAGGGGACGACGTCGATTGCCCCCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTTACATTACGCCATGCCCTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGACAATGCCGAATACAATTGGCACCACCTagaag                                                                                           <  2:276159/288‑1 (MQ=255)
                                                                                                                                                                                                                 cccagTAGCCGTAAGGGATAAGTCGAGGG‑CACCCCGCATCGGCGATTTTGTTCGCCAGATCCCGGGGGTGCAGGGGGCGTCTTCGATTGTCCGCCCCCTACGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCCGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCTagaag                                                                                           <  2:298326/296‑1 (MQ=255)
                                                                                                                                                                                                                          cggaaGGGATAAGTAGAGGA‑CACCGCGCAGTGGTGGTTTTGTTCGCAGGATCCTGGGGGTGGGGGGGGTGGTGGTGGTTGGCCTCCCCTTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgt                                                                                  <  1:617786/297‑1 (MQ=255)
                                                                                                                                                                                                                                ggATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCt                                                                            >  1:567778/1‑300 (MQ=255)
                                                                                                                                                                                                                                      gTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAGACACAGTTGTTCGAATGca                                                                      >  2:573612/1‑299 (MQ=255)
                                                                                                                                                                                                                                             g‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAgccgc                                                               >  1:224893/1‑300 (MQ=255)
                                                                                                                                                                                                                                                      aagggtaggttttGTTCGCAAGATCCTGGGGGGGAGGGGGGCGGCGTCGATTGCCCCCCCCCTCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAg                                                       <  1:592983/291‑1 (MQ=255)
                                                                                                                                                                                                                                                         cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAAAAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATTACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGACGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATa                                                    >  1:132684/1‑300 (MQ=255)
                                                                                                                                                                                                                                                         cggcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTATGCGTTCCTTGCGCCAGTGGCCATATGTTGCTTAGTTCATGAAAGGAGCGCCAAAAGATGATGAATCCACATATAAAAAAATCTTAAAAAAAGGCCTGACATTCCGACAGGCCTTCTGCGTTAATTAACTACACTTTTCCCACAGGGTCAAGCCAATCAGACCACCAATCATGCCGCCAACAATTGGCACCAACCAGAAGAACCACCATTGTTCTCATGCCCAgccgtactggcatata                                                    >  1:132514/1‑288 (MQ=19)
                                                                                                                                                                                                                                                                      cgcaaggtctggggggTGAGGGGGGCGGCGGTGGTTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGtgc                                       <  1:329284/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gagCCCGTGGGTACAGGGGGCGTCGGCGATTGGCCACCCCCTACGCGCTCCTTGCGCCAGTTGCAATATTTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACACTACGCCAGGCCCTCTGCGTTAATTAATCACTCTCCTCCCGCAGGGTCCGGTAACTCAGACCACCGATAATGCCGCCGACAATTGGCACCACCTAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTgcgc                                  <  2:389706/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gcctgggggTGGGGGGGGTGGCGGCGATGTTCCCCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTgcgcgc                                <  1:60783/296‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  ggggTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGt                           >  1:149983/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                       gCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCACTGGGAGATAGCAAACGCGGTGCTGCGCGCCGGGTTAACa                      >  1:16729/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                            gggCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCAACCAGAAGAACCACAGGTGGCCTACTGCACAGACGACATGGAAGATCGCAACCGCAGTGATGTGCGGCGGGTTAAGAGAAgt                 >  2:112304/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                           gccgccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAAt  >  1:455811/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                        
GACGAGTAAAGTTGCAATTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTT‑GCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAA‑TCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGG‑CACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAAT  >  NC_000913/914921‑915516

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑