Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,660,928	0	T	.	64.1%	0.1 / 25.2	39	coding (515/3114 nt)	mdtF	multidrug efflux pump RND permease MdtF
Reads supporting (aligned to +/- strand):  ref base T (3/11);  new base . (25/0);  total (28/11)
Fisher's exact test for biased strand distribution p-value = 2.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAAT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑TTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCT‑AAACGCCTCGATCATTGTGCAGACGCGTCTGC‑AAACGCCGGAAGAATTTGGCAAAATCCTGTTGAAAGTTCAGCAAGATGGTTCGCAAGTGCTGCTGCGTGATGTCGCTCGCGTCGA  >  NC_000913/3660657‑3661219                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                                      aGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCAC‑GGTTCAAGTGCAAAATAAAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCG‑TAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATCAAAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGAc                                                                                                                                                                                                                                                          >  1:576509/1‑324 (MQ=255)                   gCCTCTATCACTCTGACC‑TCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAAT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCA‑GGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑‑TTATTTCTGATAACGGCAGCCTC‑AACCAGTACG‑TATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCa                                                                                                                                                                                                                                                                            >  1:576513/1‑284 (MQ=255)                                       cGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAAT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATCAAAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑CGG‑TATAACTC‑TTTGGTTCCGAg                                                                                                                                                                                                                                                                                   >  1:576514/1‑264 (MQ=255)                                       cGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAAT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTTATTTCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGt                                                                                                                                                                                                                                                                                  >  1:1478896/1‑262 (MQ=255)                                                      tcCTGATATCGCACAGGTTCAAGTGCAAAAT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑TTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGtgg                                                                                                                                                                                                                                                                                                                     >  1:1478898/1‑219 (MQ=255)                                                                   caGGTTCAAGTGCAAAAT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGAcc                                                                                                                                                                                                                                                         >  1:1478899/1‑259 (MQ=255)                                                                          aaGTGCAAAAT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑TTTATTTCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCggag                                                                                                                                                                                                                                                                                  >  1:1478900/1‑223 (MQ=255)                                                                              tgcaaaT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑‑TTA‑TTCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCT‑A‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TACG‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑c                                                                                                                                                                                                                                                                        >  1:576517/4‑229 (MQ=255)                                                                                aaaaT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATCAAAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑Cg                                                                                                                                                                                                                                                                       >  1:576518/1‑233 (MQ=255)                                                                                     aaaaCTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATctgg                                                                                                                                                                                                                                                                >  1:576519/2‑235 (MQ=255)                                                                                                 cGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑TTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAATTATCAAAA‑GAACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGa                                                                                                                                                                                                                                                           >  1:576522/1‑230 (MQ=255)                                                                                                             aTTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTTATTTCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑AT‑ACCTGGTACCTTcc                                                                                                                                                                                                                 >  1:1478903/1‑251 (MQ=255)                                                                                                                aCCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑TTTATTTCTGATAACGGCAGCCTCAAACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATCAAAA‑G‑ACCCGCTAAAGCCGT‑ACCGCGGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGCAAAAAA‑CTCAAT‑AAA‑T‑ATaa                                                                                                                                                                                                                              <  1:576521/241‑1 (MQ=255)                                                                                                                    gAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑TTTATTTCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑A‑CCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATctgg                                                                                                                                                                                                                                                                >  1:1478905/1‑201 (MQ=255)                                                                                                                     aaGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTTATTTCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTCTTTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑aaact                                                                                                                                                                                                                                                >  1:576528/1‑216 (MQ=255)                                                                                                                             cagcagcagGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCGTTTTATTTCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTCTTTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCg                                                                                                                                                                                                                >  1:1478906/1‑238 (MQ=255)                                                                                                                             cagcagcagGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑TTTATTTCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGc                                                                                                                                                                                                                                                                              >  1:1478907/1‑180 (MQ=255)                                                                                                                                         aTTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑TTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATCAAAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATctggctgg                                                                                                                                                                                                                                                            <  1:576527/188‑1 (MQ=255)                                                                                                                                                                  atatCCTGATGGTAGCGGCG‑TTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATg                                                                                                                                                                                                                                                                          >  1:576529/1‑148 (MQ=255)                                                                                                                                                                  atatCCTGATGGTAGCGGCG‑TTTATTTCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGc                                                                                                                                                                                                                                                                              >  1:1478910/1‑143 (MQ=255)                                                                                                                                                                        tGATGGTAGCGGCG‑TTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTg                                                                                                                                                                                         >  1:1478912/1‑217 (MQ=255)                                                                                                                                                                          aTGGTAGCGGCG‑TTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑CGG‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATctggctgg                                                                                                                                                                                                                                                            >  1:1478913/1‑154 (MQ=255)                                                                                                                                                                                  ggcg‑TTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTcc                                                                                                                                                                                                                 <  1:1478909/184‑1 (MQ=255)                                                                                                                                                                                             tCTGATAAC‑GCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTT‑‑GG‑TACG‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAACCTCAAT‑AAA‑T‑ATTACCTgg                                                                                                                                                                                                                         >  1:1478915/1‑167 (MQ=255)                                                                                                                                                                                                gATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGTAACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑GTTACAACTC‑TTTGGTTCCGAGTAt                                                                                                                                                                                                                                                                                >  1:576531/1‑115 (MQ=255)                                                                                                                                                                                                 aTAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATCAAAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGCAAAACAACCa                                                                                                                                                                              <  1:576530/206‑1 (MQ=255)                                                                                                                                                                                                         aGCCTCGAACCAGTACGATATCGACGGAC‑TACTGTAGCGTCTAA‑TATC‑AAAGG‑ACCCGCT‑AAGCCGT‑ACCGC‑GGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCT‑AAACGCCTCGATCATTGTGCAGACGCGTCTGC‑a                                                                                      <  1:1478914/284‑1 (MQ=255)                                                                                                                                                                                                                        cGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACT‑GGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCT‑AAACGCCTCGATCATTGTGCAGACGCGTCTGC‑AAACGCCGGAAGAATTGGCGAAAATCCTg                                                          >  1:1478917/1‑294 (MQ=255)                                                                                                                                                                                                                        cGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCT‑AAACGCCTCGATCATTGTGCAGACGCGTCTGC‑AAACGCCGGAAGAATTGGCGAAAATCCTGt                                                         >  1:1478918/1‑295 (MQ=255)                                                                                                                                                                                                                           tatCG‑CGGACTTA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTTTGGGTTA‑G‑C‑G‑TACAACTC‑TT‑‑GTTCCGAGTATG‑CATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTA‑TTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCT‑AAACGCCTCGATCATTGTGCAGACGCGTCTGC‑a                                                                                      <  1:576532/264‑1 (MQ=255)                                                                                                                                                                                                                           tatCG‑CGGAC‑TACTGTAGCGTCTAA‑TATC‑‑AA‑G‑ACCCGCT‑AAGCCGT‑ACCGC‑GGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑‑TTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑TAATAACCTGGTACCTTCCGATG‑TA‑TTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCT‑AAACGCCTCGATCATTGTGCAGACGCGTCTGC‑AAACGCCGGAAGAATTTGGC‑AAATCCTGTTGAAAGTTCAGCAAGa                                         <  1:576533/305‑1 (MQ=255)                                                                                                                                                                                                                                           gTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGc                                                                                                                                                                                                                                                                              <  1:1478916/77‑1 (MQ=255)                                                                                                                                                                                                                                                     a‑tatC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑‑AAAA‑CTCAAT‑‑AA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGC‑AAACAACCAGATTTCCGGTGGTCAACT‑GGTGGCATGCCGACAGGGCGGCAGACCAGCAGCT‑AAACGCCTCGATCATTGTGCAGACGCGTCTGC‑AAACGCCGGAAG‑ATTTGGGCAAATCCTGTTGAAAg                                                   >  1:576535/1‑272 (MQ=255)                                                                                                                                                                                                                                                        atC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑‑AAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTcc                                                                                                                                                                                                     >  1:576536/1‑131 (MQ=255)                                                                                                                                                                                                                                                         tC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACcagca                                                                                                                            >  1:1478919/1‑201 (MQ=255)                                                                                                                                                                                                                                                         tC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑‑AAAA‑CTCAATAAAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGAtt                                                                                                                                                                                               >  1:576537/1‑136 (MQ=255)                                                                                                                                                                                                                                                                      gCT‑AAGCCGT‑ACCGCGGGCG‑T‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTAcc                                                                                                                                                                                                                     >  1:1478923/1‑104 (MQ=255)                                                                                                                                                                                                                                                                        t‑AAGCCGT‑ACCGTCGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TT‑‑GTTCCGAGTATGCCATG‑CGTATCTGGCTGGA‑CCGC‑‑AAAA‑CTCAAT‑AAATT‑ATAACCTGGTACCTTCCGATGTTA‑TTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCTAAAACGCCTCGATCATTGTGCAGACGCGTCTGC‑AAACGCCGGAAGAATTTGGCAAAATCCTGTTGAAAGTTCAGCAAGATGGTTCGCAAGTGCTGCTGCGTGATGTCGCTCGCGTCGa  <  1:576534/307‑1 (MQ=255)                                                                                                                                                                                                                                                                                         gCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCTAAAACGCCTCGATCATTGTGCAGACGCGTCTGC‑aaaac                                                                                  <  1:1478920/216‑3 (MQ=255)                                                                                                                                                                                                                                                                                           ttt‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTA‑TTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCTAAAACGCCTCGATCATTGTGCAGACGCGTCTGCAAAACGCCGGAAGAATTTGGCaa                                                                 <  1:1478922/230‑1 (MQ=255)                                                                                                                                                                                                                                                                                           gTT‑‑GG‑TA‑GCC‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATGCCGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTgg                                                                                                                                                          <  1:1478921/147‑1 (MQ=255)                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                                      AGTGATGCGGCGGGCAATGCCTCTATCACTCTGACCTTCGAGACTGGGACATCTCCTGATATCGCACAGGTTCAAGTGCAAAAT‑AAACTGCAACTCGCTATGCCTTCATTACCTGAAGCAGTGCAGCAGCAGGGGATTAGCGTCGATAAGTCGAGCAGTAATATCCTGATGGTAGCGGCG‑TTTATTTCTGATAACGGCAGCCTC‑AACCAGTACGATATCG‑CGGAC‑TA‑TGTAGCGTCTAA‑TATC‑AAA‑G‑ACCCGCT‑AAGCCGT‑ACCGCGGGCGTT‑‑GG‑TA‑G‑C‑G‑TACAACTC‑TTTGGTTCCGAGTATGCCATG‑CGTATCTGGCTGGACCCGC‑AAAAA‑CTCAAT‑AAA‑T‑ATAACCTGGTACCTTCCGATGTTATTTCCCAGATTAAGGTGCAAAACAACCAGATTTCCGGTGGTCAACTGGGTGGCATGCC‑ACA‑GGCGGCAGACCAGCAGCT‑AAACGCCTCGATCATTGTGCAGACGCGTCTGC‑AAACGCCGGAAGAATTTGGCAAAATCCTGTTGAAAGTTCAGCAAGATGGTTCGCAAGTGCTGCTGCGTGATGTCGCTCGCGTCGA  >  NC_000913/3660657‑3661219

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