Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,256,061	0	C	T	30.2%	55.9 / 12.4	44	E136K (GAA→AAA)	plsB	glycerol‑3‑phosphate 1‑O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (10/20);  new base T (13/0);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 3.34e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCC‑‑AGTTTCTGAGCGATAGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCG‑TGTAATAGGTGA‑ACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTC  >  NC_000913/4255818‑4256305                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            aCAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCC‑‑AGTTTCTGAGCGATAGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCTATTCGCCTTTTT‑TTCGcccc                                                                                                                                                                                                                                                      >  1:476761/1‑251 (MQ=255) aCAGACCCTGACGAGCAGGAAACCGGGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATCCGCGCACCGCGCGCC‑‑AGTTTCTGAGCGATAGTTTTACCCGGGCCGTGT‑TCA‑CCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGcccc                                                                                                                                                                                                                                                      <  2:756443/251‑1 (MQ=255)                     aaaCGGGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCC‑‑AGTTTCTGAGCGATAGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCAc                                                                                                                                                                                                                                <  1:101423/251‑1 (MQ=255)                                             cGTTGACGGGCAAAGGTTATTCGCCCCGCGCTCGCC‑‑AGTTCCTGAGCGATAGTTTCACCCGGGCCGTGT‑CCA‑TCCGCCATACGGCGCAGCGAACCGAACGG‑CGAGAAACCCACAAAACTGTCGCGACCGGGCCACTGTACAGCTAAAAATTTCTGTACGCCGTTATGCATACGCAGCTGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAC‑CATCACCGACAC‑TGGCACCATCTGCACATc                                                                                                                                                                                                       <  1:98933/251‑1 (MQ=255)                                                tgacggtcaaaggacctccgcccctcgcgcgCC‑‑ATTTTCTAAACAATAGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCCTCCGGCCCAGCTAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCTGCGCGCGACC‑AAA‑CATCACCTACAC‑TGGCACCATCTGCACATCCa                                                                                                                                                                                                     <  1:296093/225‑1 (MQ=255)                                                  gagtggtcaaggggcagcgtcgaCGCGCGCCCGAGTTTCAGAGCTATAGTTTTACCCGTGACG‑GTGTCA‑CCCGCCATACGGCGCCGGGAAACTGACGG‑CGAGAAACGCAAAAAACTTTCGAGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCATTGCGCCTATTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑g                                                                                                                                                                                                   <  2:690454/229‑1 (MQ=255)                                                     ggaaaaGTGCATCCGCCCCACGCGCGCC‑‑ATTTTCGGAGCTAGTGTTTTAGCCGTGCCGTGT‑CCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAACCGCACAAAACTGGCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCTTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTg                                                                                                                                                                                              <  1:407459/248‑1 (MQ=255)                                                             ccaacccccccaggcgcccC‑‑ATTTTCAGAGCGAAAGTTTTACCGGTGCCGT‑T‑TCAACCCGACATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAACCTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTa                                                                                                                                                                                      <  1:368279/238‑1 (MQ=255)                                                                 atacgagccctccgcctatttttttGAGCGTTAGTTTTACCTGTGCCGTGT‑TCT‑TCCGCCAAACGGCGCAGGGAAACTGACGG‑GGAAAACCGCACAAACCTGTGGCGACCGAGCCACAGTCCAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACg                                                                                                                                                                                    <  2:784337/227‑1 (MQ=255)                                                                  cgcgCCCCGCGCGCA‑‑TGTTGCTGAGCGTTATTTTTATCCTTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCTAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCCCTTCGCCTTTTT‑CACGCCCCGGCGCGCGACT‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTg                                                                                                                                                                                 <  1:515808/251‑1 (MQ=255)                                                                            gcgcC‑‑AGTTTCTGAGCGATAGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTT‑TACGCCCCCGCGCGCGACC‑AAA‑CAACACCCACAC‑TGGCACCATCTGCACCTCCA‑CATTTGGGGTGCTACGGTGCaaaatggg                                                                                                                                                                        >  1:722575/1‑245 (MQ=255)                                                                                    gTTTCTGAGCGATAGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCGCCTTTTT‑TACCCCCCCGCGCGCGGCC‑CAA‑AAACACCGACAC‑TGGCACCCTCTTCACATCCC‑GATTTGGGGTTCTTCGGTGCAAATCGTGATAGTCg                                                                                                                                                                 >  1:374587/1‑251 (MQ=255)                                                                                          gagcattttttttacccggcccgtgT‑GCA‑CCCGCCATACGGCGCAGCAAAACTGAGGG‑CGAAAACCGCAAAAACCTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGaa                                                                                                                                                            <  2:437789/230‑1 (MQ=255)                                                                                             ggATAGTTTTATCCGTGCCGGGT‑ACA‑CCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACGGTCGCGCCCGAGCCACAGTACAGCGAAAAATTGCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTGC‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGc                                                                                                                                                        <  1:101577/250‑1 (MQ=255)                                                                                             cGATAGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCACTTCGCCTTTTT‑CCCGCCCCCGCGCGCGGCC‑AAA‑CATCACCCACCC‑TGGCACCATCTGCACCTCCA‑GATTTTGGTTTCTACCGTGCAAAACGAGGTAGTCGTGGAACAGc                                                                                                                                                        >  1:8387/1‑251 (MQ=255)                                                                                               aTAGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCTAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCCCCTTTTT‑TACCCCCCGGGGCGCGGAC‑AAA‑CCTCACCG‑CACATGGCACCATCTGCACACCCA‑CATTTGAGTTGCTACAGAGCACAACAAGATAGTCGGGGaaaagatt                                                                                                                                                      >  1:529503/1‑245 (MQ=255)                                                                                                tAGTTTTATCCGTGACGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTCCAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCCCATCTCCTTTTTTAACGCCC‑AGCGCGCGACC‑AAC‑AATAAACG‑CACCTTGCACCATCTGCACAACCA‑AATATAGGCGGTCACGGTGCAAATAGAGCTAGTCGTGGAACAGcat                                                                                                                                                      >  1:132750/1‑248 (MQ=255)                                                                                                 aGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTACTTCGCCTTTTT‑CACGCCCCCGCGCGCGGCC‑CAA‑CATCACCGCCAC‑TGGCACCCACTGCACATCCA‑AATTTGGGTTGCTACGGGGCCAAACGCGATAGTCTTGGGACAGCTT‑a                                                                                                                                                    >  1:156980/1‑250 (MQ=255)                                                                                                   ttttATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCTAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGGTCCCTTTCCCTTTTT‑TACGCCCCGGCCCGCGGCC‑AAA‑CATCACCCCCAC‑TGGCACCCTCTCCACCCCCC‑GAATTTGGTTGCTACGGGGCCAAACGCGGTAGTCGTGGGACAGCtttat                                                                                                                                                   >  1:4354/1‑247 (MQ=255)                                                                                                   ttttATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCCCCTTTTT‑TCCGCCCCCGCGCGCGACC‑AAA‑CATCACCG‑CACCTCGCCCCCTCTGCCCATCCA‑AGTTTTGGTTGCTACGGTTCAAATCGAGATAGTCGTGGAACAACTT‑AAt                                                                                                                                                  >  1:187127/1‑250 (MQ=255)                                                                                                        cccgggccgtgT‑CCA‑TCCGCCATACGGCGCAGCAAAACTGACGG‑CGAGAACCGCAAAAAACTGTCGCGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGTTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGActc                                                                                                                                            <  1:709955/246‑1 (MQ=255)                                                                                                         cccgtgccggtgtaa‑TCCGCCATACGGCCCGGCGAAACTAACGG‑CGAGAAACGCCCAAAACTGACGCGACCTAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGCCC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGActc                                                                                                                                            <  2:292275/237‑1 (MQ=255)                                                                                                          cgtgccgtgT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑TACCCCCCGCCGCCCGGCC‑AAA‑CATCCCCCACAC‑TGGCCCCCTCTTCACAACCA‑GAATTTGGTTGCTCCGGTGCAAAACGAGATAGTCGTGGAACAGCTT‑AATAAACtctt                                                                                                                                          >  1:236380/1‑251 (MQ=255)                                                                                                              ccgtgT‑CCA‑CCCGCCATACGGCGCGGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTGT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCttt                                                                                                                                      <  1:257611/251‑1 (MQ=255)                                                                                                                 tgt‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑TACGCCCCCGCGCGCGACC‑AAA‑CAACACCGACAC‑CGGGACCATCTGCACATCCA‑GATTTTGGTTGCTACGGGGCAAAACGAGAAAGTCGGGGAAAAGCTT‑AATAGACTCTTCttttgg                                                                                                                                   >  1:276080/1‑248 (MQ=255)                                                                                                                 gggttacaccccccATGCGGCGCAG‑GAAATCGACGGGCGAGAACCGCACAATATTGTCGGGCCCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTGAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTc                                                                                                                                     <  2:761695/239‑1 (MQ=255)                                                                                                                      cA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGGGCGGGACC‑AAA‑CATCACCGACAC‑TGGCACCATTTGCACATCCA‑GATTTTGGGTGCTACCGAGCAAAACCAGAAAGTCGTGGAACAGCCC‑AATAAACTCCTCTTTCGGCggg                                                                                                                               >  1:315756/1‑249 (MQ=255)                                                                                                                                 aCGGCGCAGCGAAACGGACGG‑CGAGAAACGCACAAAACGGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCA‑TGTAATAGGTg                                                                                                                     <  1:105821/250‑1 (MQ=255)                                                                                                                                     cgcAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCATTTT‑TACGCCCCCGGGCGCGACC‑AAA‑CATCAACCACCC‑CGGCACCATATGACCAGCCAAGATTTGG‑TTGGTACGGTGCAAATCTGGCTAGTAGAGGAACAAGTT‑AATAAACTCTTATTTCGGCA‑TGGAATAGGGGa‑aaac                                                                                                               >  1:105513/1‑248 (MQ=255)                                                                                                                                             aaCTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCG‑TGTAATAGGTGA‑ACACACGCGGcc                                                                                                       <  1:188203/251‑1 (MQ=255)                                                                                                                                                 gACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCTCCCTTTT‑CACGCCCCGGCGCGCGGCC‑AAA‑CATCACCCACAC‑CGGCACCCTCTCCCCATCCA‑GATTTGGGGTGCTACGGGGCAAATCGCGAAAGGCGTGGGACAGCTT‑AATCGACTCTTCTTTTGGCG‑GGTAATAGGTGA‑ACACACGCcgcccgcc                                                                                                   >  1:697427/1‑251 (MQ=255)                                                                                                                                                         agaAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGATTCACTTCCCCTTTTT‑TACGCCCCCGCGCGCGACCCAAA‑CAAC‑CCCACAC‑TGGCACCCTCTGCACACCCA‑GCTTTGGGGTGCTACCGTGCAAATCCAGGCAGTCGTGGAAGCGCAT‑AATAAAATCTTCTTTCGGCG‑GGTAATAGGGGAGACAC‑CGCGGCCCGCCGagaatg                                                                                            >  1:67425/1‑251 (MQ=255)                                                                                                                                                                             cgACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑TACGCCCCCGCGCGCGACC‑CAA‑CAACACCCACAC‑TGGCACCATCTGCACATCCC‑GATTTTGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCCGCG‑GGGAATAGGTGA‑ACACACGCGGCCCGCCGCGGATGAACACAAAGCGCCGCAgg                                                                          >  1:557552/1‑248 (MQ=255)                                                                                                                                                                                     ccACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑TACGCCCCCGCGCGCGCCC‑AAA‑AAACACCGACAC‑TGGCACCATCTGCACAACCC‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGCGGGACAGCTT‑AATAGACTCTTCTTTCCGCG‑TGGAATAGGTGA‑ACACACGCGGCCCGCCGTGAAAGAAAACATAAATCGGCAGTAGCGTGCCgc                                                                >  1:380089/1‑250 (MQ=255)                                                                                                                                                                                                 gAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTTCCTTCTCCTTTTT‑CTCGCCCCCGCCCGCGCCC‑CAA‑CATCACCCCCAC‑TGGCACCATCTTCACATCCA‑GATTTTGGTTGCTCCGGGGCACATCGAGGTAACCGTGGAAAAGCTT‑TATAGAATTTTTTTTCGACG‑TGACATATGGGC‑ACAAAAGCGGCCCGCCGCGAATGGAAACATAGCGCGCCAGAAGCGTGCCGGCGATTTTCAGCg                                                    >  1:373130/1‑251 (MQ=255)                                                                                                                                                                                                           tGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑TCCCCCCCCGCGCGCGACC‑AAA‑AATCACCCACAC‑TGGGAACCTCTGCACCACCC‑GATTTGGGTTGTTACGGGGCAAAACGAGATAGGCGGAGGACAGCTTTAAATGAC‑CTTCTTTTTGCG‑GGTAATAGGTGA‑GAAAACGCGGCCCGCCCTGAAAGAAAACCTAACACGGGCGTAGCGTGCCgcccgtttctcccggcgctaaagg                                          >  2:506157/1‑228 (MQ=255)                                                                                                                                                                                                                gCCGTTAAGCATACGCAGCGGCGGGTTCAATTCCCCTTTTT‑CACGCCCCCGCGCGCGGCC‑CAA‑CATCACCCACAC‑TTGCACCCTCTGCCCCTCCA‑GATTTGGGTTGCTACGGGGCAAATCGAGATAGACGTGGAAAAGCTT‑AAAAGACTCTTCTTTCGGCGAGGTAA‑AGGTGA‑ACACACGCGGCCCCCCGTGAATGAAAACATAAGGCGGCAGAAGCGGGCCGTCGATTTTCCGCGGCCCTAACggggc                                       >  1:373824/1‑247 (MQ=255)                                                                                                                                                                                                                  cgtaaagcacccccAGCGGCGGTTTCTCTTCGCTTTTTT‑CCCGCCCCGGCGCGCGACC‑AAA‑CACCCCCGACCC‑TGCCCCCATCTGCACATCCA‑GATTTGGGTTGCTCCGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCT‑TGTAATAGGTGA‑CCACCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGc                                   <  1:559344/238‑1 (MQ=255)                                                                                                                                                                                                                    ataggcctccgcggcgtggggTTCACTTCCCCTTTTT‑CCCGCCCCGGC‑CCCCACC‑AAAACACCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAACCGAGAAAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCG‑TGTAATAGGTGA‑CCACACGGGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGc                                   <  1:80206/231‑1 (MQ=255)                                                                                                                                                                                                                       aGCATACGCAGGGGCGGTTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACA‑AAA‑CATCCCCGACAC‑TGGAACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGAAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCG‑TGTAATAGGTGA‑ACACCCGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTc                              <  1:506157/251‑1 (MQ=255)                                                                                                                                                                                                                                     cggGTTCACTTAGCCGTTTT‑CACGCCCCGGCGCGCGACC‑ACA‑CATCACCGACAC‑TGGCACCCTCTGCCCATCAC‑GATTTGGGTTGCTACGGTGCAAATCGAGATCGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCG‑TGTACTAGGTGA‑ACACACGCAGCCCGCAGTGAATGAACACATAGCGCGGCAGTCGGGTGCAGTCGAGTACCTGCGACCATACCGGGTCAGTCAAACCATGCGCCAGACACt                >  2:449182/1‑251 (MQ=255)                                                                                                                                                                                                                                       tgTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCG‑TGTAATAGGTGA‑ACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGGCAGGCAAGTCATGTGCCAGacacggg              >  1:693530/2‑248 (MQ=255)                                                                                                                                                                                                                                        gTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCG‑TGTAATAGGTGA‑ACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCg                                                    >  1:742204/1‑212 (MQ=255)                                                                                                                                                                                                                                                   cTTTTT‑CACGCCACGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCG‑TGTAATAGGTGA‑ACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGACATGTGCAAGACACTGGGCGCGCAACgac  >  2:804514/1‑249 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            ACAGATCCTGACGAGCAGGCAAACGTGGGCCTACGGCAGCCAGACGTTGACGGGCAAAGTGCATACGCGCCACGCGCGCC‑‑AGTTTCTGAGCGATAGTTTTATCCGTGCCGTGT‑TCA‑TCCGCCATACGGCGCAGCGAAACTGACGG‑CGAGAAACGCACAAAACTGTCGCGACCGAGCCACAGTACAGCGAAAAATTTCTGTACGCCGTTAAGCATACGCAGCGGCGGGTTCACTTCGCCTTTTT‑CACGCCCCGGCGCGCGACC‑AAA‑CATCACCGACAC‑TGGCACCATCTGCACATCCA‑GATTTGGGTTGCTACGGTGCAAATCGAGATAGTCGTGGAACAGCTT‑AATAGACTCTTCTTTCGGCG‑TGTAATAGGTGA‑ACACACGCGGCCCGCCGTGAATGAACACATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTC  >  NC_000913/4255818‑4256305

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