Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,411,689	G→T	E13* (GAA→TAA)	yhdJ →	DNA adenine methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,411,689	0	G	T	97.5%	123.3 / NA	40	E13* (GAA→TAA)	yhdJ	DNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (17/22);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAA  >  NC_000913/3411442‑3411919                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       gTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTTGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATGCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCTGCATGGGGAAGCGCCTTGTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAg                                                                                                                                                                                                                                     <  1:33940/251‑1 (MQ=255)                  ggACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCa                                                                                                                                                                                                                     <  1:337950/250‑1 (MQ=255)                                      cccGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCCTGAGAACCGGATGTGAACCGACCCGGTTTGGTAATTAAGCTACGAcc                                                                                                                                                                                                                             >  2:491031/1‑222 (MQ=255)                                      cccGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGAcc                                                                                                                                                                                                                             <  1:491031/222‑1 (MQ=255)                                         gTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTaa                                                                                                                                                                                             <  2:67875/251‑1 (MQ=255)                                                               gcTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCAGGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAGGTGAACCTACCCGGTTTGGTAATTAAGCTAAGACCATTATGCACGGTAATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTc                                                                                                                                                                       <  1:182674/251‑1 (MQ=255)                                                                  gatgatGGGCATCACCCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCgat                                                                                                                                                                    <  2:457614/251‑1 (MQ=255)                                                                               agccggggTACGTTGGGTAAAAAATGGAAAAAATACGGCATGACTTAATTCAGGTTAGCAAAATGCTGGATTAAAAAGGCGCCCCGCGGCAGGGGGAACCGCTTTTTTTATAGGTGTCACAAAGGGAGTGCCCATGAGACAATGACGGAAATCGACCCGGTTGGGTAATGAAGCTAAGGCTATTATTCCCGGTGATGCGCTGGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCg                                                                                                                                                       <  1:548104/245‑1 (MQ=255)                                                                               aaCCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCg                                                                                                                                                       >  2:301327/1‑251 (MQ=255)                                                                                  cGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAACGTGTCGATCTGATCTTTGCCGAcc                                                                                                                                                    >  1:219854/1‑251 (MQ=255)                                                                                    tGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAACGTGTCGATCTGATCTTTGCCGAcccc                                                                                                                                                  >  2:16601/1‑250 (MQ=255)                                                                                       tACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACccacc                                                                                                                                                >  2:431957/1‑250 (MQ=255)                                                                                          gCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCatat                                                                                                                                            >  2:207378/1‑251 (MQ=255)                                                                                             gCGTAAAAAATTGAAAAAATCCGGCATGAACGAATCCAGGTTAGCTAAATGTTGGATGAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGTTGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTGTGCCGACCCACCATATAAc                                                                                                                                         <  2:405840/251‑1 (MQ=255)                                                                                                    aaaTTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGCCCATTATTCACGGTGATGCGCTTGCCTATCTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTa                                                                                                                                  <  1:191630/251‑1 (MQ=255)                                                                                                    aaaTTGAAAAAATACGGCATGAACAAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAGAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTa                                                                                                                                  <  1:431957/251‑1 (MQ=255)                                                                                                          aaaaaaTACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAAtt                                                                                                                            >  2:136030/1‑251 (MQ=255)                                                                                                                     cATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCt                                                                                                                                                                                                                                   >  1:176980/1‑137 (MQ=255)                                                                                                                     cATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCt                                                                                                                                                                                                                                   <  2:176980/137‑1 (MQ=255)                                                                                                                      attaaCTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGa                                                                                                                <  2:471552/248‑1 (MQ=255)                                                                                                                        gAACGAATTCAGGTTAGCTAAATGCTGGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGAt                                                                                                               <  2:92272/250‑1 (MQ=255)                                                                                                                             aatccAGGTTAGCTAAATGCTTGAATAAAAAGGCGCTACTCGGCATGGGGAAGCGCTTTTTTTATAGGTGTCACAAAGGGATTGACCATGAGAACAGGATGTGATCCGACCCGGTTTGGTAATTAAGCTATGACCATTATTCTCGTTGATGCTCTTGCTGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAg                                                                                                          <  1:273910/246‑1 (MQ=255)                                                                                                                             aaTTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACGCGGCATGGGGAAGCGCTTTTTTAATAGGTGTCAAAAAGGGAGTGACCATGAAAACAGGAGGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCg                                                                                                             <  2:113950/247‑1 (MQ=255)                                                                                                                                 cAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCAGGGGGAAGCTCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTgg                                                                                                     <  2:512601/251‑1 (MQ=255)                                                                                                                                       aGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAaagaa                                                                                               <  1:207378/251‑1 (MQ=255)                                                                                                                                         caaaaTGCTTAATTAAAAAGGGGCTACGAGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGGTAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAAGTTTGATGGTCTGATCGAAGCCTGGAAagaaga                                                                                             <  2:27020/249‑1 (MQ=255)                                                                                                                                              tGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGt                                                                                        >  2:136093/1‑251 (MQ=255)                                                                                                                                                            ggCGCTACTCAGCATGGGGACGCGCCTTTTTTATAGGAGTCCAAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGATAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGCTCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAACATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTg                                                                          >  2:516090/1‑251 (MQ=255)                                                                                                                                                            ggCGCGACTAGGCCTGGGGTAGCGCCTTTTTTCTAGGTGTACCACGGGGAGTGAACATGAGAAAAGGATGTGGACCGGCCCGGTTTGGAAATTAAGCTAAGACCAATATTCACGGAGATGCGCTTGCCGCACTTAAAAAGATACCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATCACATCGGTAAAAATTTTGATGGTCTGATCCACGACTGGAAAGAAGATCTGTTTCTCGACTGGCtt                                                                          >  2:539118/1‑250 (MQ=255)                                                                                                                                                               gcTACTCGGCATGGGGATGCGCCTTTTTTATGGGTGTCAAAAAGGGAGTGACCATGAGAACTGGATTGGAACGGTCCCGGTGTGGTAATTAAGCTAAGACCATTATTCACGTTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCAGATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGAGCTGTTTATCGACTGGCTGttt                                                                       <  1:539118/251‑1 (MQ=255)                                                                                                                                                                     cGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGAGGTGAACAGCCCCGTTTTGGAAATTAAGCTAAGACCATTATTCACGGTGATGCGCTCGCCGAACATAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGAGGGTCTGATCGAAGCCTGGAAAGAAGATCTGGTTATCGACTGGCTGTTTGAAg                                                                   <  2:96583/249‑1 (MQ=255)                                                                                                                                                                             ggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGc                                                            >  1:157064/1‑248 (MQ=255)                                                                                                                                                                                              aGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTg                                         <  1:136030/250‑1 (MQ=255)                                                                                                                                                                                                            gAGTGACCATGAGAAAAGGATGTAAACCGAACCGGTCTGGCAATTAAGCTAAGAACATTATTCACGGCGATGAGATTGCCGAACTTAAAAAGATCCACGCCGAAAGCGTCGATCTGCTCTTTGCCGACCCACCATATAACATCGATAACAATTTTGATGGTCTGCTCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGCTTGCAGAGTGCCACCGAGTTCTGAAAAAGCAGGgcagc                          >  2:347734/1‑251 (MQ=255)                                                                                                                                                                                                               tGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATg                       >  2:158453/1‑251 (MQ=255)                                                                                                                                                                                                                                 gtgagccgaacaggTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATACCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCAGATAACATCGGGAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGACGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACg     >  2:12786/13‑251 (MQ=255)                                                                                                                                                                                                                                 gtgAACCGACCCGTTTTGCTAATTACGCTAAGCCCCTTATTCTCTGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAATTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTAc      <  1:155736/250‑1 (MQ=255)                                                                                                                                                                                                                                     aCCGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGaa  >  1:55626/1‑250 (MQ=255)                                                                                                                                                                                                                                       cGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATc                                                                                                              <  2:189112/140‑1 (MQ=255)                                                                                                                                                                                                                                       cGACCCGGTTTGGTAATTAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATc                                                                                                              >  1:189112/1‑140 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                       GTAGAACAGCCCCTGTTGGACATGGTGATGCAATACACCCGTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAGATCTGTTTATCGACTGGCTGTTTGAAGTGATTGCAGAGTGCCACCGCGTTCTGAAAAAGCAGGGCAGCATGTACATCATGAACAGTACGGAA  >  NC_000913/3411442‑3411919

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