Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,756,553	A→T	100%	K286* (AAA→TAA)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,553	0	A	T	98.9%	304.1 / ‑3.9	88	K286* (AAA→TAA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (41/46);  total (41/47)
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.27e-01

AGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAG  >  NC_000913/1756456‑1756649                                                                                                  |                                                                                                 agaagaAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCCGATGCTGGATTCCATGATCtaa                                                                                                >  2:1844323/1‑100 (MQ=255) agaagaAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaa                                                                                                >  1:1200236/1‑100 (MQ=255)  gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaa                                                                                               <  2:1880241/100‑1 (MQ=255)   aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaaa                                                                                              <  1:2028835/100‑1 (MQ=255)   aagaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaaa                                                                                              <  2:3503231/100‑1 (MQ=255)      aagTTATCTTCTCCCAGAAAAAAAAAAAAAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAccc                                                                                           <  1:3350607/100‑1 (MQ=255)       agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCa                                                                                          >  2:866724/1‑100 (MQ=255)         ttATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACg                                                                                        <  1:983248/100‑1 (MQ=255)           aTCTTCGCCCAGAATATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGcc                                                                                      <  1:1631837/100‑1 (MQ=255)           aTCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGcc                                                                                      <  2:2728347/100‑1 (MQ=255)             cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCg                                                                                    >  2:2248393/1‑100 (MQ=255)             cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCg                                                                                    <  2:1638189/100‑1 (MQ=255)                cGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACt                                                                                 >  1:1167607/1‑100 (MQ=255)                   ccAGAAGATGATGATCGAAAACTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGGTGCTGGATTCCATGATCTAAAACCCACGCCCGAc                                                                                  <  1:62045/96‑1 (MQ=255)                   ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgc                                                                              >  1:1494685/1‑100 (MQ=255)                    cAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgcg                                                                             <  2:1094214/100‑1 (MQ=255)                     agaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCgcgc                                                                            >  2:509109/1‑100 (MQ=255)                      gaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCa                                                                           <  1:1715739/100‑1 (MQ=255)                          atgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAg                                                                       <  2:449481/100‑1 (MQ=255)                           tgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGc                                                                      >  2:2502632/1‑100 (MQ=255)                            gatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGcc                                                                     <  1:755976/100‑1 (MQ=255)                             atgatCGAAAAATGTATCCGTGCACTTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                                    <  2:892792/100‑1 (MQ=255)                              tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                                   <  2:1176909/100‑1 (MQ=255)                              tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                                   <  2:1747623/100‑1 (MQ=255)                              tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                                   >  2:1011629/1‑100 (MQ=255)                                atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                                 <  1:163756/100‑1 (MQ=255)                                 tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                                                >  2:1961502/1‑100 (MQ=255)                                 tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCAATGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                                                >  2:2633314/1‑100 (MQ=255)                                  cGAAAAATGTGTCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAACACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                                               <  1:3365791/100‑1 (MQ=255)                                  cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                                               <  2:790619/100‑1 (MQ=255)                                   gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGAGACCCAGATGCTGGATTCCATGATCTAAAACGCACGCCCGACTCGCGCAGAAGCCGGTGACg                                                              >  1:3042348/1‑100 (MQ=255)                                    aaaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACgn                                                             <  2:2909015/100‑2 (MQ=255)                                       aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACATTGc                                                          <  1:294479/100‑1 (MQ=255)                                        aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAGCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCa                                                         <  2:252466/100‑1 (MQ=255)                                           tATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAAc                                                      <  1:1011624/100‑1 (MQ=255)                                            aTCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACg                                                     <  1:966884/100‑1 (MQ=255)                                            aTCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACg                                                     >  1:3540001/1‑100 (MQ=255)                                              ccGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGcc                                                   <  2:2388936/100‑1 (MQ=255)                                                  gCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                                               >  1:2364891/1‑100 (MQ=255)                                                  gCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                                               <  2:1403656/100‑1 (MQ=255)                                                  gCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATcc                                               <  2:3571348/100‑1 (MQ=255)                                                     cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                            <  2:88017/100‑1 (MQ=255)                                                     cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                            >  1:3121097/1‑100 (MQ=255)                                                       tAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                                            <  1:2935545/98‑1 (MQ=255)                                                       tAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                          >  2:477877/1‑100 (MQ=255)                                                       tAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                          >  1:2523230/1‑100 (MQ=255)                                                       tAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                          <  1:3545734/100‑1 (MQ=255)                                                           gTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACgg                                        >  1:2261177/1‑98 (MQ=255)                                                           gTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTa                                      >  2:1188153/1‑100 (MQ=255)                                                            tCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTAc                                     <  2:1167613/100‑1 (MQ=255)                                                                tATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc                                 >  1:880137/1‑100 (MQ=255)                                                                tATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc                                 >  2:36894/1‑100 (MQ=255)                                                                tATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGAGACGTTGCAAACGCCATCCTCGACGGTACTGAc                                 >  1:3257354/1‑100 (MQ=255)                                                                 aTCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACg                                <  1:17220/100‑1 (MQ=255)                                                                  tCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                               <  1:2502618/100‑1 (MQ=255)                                                                   cACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCa                              >  1:977724/1‑100 (MQ=255)                                                                    aCTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                             >  2:3363222/1‑100 (MQ=255)                                                                    aCTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCAGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGccg                             >  1:1420626/1‑98 (MQ=255)                                                                     cTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt                            >  2:784740/1‑100 (MQ=255)                                                                     cTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt                            >  1:75139/1‑100 (MQ=255)                                                                     cTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt                            >  2:2069413/1‑100 (MQ=255)                                                                      tGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTg                           >  2:1998119/1‑100 (MQ=255)                                                                          aCCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc                       <  1:207857/100‑1 (MQ=255)                                                                          aCCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc                       <  2:1471504/100‑1 (MQ=255)                                                                           nccAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCt                      >  1:3530506/2‑100 (MQ=255)                                                                           cccAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGca                      >  1:3227512/1‑99 (MQ=255)                                                                           cccAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc                       <  2:3227534/99‑1 (MQ=255)                                                                            ccAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTg                     >  2:3439120/1‑100 (MQ=255)                                                                              aGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTc                   <  1:36893/100‑1 (MQ=255)                                                                              aGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTc                   >  1:1313385/1‑100 (MQ=255)                                                                              aGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTc                   <  2:1567520/100‑1 (MQ=255)                                                                              aGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTc                   >  1:1530323/1‑100 (MQ=255)                                                                               gctgttGGATTTCATGATCTAAAACCTCCGTCCGAGTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGATGTc                   <  1:1601665/94‑1 (MQ=255)                                                                                 tGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg                >  2:3715551/1‑100 (MQ=255)                                                                                 tGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg                <  2:1530330/100‑1 (MQ=255)                                                                                  gttGGACTCCATGATCTAAAACCCACGCCCGACTAGGGCAGAAGCCGGTGATGTTGTAAACGCCATCCTCGACGGTAATGACGAAGTGATGCTGTCTg                 <  1:2070961/96‑1 (MQ=255)                                                                                  gCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg                <  2:467041/99‑1 (MQ=255)                                                                                  gCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGt               <  1:408972/100‑1 (MQ=255)                                                                                  gCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGt               <  1:3255501/100‑1 (MQ=255)                                                                                   cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAATCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTAATGCTGTCTg                 >  1:986332/1‑97 (MQ=255)                                                                                   cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTggcg              <  2:2523244/100‑3 (MQ=255)                                                                                     ggATTCCATGATCTAAAACCCAGGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa            >  1:2166175/1‑100 (MQ=255)                                                                                     ggATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa            <  2:3085926/100‑1 (MQ=255)                                                                                     ggATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa            <  1:868268/100‑1 (MQ=255)                                                                                      gATTCCATGATCTAAAACCCACGCCCGCCTCGCGCAGAAGCCGGTGACGTTGCCAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAAt           >  2:1054078/1‑100 (MQ=255)                                                                                       aTTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATc          <  2:2607262/100‑1 (MQ=255)                                                                                          ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGc       >  1:1069819/1‑100 (MQ=255)                                                                                          ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGc       >  2:2650042/1‑100 (MQ=255)                                                                                          ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACCCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAAAcc         >  1:3232238/1‑98 (MQ=255)                                                                                           cATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa      >  2:2089755/1‑100 (MQ=255)                                                                                            aTGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaa     <  2:3484464/100‑1 (MQ=255)                                                                                            aTGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaa     >  1:606208/1‑100 (MQ=255)                                                                                               aTCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaaan  <  2:3530532/100‑2 (MQ=255)                                                                                                  |                                                                                                 AGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAAG  >  NC_000913/1756456‑1756649

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