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	95.3%	277.0 / 3.9	85	K286* (AAA→TAA)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (0/4);  new base T (40/41);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

GAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAA  >  NC_000913/1756457‑1756648                                                                                                 |                                                                                                gaagaaGTTATCTTCGCCCAGAAGATGATTATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaa                                                                                              <  2:426558/100‑1 (MQ=255) gaagaaGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCtaaa                                                                                              >  2:2743121/1‑100 (MQ=255)   agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAc                                                                                            <  2:2321357/100‑1 (MQ=255)   agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAc                                                                                            <  1:1783536/100‑1 (MQ=255)     aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAAccc                                                                                          <  1:312639/100‑1 (MQ=255)      agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCa                                                                                         <  2:2120150/100‑1 (MQ=255)       gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCAc                                                                                        >  1:106505/1‑100 (MQ=255)         tATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACACAGATGCTGGATTCCATGATCTAAAACCCAc                                                                                        >  2:327074/1‑98 (MQ=255)             ttCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGa                                                                                  >  2:2726582/1‑100 (MQ=255)             ttCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGa                                                                                  >  1:2331320/1‑100 (MQ=255)              tCGCCTAGAAGATGATGATCGAAGAATGTACCCGTGCACGTGAAGGCGTTATCACTGCCACCCAGATGCTGGATTCCATGGTCTAAAGCCCACGCCCGAc                                                                                 <  2:2380940/100‑1 (MQ=255)                gCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTc                                                                               <  1:836594/100‑1 (MQ=255)                gCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTc                                                                               >  1:708497/1‑100 (MQ=255)                   cAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgcg                                                                            <  1:1073240/100‑1 (MQ=255)                   cAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTcgc                                                                             >  2:920763/1‑99 (MQ=255)                     gaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCa                                                                          >  1:404636/1‑100 (MQ=255)                     gaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAAAAAAAAAAACTGAGAAAAAAATGATGGATTCCATGATCTAAAACCCACGCCCGACTcgcgn                                                                           <  2:1011268/99‑2 (MQ=255)                      aagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAg                                                                         <  2:2654289/100‑1 (MQ=255)                          tgatgatCGAAAAATGTATCCGTGCACTTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGaa                                                                       >  2:2650004/1‑98 (MQ=255)                          tgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGc                                                                     >  1:687010/1‑100 (MQ=255)                          tgatgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGc                                                                     >  2:817966/1‑100 (MQ=255)                             tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                                  <  1:2248728/100‑1 (MQ=255)                             tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                                  >  2:700402/1‑100 (MQ=255)                             tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCAATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                                  <  2:1700961/100‑1 (MQ=255)                              gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                                                 >  1:195390/1‑100 (MQ=255)                              gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                                                 >  1:1597020/1‑100 (MQ=255)                              gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                                                 <  2:1581334/100‑1 (MQ=255)                               atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                                <  1:2570571/100‑1 (MQ=255)                               aaagAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                                                <  2:391915/97‑1 (MQ=255)                                tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                                               >  1:823920/1‑100 (MQ=255)                                 cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                                              >  1:2505465/1‑100 (MQ=255)                                  gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg                                                             >  1:2586663/1‑100 (MQ=255)                                  gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg                                                             >  2:1394181/1‑100 (MQ=255)                                  gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg                                                             <  2:2675785/100‑1 (MQ=255)                                  gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg                                                             >  1:949348/1‑100 (MQ=255)                                    aaaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACAtt                                                           <  1:526705/100‑1 (MQ=255)                                     aaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                                          <  2:358260/100‑1 (MQ=255)                                     aaaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTg                                                          <  2:1979232/100‑1 (MQ=255)                                      aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                                         <  2:2635398/100‑1 (MQ=255)                                      aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                                         <  1:322115/100‑1 (MQ=255)                                        tGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCaa                                                       <  1:2022613/100‑1 (MQ=255)                                          tATCCGTGCACGTAAAGTCGTTATCACTGCTACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGTTGACGTTGCAAAc                                                     <  2:2578271/100‑1 (MQ=255)                                          tATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAAc                                                     >  2:85995/1‑100 (MQ=255)                                            tCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGc                                                   >  1:915703/1‑100 (MQ=255)                                              cGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCa                                                 <  2:1090290/100‑1 (MQ=255)                                              cGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCa                                                 <  2:1092184/100‑1 (MQ=255)                                               gTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCAt                                                <  2:1268167/100‑1 (MQ=255)                                                   aCGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTc                                            >  1:2672470/1‑100 (MQ=255)                                                    cGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTACATGAGCTAAAACCCACGCCCGACTCGCGAAGAAGCCAGAGACGTTGCAAACGcc                                                  >  2:2484473/1‑93 (MQ=255)                                                      tAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                         <  1:1381034/100‑1 (MQ=255)                                                       aaaGTCGTTATCACTGCACCCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc                                         <  1:1221901/99‑1 (MQ=255)                                                       aaaGTCGTTACCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACg                                        <  1:1409830/100‑1 (MQ=255)                                                            cGTTATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt                                   >  1:552193/1‑100 (MQ=255)                                                               tATCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc                                >  2:467670/1‑100 (MQ=255)                                                                aTCACTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCAGGTGAAGTTGCAAACGCCATCCTCGACGGTACTGACg                               >  2:212941/1‑100 (MQ=255)                                                                   aCTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                            >  1:908175/1‑100 (MQ=255)                                                                   aCTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg                            >  1:889660/1‑100 (MQ=255)                                                                    cTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt                           >  1:2515325/1‑100 (MQ=255)                                                                    cTGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGt                           <  2:2515307/100‑1 (MQ=255)                                                                     tGCGACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTg                          <  1:1431517/100‑1 (MQ=255)                                                                       cGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGAt                        <  1:475720/100‑1 (MQ=255)                                                                        gACCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATg                       <  2:708495/100‑1 (MQ=255)                                                                         aCCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc                      >  2:1530526/1‑100 (MQ=255)                                                                         aCCCAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGc                      >  2:78379/1‑100 (MQ=255)                                                                          cccAGATGCTGGATTCCATGATCTAAAACCCACGCCCGTCTCGCGCAGCAGCCGGTGACGTTCCAAACGCCATCCTCGCCGGTCCTGCCGCAGTGATGc                      >  2:1692119/1‑99 (MQ=255)                                                                            cAGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGAtacggt                   >  2:1358459/1‑95 (MQ=255)                                                                             aGATGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTc                  <  1:1744173/100‑1 (MQ=255)                                                                               aTGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg                <  1:650061/100‑1 (MQ=255)                                                                               aTGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg                >  2:516767/1‑100 (MQ=255)                                                                                tGCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg               >  1:2423219/1‑100 (MQ=255)                                                                                 gCTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGt              <  2:1166410/100‑1 (MQ=255)                                                                                  ntGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg             >  2:378386/2‑100 (MQ=255)                                                                                  cTGGATTCCATGGTCTAAAACCCGCGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg             <  2:1051894/100‑1 (MQ=255)                                                                                  cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg             <  2:2296403/100‑1 (MQ=255)                                                                                  cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg             <  1:2153813/100‑1 (MQ=255)                                                                                  cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg             <  1:404963/100‑1 (MQ=255)                                                                                  cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg             >  1:1330301/1‑100 (MQ=255)                                                                                  cTGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg             <  1:1456444/100‑1 (MQ=255)                                                                                  cTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg             <  2:2702160/100‑1 (MQ=255)                                                                                   tGGATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg             <  2:1817177/99‑1 (MQ=255)                                                                                    ggATTCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa           <  2:195390/100‑1 (MQ=255)                                                                                       ttCCATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATcc        <  2:1798195/100‑1 (MQ=255)                                                                                         ncATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGc      >  2:685546/2‑100 (MQ=255)                                                                                         ccATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGCCGCAGTGATGCTGTCTGGTGAATCCGc      >  1:2192840/1‑100 (MQ=255)                                                                                          cATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa     <  2:2274536/100‑1 (MQ=255)                                                                                          cATGATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa     >  1:309540/1‑100 (MQ=255)                                                                                             gATCTAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCaaaa  >  1:916912/1‑100 (MQ=255)                                                                                                 |                                                                                                GAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCAAAA  >  NC_000913/1756457‑1756648

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