Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,219,752	C→T	S386F (TCC→TTC)	aceK →	isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,219,752	0	C	T	97.5%	112.5 / ‑7.0	40	S386F (TCC→TTC)	aceK	isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/23);  minor base A (1/0);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

TATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTC  >  NC_000913/4219505‑4219997                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      tATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTtccc                                                                                                                                                                                                                                                    >  1:311089/1‑251 (MQ=255)   tCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGCTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGTTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACTCCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCg                                                                                                                                                                                                                                                   <  2:121740/250‑1 (MQ=255)              gTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGa                                                                                                                                                                                                                                        <  1:6927/250‑1 (MQ=255)              gTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGa                                                                                                                                                                                                                                        <  2:576218/250‑1 (MQ=255)                   gAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCCCGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACt                                                                                                                                                                                                                                  >  1:357785/1‑251 (MQ=255)                     gCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTTTTTACGCTGCCGGGCTTTGATCTGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGc                                                                                                                                                                                                                                <  2:79938/251‑1 (MQ=255)                        gTTCATTGTAGCGCCGGGAATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCCAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGCTGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGCAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTACTGGAATTACTGCTTc                                                                                                                                                                                                                             >  2:282089/1‑251 (MQ=255)                                aaGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTAATGATTCAGGAagc                                                                                                                                                                                                                      >  2:317975/1‑250 (MQ=255)                                  gcgcCGGGTATTCGTGGAATGGTGATTTTGGTGTTTACGCTGCCGGTCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCTCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGTAATTACTGCTTCAGGAagcagc                                                                                                                                                                                                                   <  1:396862/251‑1 (MQ=255)                                           agtcttGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGTCAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACGGCTTCAGGACGCAGCGTaaaaaa                                                                                                                                                                                                           <  2:91425/245‑1 (MQ=255)                                                tGGAATGGTGATGTTGGTGTTTACGCTGCAGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCg                                                                                                                                                                                                     >  2:565549/1‑251 (MQ=255)                                                 ggAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGGCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGGGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTGTGTGCTGGAGAAGCGGCATATTTTCGCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGa                                                                                                                                                                                                    <  1:537641/251‑1 (MQ=255)                                                      ggTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGAtctc                                                                                                                                                                                                <  1:268364/250‑1 (MQ=255)                                                              gggTGTTTACGCTGCCGGGCTTTGATCGGGGATTCAAAGTCATCAAAGACAGGGTCGCGCCGCAGAAAGAGATGGCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACaa                                                                                                                                                                                       <  2:6107/250‑1 (MQ=255)                                                                        gcagCCGGGCTGTGATCGGGTAGTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTc                                                                                                                                                                             <  2:98367/248‑1 (MQ=255)                                                                        gCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGACCAAATTGTGATTc                                                                                                                                                                             >  1:446620/1‑251 (MQ=255)                                                                           gCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCg                                                                                                                                                                            <  1:362553/249‑1 (MQ=255)                                                                                    tGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTtat                                                                                                                                                                  >  1:364380/1‑250 (MQ=255)                                                                                         gggTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAg                                                                                                                                                            <  1:362485/251‑1 (MQ=255)                                                                                              ttCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAgcggcg                                                                                                                                                       >  1:8911/1‑251 (MQ=255)                                                                                                  aaGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCAGCTTTATATTGAGCGGCGGATg                                                                                                                                                   >  1:34534/1‑251 (MQ=255)                                                                                                         caAAGACAGGTTCGCGCCGCAGAAAGAGGTGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTTGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGc                                                                                                                                            <  2:463713/251‑1 (MQ=255)                                                                                                         caAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGc                                                                                                                                            >  1:171096/1‑251 (MQ=255)                                                                                                                                 agagATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTTGAACAAGTgg                                                                                                                    <  1:446756/251‑1 (MQ=255)                                                                                                                                             cGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGt                                                                                                        >  1:475710/1‑251 (MQ=255)                                                                                                                                             cGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGt                                                                                                        <  1:370922/251‑1 (MQ=255)                                                                                                                                                       cGTGCCTGCTATCAACTGGTGAAATAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATTGAATTACTGCGTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACg                                                                                               <  1:484199/250‑1 (MQ=255)                                                                                                                                                       cGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGGCAGCAGTTGCGCGACGc                                                                                              <  1:117084/251‑1 (MQ=255)                                                                                                                                                                     aCTGGGGAAAGAACTCGAGCGCGTGGGCCGAAGGGCGGACACCCAGGATTTTTAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCAGTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCTATCTTTATATTGAGCGGCGGATGGTGCCGCTCAAAAGCTTGCGGGAACAAGTGGAAGGTCAGCAGTTGCGCGACCCCAGCGAATAATACg                                                                                <  2:245334/251‑1 (MQ=255)                                                                                                                                                                        ggTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACGCCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTa                                                                             <  2:42877/251‑1 (MQ=255)                                                                                                                                                                                           gTGGGCCGACTGGCAGCAACCCAGGAGTTTGAACACTTTTTGCAGGAGAAGCGGCATATTTACCAGGCATTACTGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATAGTGATTCGCCATATTTCTATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGACCAAGAGGAAGGTCAGAAGTTACGCGACGCCATTGAAGAATCAGGTAACGCTCTTCGCCAGCTTGc                                                          >  2:546697/1‑251 (MQ=255)                                                                                                                                                                                                 cGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGCAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCCATATCTGGCTGGAACTAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCgctgc                                                    >  2:184217/1‑251 (MQ=255)                                                                                                                                                                                                         ggACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGAGTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCGGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCAGTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACAttt                                            <  1:333669/251‑1 (MQ=255)                                                                                                                                                                                                                       ttGAAAACTTTGTGCTTGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACAt                               <  1:317975/250‑1 (MQ=255)                                                                                                                                                                                                                                    gCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGCATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACt                 >  1:167213/1‑251 (MQ=255)                                                                                                                                                                                                                                    gCTGGAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACt                 <  2:311089/251‑1 (MQ=255)                                                                                                                                                                                                                                       ggAGAAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCg              <  1:271889/251‑1 (MQ=255)                                                                                                                                                                                                                                        gagaAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCgg             >  1:300111/1‑251 (MQ=255)                                                                                                                                                                                                                                          gaAGCGGCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATCTTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGCCGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGCCATGCTGTTTAAAAACTTCGgtg           >  1:354369/1‑251 (MQ=255)                                                                                                                                                                                                                                               ggCATATTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATTCGGTAACGCTATTCGCCAGCTTGACGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCAcc      >  1:69188/1‑251 (MQ=255)                                                                                                                                                                                                                                                   tatTTTCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTc  >  1:14985/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      TATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTC  >  NC_000913/4219505‑4219997

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