Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	623,980	A→T	18.9%	D177E (GAT→GAA)	fepB ←	ferrienterobactin ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	623,980	0	A	T	18.9%	109.9 / 2.9	49	D177E (GAT→GAA)	fepB	ferrienterobactin ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (11/27);  new base T (8/2);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 8.18e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCAGCCAGATTTTCCCCACCAAGCTGAATGATGTCATGGCGTTTACCCTGACTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCC  >  NC_000913/623745‑624201                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     cTGGAGCCAGTTTTTCCCCACCAAGCTGAAGGATTTCAGGGCGGTTACCCTGACTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCa                                                                                                                                                                                                                                   <  1:542710/239‑1 (MQ=255)        ctgtttttCCCCGCCAAGCTGAATGGTTCGAGTGGTTTTACCCTGACTTTGGCTGGTATTTACCCCTGCGGGTAACTTCGTCAGCGTAAAGCCGAGTTCTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCCGTGAT‑TGGCGCTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑Tgcg                                                                                                                                                                                                                           <  1:808795/235‑1 (MQ=255)               tCCCCACCAAGCTGAATGATGTCATGGCGTTTACCCTGACTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTTGTTTTTTCTTTTGCCGCCGCCCGTTGCTTTTCAAAC‑TTCGCA‑AT‑cc                                                                                                                                                                                                                   >  2:370717/1‑238 (MQ=255)                    cccaCGCTGAATTTTTTCATGGCGTTTTCCCTGCTTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTg                                                                                                                                                                                                             <  1:370717/238‑1 (MQ=255)                           ttaaTGGTGTCTTGGCGTTTACCCTTACTTTGGCTGGCATTTAACCCGGCGGGTAACTTCGCCAGCTTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGc                                                                                                                                                                                                         <  1:215098/234‑1 (MQ=255)                            gAATGATGTCATGGGTTTTACCCTGGTTTTGGCTGGCAATTAACCCTGCGGGTAACTTCTCCAGCGTAAAGCCGAGTTGTGCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAG‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGt                                                                                                                                                                                                     <  1:382410/239‑1 (MQ=255)                                  tgTCATGGCGTTTACCCTGACTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGGGACCGGCTGCGGCGGT‑AATTTGTTTTGTTCTTTTGCCGCCCCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCCTGTTtttc                                                                                                                                                                                                >  2:863424/1‑238 (MQ=255)                                    tttttggtttttCCCTGTTTTTGGTGGGCATTTAACCCCGGGGGTATCTTCGCCCGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATg                                                                                                                                                                                             <  1:48356/231‑1 (MQ=255)                                    ttctGGTGTTTTCCCGGTCTTTGGCCGGTATTTAACCCTGGGGGTTAGTTCGGCAGCGTAAAGCCGGGTTTTTCCAGCATCCGGCCTTGTGCTGATTCTGGCCGCCAGAGAT‑TGGCCCTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATg                                                                                                                                                                                             <  1:847071/236‑1 (MQ=255)                                       tGGCGTTTACCCTGACTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGG‑AATTTGATTTGCTCTTTTGCCGCCCCCCGGTGCTTTTCCAAC‑CGCGGC‑AT‑CCGCCCTGCCGCTCTTGTCTCATT‑ccc                                                                                                                                                                                         >  2:27810/1‑239 (MQ=255)                                       tGGCGTTTACCCTGACTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGG‑AAATTGGTTTTCCCTTTTGCCGCCCCCCGGTGCCTTTTCAAC‑CGCGCA‑ATCCCGCCTTCCCCCTTGTTtctcttccc                                                                                                                                                                                           >  2:175243/1‑234 (MQ=255)                                       tGGCGTTTACCCTGACTTTGGCTGGCATTTAACCATGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCCGCGGCGGTGAATTGGATTTGCCCTTCCGCCACCGCAAGTGGTTTTGCAAAC‑GGAGCA‑AC‑CCGccagccccgttgttttcactccc                                                                                                                                                                                           >  1:91194/1‑217 (MQ=255)                                            tttACCCTTTTTTTGGCTGGCTTCGAACCTGGCGTGTAACTCGGCTAGGGTAAAGCCCGTTTGGTCCGGTTTCTGCCCTTGTGCTGATTCTTGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGTATTACAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTTCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCGGCCGCTTGTTTCTCATG‑CCCGGTaa                                                                                                                                                                                    <  1:765963/239‑1 (MQ=255)                                            tttACCCTGACTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTTATTTGCTCTTTTGCCGCCGCCAGTTGCTTATTAAAA‑CGCGCA‑AA‑CCGCCCTGCCCCCTGTTTTCCATT‑ccccggga                                                                                                                                                                                    >  2:740399/1‑234 (MQ=255)                                                      cTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGGGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCaa                                                                                                                                                                           <  2:352374/238‑1 (MQ=255)                                                         tggctggcATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGCTTTGCTCTTTCGCCACCCCCCGCTGATTTTCAAAA‑TGCGCA‑AA‑CCCCCCCTCCGCGTGTTTTTCATGCCCCGGAAATTTcccacgttt                                                                                                                                                                        >  2:491716/1‑231 (MQ=255)                                                          ggctggcATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGG‑AATTTTGTTTGCTCTTTCGCCGCCCCCCGATGCTTTTTCAAC‑AGGGCACAT‑CC‑CTCTGCCCCTTTTTTTTCCTG‑cccccgcattttcccaagtatc                                                                                                                                                                      >  1:536210/1‑220 (MQ=255)                                                             tgggttttAACCTTGCGGGTAACTTCGCTGGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGTTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGtt                                                                                                                                                                   <  1:740399/234‑1 (MQ=255)                                                                     aaCCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCCGTTGCTTATCCAAC‑TGCGCA‑AT‑CCCCTCTGCCGCTTGTTTTTCATG‑CCCGGTAAATTCGCCCAGTTTCGTTTA‑CAACg                                                                                                                                                           >  2:708477/1‑238 (MQ=255)                                                                     aaCCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCCGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTAGTTTGCTCTTTCTCCGCCGCCAGTTGCTTTTTCAAC‑TCCGCA‑AA‑CCCCTCTGCCCCTTTTTTTTCTTG‑CCCCGGAAATTTGCCCAGTTTGGTTaaaaccg                                                                                                                                                            >  1:703459/1‑233 (MQ=255)                                                                          tGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGG‑AATTTTATTTGCTCTTTTGCCCCCGCCCGTTGCTTCTCAAAA‑TCGGCC‑AT‑CCGCCCTGCCCCTCGTTTTTCCTG‑ACCGGGACTTTTGCCACGATGCGTTTA‑CACCACCTGcc                                                                                                                                                     >  1:190581/1‑239 (MQ=255)                                                                                              tAACGCCGGGTTGTCCCAGCTTCTGCCTTTTTGCTGATTCTGGCGTCAAGAGAT‑TGGCACTGGTGGAAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATtg                                                                                                                                <  1:708477/239‑1 (MQ=255)                                                                                                                     tGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTTCTCTTTCGCCGCCGCCAGTTGCTTATTCAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTTCTG‑CCCGGGAATTTCGCCCAGATGCGTTAA‑AAGCGCCCGCCACCCTTTTTCGT‑CGTAATTGTTGGTTAATGtttgggggcgggg                                                                                                         >  2:89440/1‑227 (MQ=255)                                                                                                                      gCCTTTGTGTTGTTTCTGGCGTCCAGAGCTGGGGTATTGTG‑GCAGTGGCAGTATGGACAATGGCAGTGCCCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGGCTGTTGCTTATCAAAC‑TGCGCA‑AG‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGAtggtg                                                                                                        <  1:491716/239‑1 (MQ=255)                                                                                                                        ccTTGTGTTTGTTCTGGGGTTCAGGGAT‑TGGCACTGGTTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCGGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGa                                                                                                      <  1:646753/239‑1 (MQ=255)                                                                                                                          ttggggtgtttttggggttccgggaT‑TGGCACTGTGTGCCGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCTGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGaaa                                                                                                    <  1:819816/216‑1 (MQ=255)                                                                                                                          tgtggCTGTTTTTGGCGTCCAGAGTT‑TGGCCCGGTGTGCAGGGGCAGTATAGACAATGGCAGGGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGaaa                                                                                                    <  1:293864/235‑1 (MQ=255)                                                                                                                          ctgttgtttttCTGGGGGCCAGAGAT‑TGGCATTGTGTGCAGTGGGAGTATAGACAATGTCAGTGCCCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGaaa                                                                                                    <  1:878218/230‑1 (MQ=255)                                                                                                                           tgtgCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTTATTTTCTCTTTCGCCGCCGCCCGTTGCTTTTCAAAC‑TGCCCC‑AT‑CCCCCCTCCCCCTTTTTTCTCATG‑CCCGGGAATTTCTCCCAATTTGGTTAA‑ACACGACTTCCCCCTTTTTTCCT‑CGTAAGTGTAGTTTAAtatacgggggagggggggaaa                                                                                                   >  1:460524/1‑219 (MQ=255)                                                                                                                                 gATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTTGCCCAATTGCGTTTA‑AAACGGATGCCAGCCTTTTTTGT‑CGTAATTGTTGATTAATTTCGGGGCGATGGGGGAAAGCT‑GGTc                                                                                            >  2:646753/1‑239 (MQ=255)                                                                                                                                     cTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATAc                                                                                        <  2:447546/239‑1 (MQ=255)                                                                                                                                     cTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAAGGGCAGTGCCCGGCTGCGGCGGT‑AATTTGATTTGCTCTTGCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATAc                                                                                        <  2:536210/239‑1 (MQ=255)                                                                                                                                                             tgtgCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTTCTCTTTCGCCGCCGCCCGCTGCTTATCAAAC‑TGCGCA‑CA‑CCGCTCTGCCGCTTTTTTCTTATG‑CCCGGGAAATTTTACCAGTTGCGTTAA‑CACCAACTGCCCCACTTTGTTGTGCGTATTTG‑TGATTTATGTCGGGGGGGTGGTGGGAAACTAGATTATACataccagcacccgacacccccc                                                                  >  1:718186/1‑218 (MQ=255)                                                                                                                                                                     ggCAGTATAGCCAATGGCAGGTACCGGCTTGGGCGGG‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTgcg                                                          <  2:499478/238‑1 (MQ=255)                                                                                                                                                                     ggCAGTATAGACAAGGGCAGTGCCCGGCTGGGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAATC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTCCGCCAAGTGGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTgcg                                                          <  2:754492/238‑1 (MQ=255)                                                                                                                                                                       cAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTTCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGGAAGCT‑GGTCATACAGTGCCAGCGCCGAATCCCCCCCCGTTGCGCtc                                                       >  2:571526/1‑238 (MQ=255)                                                                                                                                                                          tataGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGa                                                                                                             <  2:468339/183‑1 (MQ=255)                                                                                                                                                                                cAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCCCTTTCGCCGCCGCCGGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAAt                                               <  2:119971/238‑1 (MQ=255)                                                                                                                                                                                 aagggCATTGGCCGGTTGGGGGGGT‑ATTTTTATTTGCTTTTTCGCCCCCGCCAGTTGTTTTTAAAAC‑TGCCCA‑AT‑CCGCTTTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCTAAGTTGTGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCa                                             <  1:518365/236‑1 (MQ=255)                                                                                                                                                                                      cggtgcgcgggggggggcgT‑AATTTTGATTGCCCTTTCCCCTCCGCCTGTTTGTTTACAAACGTGGACA‑GT‑CGGCTC‑GCCGGTTTTTTCTCAAT‑GCCGGTTATTTCCCCAAGTTTGGTTAACCGGGGATGGCCGGCTTTTGCGGT‑CGTAATTGAGGTTTAAGGCGGGGGGGATGGT‑GAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATc                                         <  1:175243/220‑1 (MQ=255)                                                                                                                                                                                                                    tgCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GACCATCCAGGGCCAGACCCGAAACCCCCCCGGTGGCGGTAATTAAAAAAAGAACCGGCCTTGGGCCGCCAACGCCtcgcgagcc           >  1:852912/1‑231 (MQ=255)                                                                                                                                                                                                                      ctctTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATAATACAGTGCCAGCGCCGAATCCCCGCCGGTTGAGCTAATTAAAATCCGATCCTGCATTTGCGCGGCAACGGCTTCGGCGCtgg         >  2:501046/1‑237 (MQ=255)                                                                                                                                                                                                                      ctctTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGCTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATAATACAGTGCCATCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCCGGTCCGCCATTTGCGCGGCAACGGCTTCGTAGCTCg         >  2:230839/1‑239 (MQ=255)                                                                                                                                                                                                                      ctctTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCAGG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGTGTGCGCTAATTAAAATCAGATCTGGCATTTGCAGGGCAACgccctccggggtgg         >  2:240073/1‑226 (MQ=255)                                                                                                                                                                                                                         tttCGCCGCCCCCATTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGtt      <  2:574133/239‑1 (MQ=255)                                                                                                                                                                                                                         tttCGCCCCCGCCAGTTGTTTACAAAAC‑TGCGCA‑AT‑CCGCTCTGCCGTTTGTTTCTATGG‑CCGGGTAATTTCGCCAAGTTGGGTTAA‑CAGCGCCTGCCAGCTTTTGTCGG‑CGTAATTGATGATTAATGGCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGtt      <  1:571526/239‑1 (MQ=255)                                                                                                                                                                                                                            cgccgccgcCATTTGTTTAACAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGc   <  2:703459/239‑1 (MQ=255)                                                                                                                                                                                                                            cccccccgccATTTGTTTATCAACC‑TGCGCA‑AC‑CCCCTTTGCCGGTTTTTTTCTATG‑CCCGGTAATTTCGCAAAGTTGCGTTAA‑CGGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGc   <  1:863424/234‑1 (MQ=255)                                                                                                                                                                                                                             gccgccgccAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGcc  <  2:340690/239‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                     CTGCAGCCAGATTTTCCCCACCAAGCTGAATGATGTCATGGCGTTTACCCTGACTTTGGCTGGCATTTAACCCTGCGGGTAACTTCGCCAGCGTAAAGCCGAGTTGTTCCAGCATCTGCCCTTGTGCTGATTCTGGCGTCCAGAGAT‑TGGCACTGTGTGCAGCGGCAGTATAGACAATGGCAGTGACCGGCTGCGGCGGT‑AATTTGATTTGCTCTTTCGCCGCCGCCAGTTGCTTATCAAAC‑TGCGCA‑AT‑CCGCTCTGCCGCTTGTTTCTCATG‑CCCGGTAATTTCGCCAAGTTGCGTTAA‑CAGCGACTGCCAGCTTTTGTCGT‑CGTAATTGATGATTAATGTCGGGGCGATGGTGGAAAGCT‑GATCATACAGTGCCAGCGCCGAATCCCCGCCGGTTGCGCTAATTAAAATCAGATCCGGCATTTGCGCGGCAACGGCTTCGGCGCTCGGTTCGCC  >  NC_000913/623745‑624201

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