Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,578,143	0	A	C	30.2%	129.5 / 20.9	64	L79V (TTA→GTA)	maeB	malate dehydrogenase (oxaloacetate‑decarboxylating) (NADP(+))
Reads supporting (aligned to +/- strand):  ref base A (38/6);  new base C (0/19);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 2.89e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAATTCATGGAAATCAAGTGCACTTTGTTTTAACTGGTCATCCATTT  >  NC_000913/2577894‑2578379                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              cgtcgtGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAAc                                                                                                                                                                                                                                              >  1:596130/1‑251 (MQ=255)  gtcgtGGAATCCCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCACA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAAcc                                                                                                                                                                                                                                             >  2:47395/1‑251 (MQ=255)  gtcgtGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAAcc                                                                                                                                                                                                                                             >  2:240683/1‑251 (MQ=255)    cgtGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAAccc                                                                                                                                                                                                                                            <  2:145291/250‑1 (MQ=255)       ggAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGc                                                                                                                                                                                                                                       >  1:313454/1‑251 (MQ=255)          atagcGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAGCATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CAAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTTGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCGTCC‑ATCACCGGTTT‑GAC‑AGCCAGCGCGTCAAT‑G‑T‑TGCCTAACCCC‑AGctcc                                                                                                                                                                                                                                    <  2:185458/247‑4 (MQ=255)          aTACCGGAATGTTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAACCATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCTAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCATAACCCT‑AGCAcc                                                                                                                                                                                                                                    <  2:362747/251‑1 (MQ=255)                            gCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTag                                                                                                                                                                                                                       >  2:434982/1‑246 (MQ=255)                               cgcgcgttttCTGTTAAATAGAAAAACATTCGGGCGCTTTAGTGTCGTAGAGGTTGATCCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGAAAACTTAAATAAATTTGCCGGGGTCGAGTCCGTCAACTTCAATGTAAAA‑TCCATCACCCCCGGCAAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ACCCCCGTTTT‑GCC‑GGCCAGCGCGCCAAT‑G‑T‑TGCCAACCCCC‑AGCACCGCCGTACCGTTAGAGATCAcc                                                                                                                                                                                                               <  2:19689/244‑1 (MQ=255)                                              cAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCt                                                                                                                                                                                                >  2:131359/1‑251 (MQ=255)                                               aaTATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTc                                                                                                                                                                                               >  2:353212/1‑251 (MQ=255)                                                      aaaCATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                                                                                                                        >  2:670692/1‑251 (MQ=255)                                                        agccgtttgctttcccgatctttCGGGGTTGTTCCC‑GC‑TGAGGGTTGTTTCGAGCGCGGCGACAATTTAAATAATTTTCCCGGGGTCGATTTGTTCAACTCAAATGTAAAA‑TACATCAACCCCGGCGATTTTTTTAAA‑AA‑GAA‑C‑CCCCTGGCCTCCC‑ACCCCCGTTTC‑CCC‑TCCCAGCGGGCCATG‑G‑T‑TCCCTACCCCC‑ACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                                                                                                                      <  1:608884/230‑1 (MQ=255)                                                        aCATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                                                                                                                      >  2:477690/1‑251 (MQ=255)                                                         cATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                                                                                                                      >  1:168897/1‑250 (MQ=255)                                                          aTTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                                                                                                                    >  1:372949/1‑251 (MQ=255)                                                           ttcgggCGTTTTATTGTCTT‑GAGGTTGGAGCCCGC‑GGAAGGTTGGTCCGAGCGGGGCACCAATTTAAATAATTTTGCCGGGGTCGAGTTCGCAAACTTCAATTTAAAA‑TACTTCAACCCCGGCGATTTTCTAAAA‑AA‑GAC‑C‑GCCCTTCCTTTCC‑TTCCCCGGTTT‑GCC‑TGCCAGCGCGCCATT‑T‑T‑TCCCTACCCCC‑GGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                                                                                                                    <  1:594769/246‑1 (MQ=255)                                                                gcgttttATTGTCTGCGGGGTTGATGCC‑GC‑GGAGGTTTGGTTCGAGCGCGGCAACAATTAAAATAATTTGGTCCGGTTCGAGTCCGTAAACTTAAATGAAAAA‑TCCATCAACCCGGGCGAATTTCT‑AAA‑AA‑GAA‑CCGCCCT‑GCCCTCC‑CACCCCGGGTT‑GCCCGGCCGGCGC‑CCAATTG‑T‑TGCCGACCCCA‑AGGACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAg                                                                                                                                                                              <  2:465631/247‑1 (MQ=255)                                                                   cTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCt                                                                                                                                                                           >  1:381205/1‑251 (MQ=255)                                                                     ttAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                                                                                                                         >  2:77844/1‑251 (MQ=255)                                                                      aaaTGCCTTCGAGGTTGA‑GCC‑CCGCGGAGGTGGGTTGGAGCGCGGCACCATCTTAAAAAATTTTGTCGGGGTCGGTTTCTACACCTACAATGTCAAA‑ACCATCATTCCCGGCAATTTTTTTAAA‑C‑‑GAAGC‑CCCCTTCCCTCCC‑ATCACGGGTTT‑GCC‑TCCCAGCCCCCCATT‑G‑T‑TCCCAAACCCC‑ACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                                                                                                                         <  1:692129/249‑1 (MQ=255)                                                                             ttCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                                                                                                                      >  2:608884/1‑230 (MQ=255)                                                                               cgagtttaaTGCC‑GC‑CGAGGGTTGGTTCGGGCGCGGCAACAACTTAAATAAATTTGTCCGGGTCGGGTTCTTCACCTTCAATGCCAAA‑TACACCAACCCCGGCGAATTTCTAAAA‑CA‑GAC‑C‑CCCCTTGCCTTCC‑ACCCCCGTTTG‑GCC‑TGCCAGCGCGCAATT‑G‑T‑TCCCTACCCCC‑AGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTc                                                                                                                                                                <  1:249720/242‑1 (MQ=255)                                                                                    ttGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                                                                                                          >  2:613266/1‑251 (MQ=255)                                                                                      gATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                                                                                                                         >  2:33597/1‑250 (MQ=255)                                                                                       gtccc‑gc‑cgAGGGTTGGTTCGGGCGCGGCGACAACTTAAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTAAAA‑TACATCAACCCCGGCGAATTTTTAAAA‑CA‑GAC‑C‑GCCCTTCCCTCCC‑ACCCCCGGTTT‑GCC‑TGCCGGCGCGCCATT‑G‑T‑TGCCAACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                                                                                                                       <  1:670692/248‑1 (MQ=255)                                                                                           cccgcggagggttggttCGAGCGGGGCGAAAACTTAAAAAAATTTGCCCGGGTAAAGTCCGACAACTTCAAT‑TCAAAATCCACCACCCCGGGCGATTTTCTAAAA‑CA‑AAC‑C‑CCCCTTGCCTCCC‑ACCACCGTTTT‑GCC‑TGCCGGCGC‑CCAAT‑G‑G‑TGCCAAACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt                                                                                                                                                     <  2:168897/242‑1 (MQ=255)                                                                                             cgcagAAGGTTGGTTCGGGCGCGGCAACATCTTCAATAAATTTGTCCGGGTGGAGTTCGTCAACTTCAATGTAAAA‑TACATCAACCCCGGCGAATTTCTAAAA‑CA‑GAC‑C‑GCCTTTGCCTTCC‑ATCCCCGTTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                                                                                                                   <  1:6175/247‑1 (MQ=255)                                                                                                       tggttggggcgcGGCGACAACTTAAATAAATTTGGCGGGGTCGAGTTCATCAACTTCAATTTCAAA‑TACATCAACCCCGGCAAATTTCTAAAA‑CA‑GAC‑C‑GCCCTTCCCTTCC‑ATCCCCGTTTT‑GCC‑GGCCAGCGCGCAAAT‑G‑T‑GGCCAACCCCC‑GGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                                                                                                                          <  1:34936/242‑1 (MQ=255)                                                                                                         gttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCg                                                                                                                                       >  1:407885/1‑251 (MQ=255)                                                                                                            cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAgg                                                                                                                                                                             >  1:446387/1‑210 (MQ=255)                                                                                                             caGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCa                                                                                                                                    >  2:240559/2‑250 (MQ=255)                                                                                                                  cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCt                                                                                                                                                               >  2:173420/1‑218 (MQ=255)                                                                                                                         aacttaaaaaattttgccggggTCGATTTCGTCAACTTAAATGTCAAA‑TACATCACTCCCGGCAAATTTCTAAAA‑CA‑GAC‑C‑GCCCTT‑CCTTCCCATCCCCGTTTT‑CCC‑TCCCAGCGCGC‑AAT‑GGT‑TGCCTACCCCC‑AGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                                                                                                                       <  1:477690/232‑1 (MQ=255)                                                                                                                            ttCAATAAATTGTTCCGGGTCGAGTCCTTCACTTTAAATTTAAAA‑TACATCAATCCCGGC‑AATTTTTTAAAACG‑AAA‑C‑CCCCTT‑CCTTCC‑CACCCCGGTTTTCCC‑TGCAGGCCCGCAAAT‑G‑T‑TGCCAACCCCC‑AGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAg                                                                                                                     <  1:79372/250‑1 (MQ=255)                                                                                                                               aaTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc                                                                                                                  >  1:508075/1‑250 (MQ=255)                                                                                                                                   aaTTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcgc                                                                                                             >  1:389399/1‑251 (MQ=255)                                                                                                                                      tggccggggTGGAGTTCGTAAACTTAAATGAAAAA‑TCCATCAACCCGGGCAATTTTCTAAAA‑AA‑GAC‑C‑GCCTTTCCCTCCA‑ACCACGGGTTT‑GCC‑GGCCAGCGCCCAAAT‑T‑T‑GGCCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCa                                                                                                           <  1:92273/244‑1 (MQ=255)                                                                                                                                        gTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCAATTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑G‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGACACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTGGATTTCAAGACAAGGTGCGTAAACGCCTGATGAGTAGGCCAGCGcc                                                                                                            >  1:624913/1‑247 (MQ=255)                                                                                                                                           cggttCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgc                                                                                                     >  1:177360/5‑251 (MQ=255)                                                                                                                                           cGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgc                                                                                                     >  1:177357/1‑251 (MQ=255)                                                                                                                                            gggTCGATTTCGTCAACTTCAATGTAAAA‑TACATCAACCCCGGCGAATTTCTAAAA‑CA‑GAC‑C‑CCCCTTGCCTCCC‑A‑CCCCGGTTTTGCC‑TCCCGGCGCGCCATT‑G‑T‑TGCCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg                                                                                                    <  1:162381/251‑1 (MQ=255)                                                                                                                                             ggTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCgcgc                                                                                                   >  1:75354/1‑251 (MQ=255)                                                                                                                                                   gTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGtt                                                                                             >  1:346580/1‑251 (MQ=255)                                                                                                                                                       gTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGcc                                                                                          >  2:429091/1‑250 (MQ=255)                                                                                                                                                           aCTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAgg                                                                                     >  2:333407/1‑251 (MQ=255)                                                                                                                                                              aaaaTGACAAA‑TACACCAACCCCGGCAAATTTCT‑AAA‑CAAGAC‑C‑CCCCTTGCCTTCC‑ACCCCCGGTTT‑CCC‑TGCCCGCGCGCCATT‑T‑T‑TGCCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGCTGCCAGAGGCtt                                                                                  <  1:497072/249‑1 (MQ=255)                                                                                                                                                                 aTGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCttggt                                                                               >  1:210385/1‑251 (MQ=255)                                                                                                                                                                  tGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCttggtt                                                                              >  2:585520/1‑251 (MQ=255)                                                                                                                                                                   gtaaaa‑TACATCAACCCCGGCAATTTTCTTAAA‑CA‑GAC‑C‑CCCTTTCCCTCCC‑ACCCCCGTTTT‑CCC‑TGCCAGCGCCCCATT‑G‑T‑TGCCAACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCttggttg                                                                             <  1:481046/248‑1 (MQ=255)                                                                                                                                                                             aTCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                                                                                                                        >  1:546986/1‑135 (MQ=255)                                                                                                                                                                             aTCAACCCCGGCAATTTTCTTAAA‑CA‑GAA‑C‑GCCTTTCCCTTCT‑ACCCCCGGTTT‑TCT‑TCCCAGCGCGC‑AAA‑G‑TGTGCAAAGCCCC‑AGCTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCt                                                                    <  1:153353/251‑1 (MQ=255)                                                                                                                                                                                    ccGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGAtttt                                                             >  1:91412/1‑251 (MQ=255)                                                                                                                                                                                                   aa‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCctggaactggaa                                               >  1:327351/1‑250 (MQ=255)                                                                                                                                                                                                   aa‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCctggaactggaa                                               >  2:257209/1‑250 (MQ=255)                                                                                                                                                                                                         cagaccccccTTGCCTCCC‑ACCGCCGGTTG‑CCC‑TGCCGGCGCCCCAAT‑G‑T‑TCCCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAAtt                                            <  1:131359/244‑1 (MQ=255)                                                                                                                                                                                                          aa‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGGTGGAGAAACATGGATTTTCCCTGGAACTGGAAATTCAt                                         >  1:216266/1‑251 (MQ=255)                                                                                                                                                                                                                              tCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAATTCATGGAAATCAAGTGCACtt                        >  1:228723/1‑251 (MQ=255)                                                                                                                                                                                                                               cccgggTTT‑GCC‑CGCCAGCGCGCCATT‑G‑T‑TGCCAACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAATTCATGGAAATCAAGTGCAcn                         <  2:327351/245‑2 (MQ=255)                                                                                                                                                                                                                               cccgggTTG‑GCC‑CGCCAGCGCGCCATG‑G‑T‑CGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAATTCATGGAAATCAAGTGCACttt                       <  2:210385/247‑1 (MQ=255)                                                                                                                                                                                                                                aCCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAATTCATGGAAATCAAGTGCACTTTg                      >  1:83399/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggttttccctcccagcgcacaatt‑G‑T‑TGCCAACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAATTCATGGAAATCAAGTGCACTTTGtt                    <  1:333407/228‑1 (MQ=255)                                                                                                                                                                                                                                                    gcgcgcaatttt‑tGCCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAATTCATGGAAATCAAGTGCACTTTGTTTTAACTGGTCATcc      <  1:240683/240‑1 (MQ=255)                                                                                                                                                                                                                                                      ccaaaT‑G‑T‑TGCATACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAATTCATGGAAATCAAGTGCACTTTGTTTTAACTGGTCATCCAttt  <  1:33597/248‑1 (MQ=255)                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              CGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCC‑GC‑CGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAA‑TACATCAATCCCGGCGAATTTCTTAAA‑CA‑GAA‑C‑GCCCTTGCCTTCC‑ATCACCGGTTT‑GCC‑TGCCAGCGCGCCAAT‑G‑T‑TGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTGCCAGAGGCTTGGTTGGAGAAACCTGGATTTTCCCTGGAACTGGAAATTCATGGAAATCAAGTGCACTTTGTTTTAACTGGTCATCCATTT  >  NC_000913/2577894‑2578379

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