Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	26.3%	74.5 / 14.0	39	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase subunit A
Reads supporting (aligned to +/- strand):  ref base T (3/25);  new base G (10/0);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 6.05e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAAT  >  NC_000913/941471‑941960                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                             gCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACt                                                                                                                                                                                                                                                      <  2:154957/251‑1 (MQ=255)  cTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGGAACGCGAATGAGTGGCGGTGGGGTGGCttcct                                                                                                                                                                                                                                                      >  1:87771/1‑247 (MQ=255)                gACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCAc                                                                                                                                                                                                                                      <  2:43152/251‑1 (MQ=255)                gACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCAc                                                                                                                                                                                                                                      <  2:324048/251‑1 (MQ=255)                              cgccgGGAAATGCCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAAt                                                                                                                                                                                                                        <  2:378701/251‑1 (MQ=255)                                         aCCCTGGGCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGat                                                                                                                                                                                                             <  2:240320/251‑1 (MQ=255)                                                 cgcgcgGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                                                                                                                                      <  1:490657/250‑1 (MQ=255)                                                              aaCTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAc                                                                                                                                                                                        <  2:513495/251‑1 (MQ=255)                                                                        gcggcTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                                                                                                                              <  2:470067/251‑1 (MQ=255)                                                                                cctaaCCATTACGGCGACTTAC‑CCTCCGAGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                                                                                                                                       <  1:236635/247‑1 (MQ=255)                                                                                          aCGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTAATATCTTTAACAAGGAAC‑CCAGAAAACTAATGCCCGCATGATCATAATCCATCCGCCCTCTACCGACACCGGAGCCGGGCGCGGAGATg                                                                                                                                                              >  2:4197/1‑249 (MQ=255)                                                                                           cGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCGCGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGt                                                                                                                                                           <  2:371534/251‑1 (MQ=255)                                                                                            ggCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGc                                                                                                                                                                            <  1:528497/233‑1 (MQ=255)                                                                                                           cgcgCAAATTGCGGAAGGTTTGAACTATGCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCTAATTAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                                                                                                                            <  2:241665/251‑1 (MQ=255)                                                                                                           cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                                                                                                                            <  1:213080/251‑1 (MQ=255)                                                                                                            gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                                                                                                                           <  2:532042/251‑1 (MQ=255)                                                                                                                                 aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGAGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGCGGCGGGGGGGGGACCTACCATCTTGAAAA‑GGCACGCCCCAAAACTTATGGCCCCCTGATCATCCTCGAACCCCCCCCTAACGAAACCGGTGCCGGGGGCGGAGGTGAGGGGATCCCTAATTGTCCTGGGAAC‑GATGCCGAAC‑Tggg                                                                                                                    >  1:563760/1‑250 (MQ=255)                                                                                                                                               cggcTGGGCAGATGGCAACAGCCCGTCGTATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑t                                                                                                     <  2:312289/251‑1 (MQ=255)                                                                                                                                                         gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACA‑GGCACGCCAGAAAACCAATGCCCCCCTTATCATCATAGAACAGCGCTTTTAAGCCACCGGGGGCGGGCGCGGAGATGTGTGGGGCCCTTTTCGGCC‑GGGGCAAGATGCCGCTC‑TGGGTACAGGTCTgggggactggctgat                                                                                            >  2:331941/1‑238 (MQ=255)                                                                                                                                                             gCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATTTTGAACA‑GGCACGCCAGAAATCTAAAGGCCCCGTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGGTCCCTATTCTTCC‑GGGGACAGATGCCCCAC‑TGGGTAACGGTCTGGCG‑GACG‑TCATGATCAcc                                                                                      >  1:59140/1‑250 (MQ=255)                                                                                                                                                               aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTGAACA‑AGCACGCCAGAAATCTAATTCCCGCCTGATCATTATCGAAACCCGCTATACCGACACCCGGGGCGGGGGCGAAGATGAGGGGATCCCTATTCTTCC‑GGGGACAGATGCCCCCC‑CGGGTAAAGGTCTGGCG‑TCAG‑TCATGATCACTg                                                                                     >  1:130688/1‑250 (MQ=255)                                                                                                                                                                                  aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACA‑GGCACGCCCGAAAACTAATTGCCGCATGATCATCATCGATCCCCGCTATAACGACACCGGGGCCGGGGGCGAAGATGAGTGGGTCCCTATTCGTCC‑GGGGACAGATGCCGCCC‑TGGGTAACGGTCTGGCG‑TACG‑GCATGAACACTGAAAACCTGGTGGATCAGGcc                                                                 >  1:454353/1‑250 (MQ=255)                                                                                                                                                                                      cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGCTCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCCATTCGTCC‑GGGGACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCCTGATCACTGAAAACCTGGTGGATCAGGCATTcc                                                             >  1:540080/1‑251 (MQ=255)                                                                                                                                                                                       aGTAAGCTGGTAGTGCTGTTTGGGAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACA‑AGCGCGCCAGAAATCTAAAGGCCGCATGATCATCATCGAACCGCGCCCTACCCACACCGGGGCGGGGGGCGGAGAAGAGTGGGGTCCTATTCTTCC‑CGGGACAGATGCCGCACAGGGTTAA‑GGTCCGGCG‑TGCGCTCA‑GATCACTGAAAACCTGGGGGAACAGGGATTccc                                                            >  2:110277/1‑250 (MQ=255)                                                                                                                                                                                                gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCTC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAata                                                   <  1:331941/251‑1 (MQ=255)                                                                                                                                                                                                             gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTAc                                      <  1:567390/251‑1 (MQ=255)                                                                                                                                                                                                             gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACATACTAACTTGAACA‑CGGACGCCCCAAATCTAAAGGCCGCATGGCCCTCAGAGAATCGCGGCATACAGACACCGGGGCCGGGGGCGAAGATGGGTGGAGCACTCTTCGTCC‑GGGAACAGATGACGCCC‑TGGTTTACGGGCAGGCGGTAGT‑TCA‑GATCAATGGAAACATGGTGGGTCAGGGATTGCTAGATAAATATTTTGTTGGCtgc                                      >  2:493013/1‑249 (MQ=255)                                                                                                                                                                                                               aataaCCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACTTACTATCTTGAACA‑GGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGCCACCCGTGCCGGGCGCGGAGATGAGTGGGTGCCCATTTGTCC‑GGGGACAGATGCCGCAC‑TGGTTAACGGGCTGGCG‑TCCG‑GCATGATCACTGAAAACCTTGTGGTTCAGGCATTCCTCGATAAATAATGCGTTgggta                                       >  1:202781/1‑245 (MQ=255)                                                                                                                                                                                                                        ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAAc                           <  1:110277/251‑1 (MQ=255)                                                                                                                                                                                                                                    aTGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGccagccag               <  1:199560/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGcc         <  1:65665/251‑1 (MQ=255)                                                                                                                                                                                                                                          ggCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGcc         <  2:454353/251‑1 (MQ=255)                                                                                                                                                                                                                                                ggggTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTc                                                           >  2:115118/1‑195 (MQ=255)                                                                                                                                                                                                                                                ggggTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTc                                                           <  1:115118/195‑1 (MQ=255)                                                                                                                                                                                                                                                ggggTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTgcc                     <  1:104480/233‑1 (MQ=255)                                                                                                                                                                                                                                                ggggTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCATGCCAGCCAGTGCGCCGaaaa    >  1:564360/1‑250 (MQ=255)                                                                                                                                                                                                                                                ggggCGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTgcc                     >  2:104480/1‑233 (MQ=255)                                                                                                                                                                                                                                                 gggTGACTTACTATCTTGAACA‑GGCACGCCATAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTTCCAGCCAGTGCGCCGAAAAAt  <  1:344776/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gggTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAAt  <  2:540080/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                             GCTGGGCGGCACCATGACCCGCTCCTGGCCGCCGGGAAATACCCTGGTCGCGCGGCTGATGAACTGCTGCGGCGGCTATCTGAACCATTACGGCGACT‑ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACA‑GGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCC‑GGGAACAGATGCCGCAC‑TGGTTAACGGTCTGGCG‑TACG‑TCATGATCACTGAAAACCTGGTGGATCAGGCATTCCTCGATAAATATTGCGTTGGCTACGATGAGAAAACCCTGCCAGCCAGTGCGCCGAAAAAT  >  NC_000913/941471‑941960

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: ‑