Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	30.4%	79.8 / 13.4	58	G566G (GGT→GGG)	gyrA	DNA gyrase subunit A
Reads supporting (aligned to +/- strand):  ref base A (32/7);  new base C (0/17);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 3.53e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTCCG  >  NC_000913/2337480‑2337964                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  gAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCCGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCCCGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGATGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑GTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTcc                                                                                                                                                                                                                                            >  1:49640/1‑251 (MQ=255)   aGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGTCAG‑ATTT‑ACCTTTCCCg                                                                                                                                                                                                                                          <  2:212423/251‑1 (MQ=255)    gACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGTGGCAG‑ATTT‑ACCTTTCCCGc                                                                                                                                                                                                                                         <  1:326778/251‑1 (MQ=255)     aCTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGGTACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCGGGAGAAGCACAGAATATGGTCGTGAGGGTTCGCCACCAGCAGT‑CGGGCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGcc                                                                                                                                                                                                                                        >  2:130708/1‑251 (MQ=255)     aCTTTCACGCCTGCTTCAAACTCGGTCACTGGCAGGATCGCAGTTATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCAGCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGcc                                                                                                                                                                                                                                        <  2:274106/251‑1 (MQ=255)          cACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGCGATACGTTCGTCCTGCTCCAGCGGCAACAGGATGACG‑ATCGGACGACCGCGCGCGCCACGAGAGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAATGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCacgac                                                                                                                                                                                                                                   >  1:409068/1‑251 (MQ=255)             gCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTCCCCGCCACGACGCt                                                                                                                                                                                                                                >  2:112004/1‑251 (MQ=255)                  cttcAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGCGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGGGCGTGTTCGCCACCAGCAGT‑CGGTCGCTAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑TCCTTTCCCGCCACGACGCTGCGCt                                                                                                                                                                                                                           >  2:166892/1‑251 (MQ=255)                                     cAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGCACGACCGCGCGCGCCCCGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACAAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTAGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCt                                                                                                                                                                                                        >  1:368806/1‑251 (MQ=255)                                      aggaccgcattgatcctttcccccTGCCCCAGGGGCAGCGGGTTGACG‑ATCGGACGCCCGCGCGCGCCACGAGGGGCTTCCGGAAACTGAAAAC‑CTTTCATCGAATAGACGCGCCCGC‑GGCTGGAAAAGCAAAAAATATGGTGGGGAGTGTTCGCCCCCAGC‑GTGCGGTCAAAAAATTC‑TTCTTTTTTAATCCG‑GGCGGCAG‑ATTT‑CCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCgg                                                                                                                                                                                                          <  1:508225/226‑1 (MQ=255)                                      aGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGCCG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCgg                                                                                                                                                                                                          >  2:508225/1‑248 (MQ=255)                                            gCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACtt                                                                                                                                                                                                 >  1:34023/1‑251 (MQ=255)                                              aGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTa                                                                                                                                                                                                >  1:452126/1‑250 (MQ=255)                                               gggatgcGTTCGTCCTCCGCCGGCGGAGGCGGGTGGACG‑ACCGGACGCCCGCGCCCGCCACGGGTGGCTTCCGGCA‑CTTATAAAACTTTCATGAAATAGCCGCGCCCAC‑GGCTGGAGAAGCACAGAATATGGTGGGGAGTGTCCCCCCCCAGCAGT‑CGGTCGAAAAAGTC‑TTCTTCTTTAATCCG‑TGCG‑CAA‑ATTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc                                                                                                                                                                                              <  2:121912/245‑1 (MQ=255)                                               gTGATGCGTTCGCCCTGCCCCAGCGCCAGCAGGTT‑ACGGATCGGCCGCCCGCGCGCGCCACGAGGGGCTTCCGGAAACTTATAAC‑TTTTCATCGAATAGACGCGACCAG‑GTCTGGAAAAGCACAAAATATGGTCGTGAGTGTCCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTTTTTAATACG‑TGCGGCAG‑ATTT‑CCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAAc                                                                                                                                                                                              <  1:437321/251‑1 (MQ=255)                                                   tACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTCACGt                                                                                                                                                                                            >  2:465968/1‑249 (MQ=255)                                                        tCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATCAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGCGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTg                                                                                                                                                                                     >  2:200073/1‑251 (MQ=255)                                                          gTCCTGCTCCAGCGGCAGCAGGCTGACG‑ATCGGCCGACCGCGCGAGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGCATAGACGCGACCAC‑GGATGGCGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCCGT‑CGGTCGCTAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTTATCAGCAAGCGGCTGATACTTACCGTAGCCCTggg                                                                                                                                                                                   >  2:495409/1‑250 (MQ=255)                                                                       ggCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCAc                                                                                                                                                                      >  1:344071/1‑251 (MQ=255)                                                                                  tgacactcggccgcccGCGCGCGCCACGAGTGGTTTCCGGAACCTAAAAAA‑CTTTCATCGAATAGCCGCGCCCAC‑GGCGGGAGAAGCAAAAAATGTGGTGGGGAGTTTCCGCCCCCAGCGGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑GGCGGCAA‑TTTT‑CCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATc                                                                                                                                                             <  1:524345/236‑1 (MQ=255)                                                                                   aCG‑ATCGGACGACCGCGCGCGACCCGAGTGGCTTCCGGCCACTGATAAC‑GTTTAATCGAATAGAAGCGACCAC‑GGCTGGAGAAGCACAGTATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTATTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATCCTTAACGTCGCCCTGGTGGGAGAGCGTc                                                                                                                                                                        >  2:297965/1‑237 (MQ=255)                                                                                         gggACGACCGCGCGCCCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAAGAGCCGCGACCAC‑GTCTGGAGAAGCACAGAATATGGTCGGGAGTGTTCCCCACCAGCAGT‑CGGTCAAAAAATTC‑TTCTTTTTTAATTCG‑TGCGGCAT‑ATTC‑CCCTTGCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTgg                                                                                                                                                      <  1:77986/249‑1 (MQ=255)                                                                                             cgacCTCGCGCGCCCCGATTGGCTTCCGGCAACTGATAAC‑CTTTCATCGAAGAGACG‑GACCCCCGGCTGGAAAAGCAAAAATTATGGTCGTGAGTTTTCGCCCCCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATGCG‑GGCGGCAG‑ATTT‑CCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGAt                                                                                                                                                 <  1:130708/251‑1 (MQ=255)                                                                                              gcccGCGCGCGCCAGGGGTGCCTTCCGGAAACTAATAAA‑CTTCAACCAAATAGACGCGACCCC‑GGCTGGAAAAGCAAAAAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑GGGTCGAAAAAGTC‑TTCTTCTTAAATACG‑TGCGGAAG‑TTTT‑CCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                                                                                                                <  1:479967/249‑1 (MQ=255)                                                                                              gacCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                                                                                                                >  1:326458/1‑251 (MQ=255)                                                                                                  gcgcgcgcCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGGATAGACGGGACCAC‑GGCTGGAGAAGCACAGCATATGGCAGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAACGCGGCTGATACTTAAc                                                                                                                                                                                              >  1:73491/1‑201 (MQ=255)                                                                                                  gcgcgcgcCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGGCACGACCACATCTTCCTGGGTGATc                                                                                                                                                >  2:422539/1‑247 (MQ=255)                                                                                                          cACGGGTGGCTTCCGGAAACTGAAAAC‑CTTTCATCGAATAGACGCGACCCC‑GGCGGGAAAAGCAAAGAATAGGGTCGGGAGTTTCCGCCACCAGAAGT‑GGGTCAATAAATTC‑TTCTTTTTTAATCCG‑TGCGGCAG‑TTTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAgg                                                                                                                                    <  1:390305/251‑1 (MQ=255)                                                                                                                 ggtttCCGGCAATTGATATA‑TTTTCATAGAACAGACGCTCCCAC‑GGCTGGAGGAGCACAAAAAAAGTTCGTAAGTTTCCCCCCCCCGCCGT‑CGTTCTAAAAATTT‑TTCTTTTTTACTCGG‑TGCGAAAG‑TTTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGt                                                                                                                             <  1:39955/248‑1 (MQ=255)                                                                                                                   cTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCCCGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCt                                                                                                                           >  2:442813/1‑251 (MQ=255)                                                                                                                            acttaTAAA‑CTTTCATCAAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTTTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑GGCGGCAG‑ATTT‑CCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                                                                                                                  <  1:97358/247‑1 (MQ=255)                                                                                                                             acggaaaaccTTTCATCGAATGGACGCGACCAG‑GGCGGGAGAAGCAAGGAAGATGGTCGTGGGTGTTCGCCCCCAGCAGG‑CGGTCGATAAAGTC‑TTCTTCTTAAAGACG‑TGGGGCAG‑ATTT‑GCTTGTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACTACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                                                                                                                  <  1:303252/242‑1 (MQ=255)                                                                                                                              attaaaatcTTTCATCGAATAGACGCCACCAC‑GGCTGAAGAAGCACAAAATATGGGCTGGGGTGTTCGCCCCCAGCAGC‑CGGTGAAAAAAGTC‑TTCTTCTTTAATACGGTG‑GGAAG‑ATTT‑CCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTgg                                                                                                                 <  2:34023/243‑1 (MQ=255)                                                                                                                                    a‑CTTCCATCAAATAGCCGCGCCCCC‑GGCGGGAGAAGCAGAGAATATGGTCGTGAGGTTTCGCCACCAGCAGG‑CGTACGATAAATTC‑TTCTTCTTAAATCCG‑TG‑GGGAGTATTT‑CCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAt                                                                                                          <  1:224398/251‑1 (MQ=255)                                                                                                                                        ttCATCAAAGAGACGCGCCCAC‑GGCTGGAGAAGCACAGAATATGGTGGTGAGTGTTCGCCCCCAGCAGT‑GGGTCGATAAAGTC‑TTCTTCTTTAATGCG‑TGGGGCAG‑ATTT‑CCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAtt                                                                                                         <  1:377633/249‑1 (MQ=255)                                                                                                                                          cATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTc                                                                                                       >  1:90256/1‑249 (MQ=255)                                                                                                                                               aaTAGACGCGACCAC‑GGCTGGAGAAGCAAAGAATATGGTCGTGAGTGATCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGATGATGTCTGCGCTGTTGGCGGTGATTTCAGTacg                                                                                                 >  1:214966/1‑250 (MQ=255)                                                                                                                                                aatagccgcgcccggggCTGGAGAAGCACAAAATGTGGTCGTGAGTGTCCGCCACCGGCGGG‑CGGT‑GAAAAAATCCTTCTTCTTAAATGCG‑TGCGGCAG‑ATTT‑CCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTacga                                                                                                <  1:114280/236‑1 (MQ=255)                                                                                                                                                  aGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCt                                                                                                                      >  1:257397/1‑226 (MQ=255)                                                                                                                                                                   gagaaaccacaGAATAGGGTCGGGAGTGTCCCCCACCAGCGGT‑CGGTCGATAAAGTC‑TTCTTCTTAAATACG‑GGCGTCAG‑ATTT‑CCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGtt                                                                             <  1:260324/244‑1 (MQ=255)                                                                                                                                                                    gagaAGCACAGAATATGGTCGTGAGTGTTCGCCAACAGCAGT‑CGGTCGCTAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAACGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACAACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGGTTTCAGTACGCCGTTTGTCACCGAACGGTTc                                                                            >  2:524567/1‑251 (MQ=255)                                                                                                                                                                      aaaGCACAAAAAAGGGTCGTGCGTGTTCCCCCCCAGCGGC‑CGGTCGATAAAGTC‑TTTTTTTTTAATGCG‑TGGGGAAG‑TTTT‑CCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCAc                                                                          <  1:200073/250‑1 (MQ=255)                                                                                                                                                                             aGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAAcc                                                                     >  2:205538/1‑249 (MQ=255)                                                                                                                                                                                           aGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCCCGACCCAGCTCCAGCTCTTCAc                                                     >  2:29393/1‑251 (MQ=255)                                                                                                                                                                                              gtTCGCCACCAGCAGT‑CGGTAAATAAAGTC‑TTTTTCTTTATTTCG‑TGCGCAAG‑TTTT‑CCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGa                                                  <  1:509885/251‑1 (MQ=255)                                                                                                                                                                                                           agT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATc                                       >  1:192201/1‑249 (MQ=255)                                                                                                                                                                                                                   cgATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGACACGCCGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGCAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGCGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCCGTCCGACGTTTGTCACCGAACTGTTCACGGACAAGCTCCCGCTCTTCACGGATCACTTCCAACAGACGATCg                              >  2:357507/1‑251 (MQ=255)                                                                                                                                                                                                                          agttttcttcttTAATCCG‑TGCGGCAG‑ATTT‑CCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTa                        <  1:422539/247‑1 (MQ=255)                                                                                                                                                                                                                            c‑ttcttcttTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACTCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGATGATGTCTGCGCTGTTGGCGGTGATTTTAGTAAGACGTTTGTCACCGACATGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGAACGCCGCTAcc                      >  2:264234/1‑250 (MQ=255)                                                                                                                                                                                                                                  tcttTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGaa                 >  2:132570/1‑250 (MQ=255)                                                                                                                                                                                                                                         tACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAAc         >  1:511803/1‑251 (MQ=255)                                                                                                                                                                                                                                          agggtgaggCAG‑TTTT‑ACCCGTCCCGCCACGACGCTGCGCTTCGGATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCCCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAAc         <  1:495409/244‑1 (MQ=255)                                                                                                                                                                                                                                           cG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAAc         >  1:563021/1‑249 (MQ=255)                                                                                                                                                                                                                                           cG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAAc         >  2:279565/1‑249 (MQ=255)                                                                                                                                                                                                                                             gggCGGCAG‑ATTT‑CCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGt      <  1:510222/249‑1 (MQ=255)                                                                                                                                                                                                                                              gggggCAG‑ATTT‑ACTTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGtt     <  1:290910/248‑1 (MQ=255)                                                                                                                                                                                                                                                 ggCAG‑CTTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCCCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCCAGAATACGCAACAGTTCcg  >  2:482148/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ggCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGCCGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCCTCAGACGAACGGCCCCACCAAGAATACGCAACAGGTccc  >  2:87289/1‑250 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                  GAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACG‑ATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGCAACTGATAAA‑CTTTCATCGAATAGACGCGACCAC‑GGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGT‑CGGTCGATAAAGTC‑TTCTTCTTTAATACG‑TGCGGCAG‑ATTT‑ACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTCCG  >  NC_000913/2337480‑2337964

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