Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,577,937	T→C	E147E (GAA→GAG)	maeB ←	malate dehydrogenase (oxaloacetate‑decarboxylating) (NADP(+))

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,577,937	0	T	C	97.6%	133.0 / NA	41	E147E (GAA→GAG)	maeB	malate dehydrogenase (oxaloacetate‑decarboxylating) (NADP(+))
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (17/23);  minor base A (1/0);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CACGGCCCTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAG  >  NC_000913/2577690‑2578182                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      cACGGCCCTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCa                                                                                                                                                                                                                                                    <  1:340956/251‑1 (MQ=255)        cTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGTATATTCATCCGCTCGCGCAGTTTCTGCTCAatat                                                                                                                                                                                                                                                >  1:627883/1‑248 (MQ=255)                      cTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCt                                                                                                                                                                                                                               >  1:509037/1‑251 (MQ=255)                           gAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTaa                                                                                                                                                                                                                           <  1:422174/250‑1 (MQ=255)                                 cAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTc                                                                                                                                                                                                                    <  1:593710/251‑1 (MQ=255)                                  aaaCCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGAGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCg                                                                                                                                                                                                                   >  2:447425/1‑251 (MQ=255)                                      cACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGt                                                                                                                                                                                                               >  1:420334/1‑251 (MQ=255)                                                            cccAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTc                                                                                                                                                                                         <  1:232465/251‑1 (MQ=255)                                                                    tACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCCTCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCgcgg                                                                                                                                                                                 >  1:645503/1‑251 (MQ=255)                                                                         gcagGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACa                                                                                                                                                                            <  1:561753/251‑1 (MQ=255)                                                                             gTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGTTGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACtt                                                                                                                                                                        >  2:407444/1‑251 (MQ=255)                                                                             gTTCATACAAGCGATTACTGCGGAACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGAGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATCCCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGGTGCCGCCGAAGGTTGGTTCGAGCGCGGCGaagactt                                                                                                                                                                        >  2:422610/1‑245 (MQ=255)                                                                                   ccaggctctggcggcagctcccGCTCCGGAAACCGCCATCCGCACGCCGGAGATGTTTTTCGCCACCACGCGCAAGCCGTTGAGGATGTCGGCAGTGATGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCTCGCAGTTTCTGCTCCATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTctatca                                                                                                                                                                  <  2:308003/232‑6 (MQ=255)                                                                                          aTTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTc                                                                                                                                                           >  1:679055/1‑251 (MQ=255)                                                                                                   gCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAg                                                                                                                                                  <  1:247188/251‑1 (MQ=255)                                                                                                      cccGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTc                                                                                                                                               <  1:642954/251‑1 (MQ=255)                                                                                                              ggAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACtt                                                                                                                                       >  2:484752/1‑251 (MQ=255)                                                                                                               gAAACCCCCATCCGAACTTCGGAGATGTTTTTCGCCCCCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGGCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAGCATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCTCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTc                                                                                                                                      <  2:393368/251‑1 (MQ=255)                                                                                                                                  cGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAAt                                                                                                                    <  1:281181/250‑1 (MQ=255)                                                                                                                                   ggAGTTGTTTTTCTCCACCACGCGCAAGCCGGTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGAGCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCATATACATCAATc                                                                                                                   <  2:548689/250‑1 (MQ=255)                                                                                                                                                  caccacGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTaa                                                                                                   <  1:407444/251‑1 (MQ=255)                                                                                                                                                  caccacGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCATCGTGGAATACCGGACTATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTAAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTa                                                                                                    >  1:649517/1‑250 (MQ=255)                                                                                                                                                    ccacGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAAt                                                                                                                    <  1:489764/232‑1 (MQ=255)                                                                                                                                                    ccacGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAAt                                                                                                                    >  2:489764/1‑232 (MQ=255)                                                                                                                                                          gcAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAAGACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTTATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACg                                                                                           <  1:654035/251‑1 (MQ=255)                                                                                                                                                                      ggATGGCGGCAGTGCTTATAAATGCCGTGCCGTGCGTATCGACGTGGAATCCCTGAATATTCTTCCGCTCGCGCAGTTTCTGCTCAATATAGAATCATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATTCATCAATCCCGTTGAATTTCTTAAACAGAACGCCCTTGCCTTcc                                                                               <  2:420334/251‑1 (MQ=255)                                                                                                                                                                             ggcAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCg                                                                        <  1:309238/251‑1 (MQ=255)                                                                                                                                                                              gcAGTTCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGTTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTAGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCTACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCgg                                                                       <  1:148782/251‑1 (MQ=255)                                                                                                                                                                              gcAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCGCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAATGTTGGTGCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCgg                                                                       <  1:447425/251‑1 (MQ=255)                                                                                                                                                                                       aTAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGTCGCTTTAATGTCTTCGAGGTTGATGCCGCCGaa                                                                                                                                                                                                  <  1:27599/119‑1 (MQ=255)                                                                                                                                                                                       aTAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGaa                                                                                                                                                                                                  >  2:27599/1‑119 (MQ=255)                                                                                                                                                                                         aaTTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGAAAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTttcctgccc                                                            >  2:474324/1‑250 (MQ=255)                                                                                                                                                                                                       gcTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGAGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCACGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTg                                              >  1:524047/1‑251 (MQ=255)                                                                                                                                                                                                          gATCTTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTGAATGTCTTCGATGTTGATGCCGCCTAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACGGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGcc                                            <  1:648727/250‑1 (MQ=255)                                                                                                                                                                                                          gATCGTCGAGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGATCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTAGAGCGCGGCGACAACTTCAATAACTTTGTCCGGGTAGAGTTCGTCAACTTCAATGTCAAATACCTCAATCCCGGCGAATTTCTTAAAAAGAACGCACTTGCCTTCAATCACAGGTTTGCCTGCCAGCGCGCCATTGTTGcc                                            >  2:505055/1‑250 (MQ=255)                                                                                                                                                                                                                cggtgAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCCTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTATcccc                                     <  2:509037/247‑1 (MQ=255)                                                                                                                                                                                                                      aTACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAcc                               <  1:484752/251‑1 (MQ=255)                                                                                                                                                                                                                      aTACCGGAATATTCATCCGCTCGCGCAGATTCTGCTCAATATAGAAACATGCTGGCGCTTTAATGTCTTCTAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTTTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAAGCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGccct                               <  2:443394/251‑4 (MQ=255)                                                                                                                                                                                                                                              gcAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTCCCGTTAGAGATCACCGcc       >  1:261142/1‑251 (MQ=255)                                                                                                                                                                                                                                               cAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCaa                                                   >  1:174118/1‑206 (MQ=255)                                                                                                                                                                                                                                                   ttCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGCTGCCGCCGAAGGTTAGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGTGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGAGACAATGTTGCCTAACCCCAGCACCGCCTTACCGTTTGAGATCACCGccacccg  >  2:111822/1‑249 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      CACGGCCCTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAG  >  NC_000913/2577690‑2578182

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