Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,756,495 (A)5→6 coding (798/1413 nt) pykF → pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,756,4901.A100.0% 227.1 / NA 58E265K (GAA→AAA) pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (26/32);  total (26/32)

TCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGG  >  NC_000913/1756341‑1756635
                                                                                                                                                      |                                                                                                                                                 
tCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg                                                                                                                                                    >  2:844609/1‑150 (MQ=255)
tCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCg                                                                                                                                                    >  2:1636/1‑150 (MQ=255)
       catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaaaatg                                                                                                                                           <  1:1027347/151‑4 (MQ=255)
       catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaaaat                                                                                                                                            <  1:210558/150‑3 (MQ=255)
       catcTCCAAAATCGAAAACCAAGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaaaat                                                                                                                                            <  1:801892/150‑3 (MQ=255)
            ccAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGaaaaaatgta                                                                                                                                         <  2:502226/148‑6 (MQ=255)
                  tCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCAc                                                                                                                                 >  1:654092/1‑151 (MQ=255)
                  tCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCAc                                                                                                                                 >  1:801013/1‑151 (MQ=255)
                   cGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACg                                                                                                                                <  2:1051790/151‑1 (MQ=255)
                    gAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGt                                                                                                                               >  1:344135/1‑151 (MQ=255)
                      aaaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTaa                                                                                                                             <  1:479258/151‑1 (MQ=255)
                      aaaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTa                                                                                                                              >  2:567531/1‑150 (MQ=255)
                       aaCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTaaa                                                                                                                            <  1:143234/151‑1 (MQ=255)
                        aCCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTa                                                                                                                              >  1:783735/1‑148 (MQ=255)
                         ccAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAg                                                                                                                           <  2:813578/150‑1 (MQ=255)
                           aGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCg                                                                                                                        >  1:267577/1‑151 (MQ=255)
                            ggAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGt                                                                                                                       <  2:682975/151‑1 (MQ=255)
                             gAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGtt                                                                                                                      <  2:929656/151‑1 (MQ=255)
                              aaGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGtt                                                                                                                      >  2:211274/1‑150 (MQ=255)
                                ggCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATc                                                                                                                   >  1:769581/1‑151 (MQ=255)
                                 gCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATc                                                                                                                   <  2:1083926/150‑1 (MQ=255)
                                       acaacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGa                                                                                                            <  2:22161/151‑1 (MQ=255)
                                         aacTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGAcc                                                                                                          >  1:993487/1‑151 (MQ=255)
                                           cTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGAccc                                                                                                         >  2:138291/1‑150 (MQ=255)
                                                 cgaAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCt                                                                                                  <  1:25731/151‑1 (MQ=255)
                                                   aaaTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTgg                                                                                                >  1:292183/1‑151 (MQ=255)
                                                   aaaTCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTgg                                                                                                <  1:916072/151‑1 (MQ=255)
                                                        cTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTcc                                                                                           >  1:1081694/1‑151 (MQ=255)
                                                            aaGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATg                                                                                        >  2:574261/1‑150 (MQ=255)
                                                                     aCGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAAcc                                                                              <  1:96053/151‑1 (MQ=255)
                                                                         catcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACg                                                                          <  1:976837/151‑1 (MQ=255)
                                                                         catcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACg                                                                          <  1:1067076/151‑1 (MQ=255)
                                                                         catcatGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACg                                                                          <  1:270180/151‑1 (MQ=255)
                                                                            catGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGccc                                                                       <  2:245922/151‑1 (MQ=255)
                                                                              tGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCg                                                                      >  2:865899/1‑150 (MQ=255)
                                                                                gTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACt                                                                   <  1:579333/151‑1 (MQ=255)
                                                                                gTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACt                                                                   <  1:131765/151‑1 (MQ=255)
                                                                                   gcgcGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTcgc                                                                <  1:303141/151‑1 (MQ=255)
                                                                                            gACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGcc                                                       <  2:930253/151‑1 (MQ=255)
                                                                                            gACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGcc                                                       <  1:1110956/151‑1 (MQ=255)
                                                                                             aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                      >  2:109553/1‑151 (MQ=255)
                                                                                             aCCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCg                                                      <  1:674338/151‑1 (MQ=255)
                                                                                              ccTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCgg                                                     <  1:416763/151‑1 (MQ=255)
                                                                                                    tgtAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGt                                               >  2:485759/1‑151 (MQ=255)
                                                                                                       aGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                            >  1:785892/1‑151 (MQ=255)
                                                                                                        gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                            >  1:150672/1‑150 (MQ=255)
                                                                                                        gAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGc                                            >  1:572282/1‑150 (MQ=255)
                                                                                                         aaaTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCaa                                          <  1:207766/151‑1 (MQ=255)
                                                                                                           aTCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAAc                                        >  1:1010516/1‑151 (MQ=255)
                                                                                                             cccGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACg                                       <  1:467681/150‑1 (MQ=255)
                                                                                                              ccGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGcc                                     <  1:983048/151‑1 (MQ=255)
                                                                                                                          gTTATCTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGt                         >  1:786764/1‑151 (MQ=255)
                                                                                                                               cTTCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTg                     >  1:954866/1‑150 (MQ=255)
                                                                                                                                 tCGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACt                      <  2:375398/147‑1 (MQ=255)
                                                                                                                                  cGCCCAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGc                 <  2:647583/151‑1 (MQ=255)
                                                                                                                                     ccAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAg               >  2:645249/1‑150 (MQ=255)
                                                                                                                                      cAGAAGATGATGATCGAAAAAATGTATCCGTGCGCGTAAAGTCGTTCTCACTGCGACGCAGATGCTGGTTTCCATGATCAAAAACCCCCGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTg             <  2:964701/151‑1 (MQ=255)
                                                                                                                                      cAGAAGATGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTg             <  2:299294/151‑1 (MQ=255)
                                                                                                                                        gaagaTGATGATCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGAt           <  1:913582/151‑1 (MQ=255)
                                                                                                                                          agaTGATGATCGAAAAAATGTATCCGTGCCCGTAAATTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGa            <  2:23319/148‑1 (MQ=255)
                                                                                                                                            atgatgatCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCt        >  2:575191/1‑150 (MQ=255)
                                                                                                                                               atgatCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCt    >  2:979049/1‑151 (MQ=255)
                                                                                                                                                tgatCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg   >  2:582180/1‑151 (MQ=255)
                                                                                                                                                 gatCGAAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg  <  1:142641/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
TCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAAGTCGTTATCACTGCGACCCAGATGCTGGATTCCATGATCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGG  >  NC_000913/1756341‑1756635

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