Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,652,331	T→C	intergenic (+2346/+397)	ydfD → / ← ynfP	lysis protein/protein YnfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,652,331	0	T	C	100.0%	132.1 / NA	40	intergenic (+2346/+397)	ydfD/ynfP	lysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/18);  total (22/18)

TATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAA  >  NC_000913/1652103‑1652559                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       tataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATAGCAGTGATGTTCACCTCAGAGAGTATGGcagaacag                                                                                                                                                                                                                   <  1:151505/250‑1 (MQ=255)                      gAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGAGTTTCACCTCAGAGAGTATGGCAGAACATAAAAACAAAAATG‑CACCTTTGt                                                                                                                                                                                            <  2:56439/251‑1 (MQ=255)                               gCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTc                                                                                                                                                                                       <  2:139459/247‑1 (MQ=255)                               gCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGGATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGGTTTTc                                                                                                                                                                                       >  1:139459/1‑247 (MQ=255)                                cacgtGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCGCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CTCCTTTGTTTTTCAACCt                                                                                                                                                                                  <  2:293126/247‑1 (MQ=255)                                    tGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCAATCATCTTAATAAGACAGCAAACAACGCTCTACGGAGGCAGATGAAACTTAA‑GAGAATGTGCAAAAAGCATAACACTGATGGTCACCGCAGAGAGAATGGCAGAACAGAAAAACAAAAAAG‑CACCTTTGTTTTTCAATCTgggg                                                                                                                                                                              >  2:229866/1‑248 (MQ=255)                                     ggATGGGAGGACATAGATCTGAAAAATGGAACACTGAAGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGAAACCCATCGGGTCCTTCATTTTATTAAGCCAGCAATCCACGCTCTCCCGAGTCAGATGACA‑TTAACGAGACTGGGCAAAGAGCCTATCATTGATGTTCACCTCAGGGAGTATGGAAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGa                                                                                                                                                                                >  1:407692/1‑248 (MQ=255)                                          ggaggaCATAGATCGGAAAAATGTAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATGAAGCCAGCAAGCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCAGTGATGTGCACCTCAGAGAGTATGGCAGTACAGAAAAACAAAAAGG‑CACCTTTGGTTTTCAACCTGAAGTGTCAg                                                                                                                                                                        <  2:396559/251‑1 (MQ=255)                                               aCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAACCGGGTCAGTCATCTTATTAAGCCAACCATCGAAGCTATCCGGAGACCAAAGAAA‑TTAACGAGACTGAGCAAAGAGCATATCATTGAGGATCACCTCAGAGAGTCTGGCAGAACAGAAAAAACAAAATGTCAC‑TTTGTTTTTCAACCTGAAgtggaagcg                                                                                                                                                                      >  1:316257/1‑242 (MQ=255)                                                cATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGagag                                                                                                                                                                  <  1:36855/251‑1 (MQ=255)                                                 aTAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGt                                                                                                                                                                 <  1:232754/251‑1 (MQ=255)                                                 aTAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGCTGACA‑TTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGGGTCCACGagagg                                                                                                                                                                 >  2:412065/1‑250 (MQ=255)                                                  tAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACCCATAGGGTCATTTATCCTATTTAGCCCGCAATCCACGCCCCCCCGAGTCCGACGACC‑TTAGCGAGACTGAGCAAAAAGCCTTTATTTGATGTTCACCCCCGAGAGAATGGCAGAACAAAAAAACAAAAATA‑CACCTTTGTTTTTCAACCTGAAATGTCAGCgagcga                                                                                                                                                                 >  2:111299/1‑247 (MQ=255)                                                         gAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaaaa                                                                                                                                                            <  1:412065/248‑1 (MQ=255)                                                          aaaaaTGGAACAATGATGATCAGCAGAAATTTAGCCAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGCACCAATCCGGTAATTCATCCTATTTAGCCAGCAATCGACGCTCTCCGGAGACAGATGACC‑TTAACGAGACTGAGCAAAGAACATATCATTGATGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAAAG‑CACCTTTGTTTTTCAACCTGAAGAGTAAGCGAGAGTAAAAAAtt                                                                                                                                                         >  1:401790/1‑250 (MQ=255)                                                                              cAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCTACGCTCTCCGTAGTCAGATGACA‑GTTACGATTCTGAGTAAAGAGCGTATCATTGATGTTCAGCTCAGAGTTTATGGCATAACAGAAAAACAAAATGG‑CACCTTTGGTTTTCAATCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGtt                                                                                                                                    <  1:287884/250‑1 (MQ=255)                                                                              cAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACAAATCGGGGCATTAATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCCGAAGACA‑TTAACGAGACTGAGCAAAGAGCCAATAATTGGTGTTCACCTCAGAGAGGATGGCAGAACAAGAAAACAAAAATG‑CACCTTTGTTTTTCAAACTGAAATGACAGCGAGAGAAAAAAAATATGGTGACCA‑TTTTACCGttt                                                                                                                                   >  2:90163/1‑250 (MQ=255)                                                                                 gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCCCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCCGCGAGAGTAAAAAATTATGGTGAACA‑TTTTTCCGTTGACt                                                                                                                                >  1:288455/1‑251 (MQ=255)                                                                                      attgaGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATa                                                                                                                           <  2:315653/247‑1 (MQ=255)                                                                                            cAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGTAACCAAACAGGTCATTCCACTTATTTAGCCCGCAATCGACGGCCTCCGGAGTCCGAAGGCA‑ATAACGAGACTGGGCAAAGAGCATATCATTGATGTTCACCTCAGAGGGAATGGCAGAACAGAAAAAAAAAAATT‑CACCCTTGTTTTTCAACCCGCAATGTCAGCGAGAATTAAAAAATATGGTTACCAATTTTAcgttgtatcaattaggg                                                                                                                        >  2:75833/1‑233 (MQ=255)                                                                                             aaaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTTATCTTATTAAGCCCGCAATCGACGCTCTCCGGAGACAGATGAAA‑TTAACTAGACTGAGCAAAAAGCATATCATTTATGTTTACCTCAGAGAGGAAGGCATAAAAGAAAAACCAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGGAAAAAATTATGGGGACCA‑TTTTTCCCTTGGATCAATAAGGCAgg                                                                                                                    >  1:250799/1‑250 (MQ=255)                                                                                                aGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATTGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGAtgt                                                                                                                 <  1:451518/251‑1 (MQ=255)                                                                                                    cGTTTCACGGTACCAAAAACACAGGCGGGAACCAAAAGGATCATTCAAACTACAACACCAGCCATCGCAGCTTTCCGGGGTCAAATGACA‑TTAACGAAACTCAGAAAAGAGCACATCATTGATGTTCACCTCAGAGAATATGGCAGAAAAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAACA‑TTTTACCGTTGAATCAATAAGGCAGAttttggg                                                                                                             >  1:224278/1‑245 (MQ=255)                                                                                                         cACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCCTTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTCAGCCCAGAGCATATCATTTATGTTCACCTCCGAGCGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGCGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATTGGGGACGccc                                                                                                        >  1:34559/1‑249 (MQ=255)                                                                                                                ccAAAAACACAGGCGGGGAACAATAGGGTCATTCATCTTATTAAACCAACAAACGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTTAGCAAAAAGCATAACATTGATGGTTACCCCCGAGAGGATGGCAGAAAAAAAAAAAAAAAATG‑CACCTTTTTTTTTCAACCCGAAGTGTCAGCGAAGGTAAAAAATTTTGGGTAACA‑TTTTACCGTTTACTCATAAAAGCAGATGTGGGACACAGCGATaa                                                                                                  >  2:174153/1‑250 (MQ=255)                                                                                                                     aaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTg                                                                                            <  2:141135/251‑1 (MQ=255)                                                                                                                         caGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CTCCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCgg                                                                                        <  2:147661/251‑1 (MQ=255)                                                                                                                           ggCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGAATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACAGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGAc                                                                                      >  2:350477/1‑251 (MQ=255)                                                                                                                                gAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGCAGTCAGATGACA‑TTAACGAGCCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTATCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGc                                                                                 >  1:418709/1‑251 (MQ=255)                                                                                                                                gAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATAGCCGCTCTCCGGAGTCAGATGACA‑TTAACGAAACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGCAAAACAACAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCC‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGc                                                                                 >  1:413075/1‑251 (MQ=255)                                                                                                                                 aaccaaccGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGcc                                                                                >  1:126823/1‑251 (MQ=255)                                                                                                                                                                      gACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctg                                           <  1:62787/251‑1 (MQ=255)                                                                                                                                                                        cGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGTCCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGG‑ACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgctg                                        <  2:396881/251‑1 (MQ=255)                                                                                                                                                                                     cAGATGAAA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGCTCACCTCAGAGAGTATGGCAGACCAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGCCTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGCCGCTGACAGc                            >  2:194363/1‑251 (MQ=255)                                                                                                                                                                                                    gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAg                             <  2:118817/236‑1 (MQ=255)                                                                                                                                                                                                    gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAg                             >  1:118817/1‑236 (MQ=255)                                                                                                                                                                                                        cTGAGCAAAGAGCATATCATTGATGTTCACCTCAGCGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCCGCGAGAGTAAAAAATTATGGTGACCA‑TTTTAACGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGAGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCAtt          >  1:222814/1‑251 (MQ=255)                                                                                                                                                                                                             cAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAg     <  1:350477/251‑1 (MQ=255)                                                                                                                                                                                                               aaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCa   <  2:288455/251‑1 (MQ=255)                                                                                                                                                                                                                 gagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAACACGTGCCGGACTCCGCCATCGAAACTCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCaa  >  2:210463/1‑250 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       TATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATG‑CACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCA‑TTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAA  >  NC_000913/1652103‑1652559

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