Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,529	A→T	T278S (ACC→TCC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,529	0	A	T	95.6%	301.8 / 4.9	90	T278S (ACC→TCC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base T (39/47);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

GACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCA  >  NC_000913/1756433‑1756621                                                                                                 |                                                                                             gACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTccc                                                                                           >  2:1193248/1‑100 (MQ=255)  aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCa                                                                                          <  2:635458/100‑1 (MQ=255)  aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCa                                                                                          <  2:1737339/100‑1 (MQ=255)  aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCa                                                                                          >  2:1379635/1‑100 (MQ=255)  aCCTGGATGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCa                                                                                          >  1:1256182/1‑100 (MQ=255)     tGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGAt                                                                                       >  2:408037/1‑100 (MQ=255)         tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTg                                                                                   <  1:1897076/100‑1 (MQ=255)          gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTgg                                                                                  <  2:2525200/100‑1 (MQ=255)          gtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTgg                                                                                  <  2:1967588/100‑1 (MQ=255)           tAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGa                                                                                 <  2:990240/100‑1 (MQ=255)           tAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGa                                                                                 >  1:924212/1‑100 (MQ=255)            aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAt                                                                                <  2:2175658/100‑1 (MQ=255)            aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGAt                                                                                >  1:300865/1‑100 (MQ=255)               aaTCCCGGTAGAAGAAGTTATCTTCGCCCAGATGTTTTTGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTcc                                                                             <  2:160090/100‑1 (MQ=255)               aaTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTcc                                                                             >  1:398804/1‑100 (MQ=255)               aaTCCCCGTAGAAGAAGTTAGCTTCGCACAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTcc                                                                             <  1:1585451/100‑1 (MQ=255)                 tCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCAt                                                                           <  2:493842/100‑1 (MQ=255)                    cGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGAt                                                                        <  2:401291/100‑1 (MQ=255)                     ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATc                                                                       <  1:2059952/100‑1 (MQ=255)                     ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATc                                                                       <  1:2555546/100‑1 (MQ=255)                     ggTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATc                                                                       <  2:1389182/100‑1 (MQ=255)                      gTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCa                                                                      <  2:1737311/100‑1 (MQ=255)                      gTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCa                                                                      <  1:997516/100‑1 (MQ=255)                       tAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCaa                                                                     <  1:753198/100‑1 (MQ=255)                           agaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAc                                                                 <  1:609849/100‑1 (MQ=255)                            gaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                                                >  1:984668/1‑100 (MQ=255)                            gaagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                                                <  2:158541/100‑1 (MQ=255)                             aagTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAAccc                                                               <  2:1841876/100‑1 (MQ=255)                              agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCa                                                              <  1:438383/100‑1 (MQ=255)                              agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCa                                                              >  2:1287575/1‑100 (MQ=255)                              agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCa                                                              <  2:924215/100‑1 (MQ=255)                              agTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCa                                                              >  2:1934151/1‑100 (MQ=255)                               gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCAc                                                             <  1:1232143/100‑1 (MQ=255)                               gTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCAc                                                             >  1:1896067/1‑100 (MQ=255)                                   tCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGccc                                                         >  1:625163/1‑100 (MQ=255)                                   tCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGccc                                                         <  1:2360502/100‑1 (MQ=255)                                    cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                        <  2:52250/100‑1 (MQ=255)                                    cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                        <  1:596515/100‑1 (MQ=255)                                    cTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                        >  2:2083569/1‑100 (MQ=255)                                    cTTCGCCCAGAAGATGATGAACGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                        >  1:1593260/1‑100 (MQ=255)                                         cccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcg                                                   >  1:1289265/1‑100 (MQ=255)                                          ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcgc                                                  <  2:464142/100‑1 (MQ=255)                                          ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcgc                                                  <  1:2360506/100‑1 (MQ=255)                                          ccAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcgc                                                  <  2:1896077/100‑1 (MQ=255)                                            agaagaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCgcgc                                                >  2:962540/1‑100 (MQ=255)                                              cttccgatcttCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAg                                              <  1:1245707/90‑1 (MQ=255)                                               agaTGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCACAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGa                                             >  1:1278041/1‑100 (MQ=255)                                                    atgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                        >  2:1686816/1‑100 (MQ=255)                                                     tgatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATACCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCgn                                       <  2:2382390/100‑2 (MQ=255)                                                      gatCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGt                                      >  2:2246217/1‑100 (MQ=255)                                                       atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGatcggaaga                                     >  1:2422784/1‑92 (MQ=255)                                                       atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                     >  2:2069696/1‑100 (MQ=255)                                                       atCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTg                                     >  1:851483/1‑100 (MQ=255)                                                        tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGagatcggaag                                    >  2:1245711/1‑91 (MQ=255)                                                        tCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGa                                    >  1:2154936/1‑100 (MQ=255)                                                         cGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGAc                                   <  2:2259433/100‑1 (MQ=255)                                                          gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg                                  >  2:2010513/1‑100 (MQ=255)                                                          gAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACg                                  >  1:2254652/1‑100 (MQ=255)                                                              aaTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                              <  1:633629/100‑1 (MQ=255)                                                               aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCa                             >  2:24380/1‑100 (MQ=255)                                                               aTGTATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCa                             <  2:1985018/100‑1 (MQ=255)                                                                  tATCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAAc                          >  2:598483/1‑100 (MQ=255)                                                                   aTCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACg                         >  1:719319/1‑100 (MQ=255)                                                                    tCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGc                        >  1:1158230/1‑100 (MQ=255)                                                                    tCCGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGc                        >  1:2282658/1‑100 (MQ=255)                                                                     ccGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGcc                       <  2:172376/100‑1 (MQ=255)                                                                     ccGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGcc                       <  1:2069684/100‑1 (MQ=255)                                                                     ccGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGcc                       <  1:1379629/100‑1 (MQ=255)                                                                     ccGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGcc                       <  1:2543868/100‑1 (MQ=255)                                                                      cGTGCACGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCa                      >  2:566583/1‑100 (MQ=255)                                                                           aCGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTc                 <  1:24380/100‑1 (MQ=255)                                                                            cGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                >  2:1003642/1‑100 (MQ=255)                                                                            cGTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCg                <  1:254030/100‑1 (MQ=255)                                                                             ntAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGa               >  2:1362663/2‑100 (MQ=255)                                                                             gTAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGa               >  1:2228987/1‑100 (MQ=255)                                                                              tAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc              <  2:387980/100‑1 (MQ=255)                                                                              tAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc              >  1:334026/1‑100 (MQ=255)                                                                              tAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc              <  2:2254666/100‑1 (MQ=255)                                                                              tAAAGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATAATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc              <  1:566582/100‑1 (MQ=255)                                                                              cgatctCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAc              <  2:1391083/95‑1 (MQ=255)                                                                               aaaGTCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACg             >  2:1923780/1‑100 (MQ=255)                                                                                   tCGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTAc         <  1:1237593/100‑1 (MQ=255)                                                                                    cGTTATCACTGCGTCCCAGATGCTGGATTCCATGATCCAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGAcagatcg        >  1:1391077/1‑94 (MQ=255)                                                                                     gTTATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTg       <  1:962537/100‑1 (MQ=255)                                                                                     gTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTg       >  2:2101351/1‑100 (MQ=255)                                                                                       tATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc     <  2:2282675/100‑1 (MQ=255)                                                                                       tATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc     >  2:672012/1‑100 (MQ=255)                                                                                       tATCACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc     >  2:1720370/1‑100 (MQ=255)                                                                                          cACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCa  <  1:1193871/100‑1 (MQ=255)                                                                                          cACTGCGTCCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCa  <  1:1720360/100‑1 (MQ=255)                                                                                                 |                                                                                             GACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCA  >  NC_000913/1756433‑1756621

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