Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,756,894 +TG coding (1197/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,8921.T92.1% 272.4 / 15.4 76coding (1195/1413 nt)pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base T (31/39);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
*NC_0009131,756,8922.G92.1% 283.0 / 16.7 76coding (1195/1413 nt)pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base G (31/39);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

ACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑TGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAG  >  NC_000913/1756743‑1757033
                                                                                                                                                      ||                                                                                                                                             
aCCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑t                                                                                                                                              >  2:942484/1‑151 (MQ=255)
  cGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑Tgt                                                                                                                                            <  2:989150/151‑2 (MQ=255)
     aaaCTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑Tgtga                                                                                                                                          <  2:398953/150‑4 (MQ=255)
      aaCTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑Tgtgacc                                                                                                                                        >  2:205270/1‑146 (MQ=255)
       aCTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑Tgtgacc                                                                                                                                        >  1:151705/1‑145 (MQ=255)
           cgcATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑Tgtgaccacca                                                                                                                                    >  2:375397/1‑141 (MQ=255)
            gcATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAAc                                                                                                                                    >  1:778415/1‑151 (MQ=255)
              aTTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACg                                                                                                                                   >  1:867405/1‑150 (MQ=255)
                     aaGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACg                                                                                                                           <  2:598193/151‑1 (MQ=255)
                        cGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCt                                                                                                                        <  1:327503/151‑1 (MQ=255)
                        cGGTATGCCGCGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGGGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGc                                                                                                                         >  2:267791/1‑150 (MQ=255)
                          gTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCtca                                                                                                                      <  2:447998/151‑1 (MQ=255)
                          gTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCtc                                                                                                                       >  2:223271/1‑150 (MQ=255)
                              gCCGTGGTTCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCtcatc                                                                                                                    >  2:375058/1‑149 (MQ=255)
                               ccGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGt                                                                                                                 <  2:115939/151‑1 (MQ=255)
                               ccGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGt                                                                                                                 <  1:237727/151‑1 (MQ=255)
                                 gtggtgCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGtt                                                                                                                <  1:188553/150‑1 (MQ=255)
                                 gtggtgCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTg                                                                                                               >  1:556712/1‑151 (MQ=255)
                                 gtggtgCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTg                                                                                                               <  1:14569/151‑1 (MQ=255)
                                   ggtgCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGt                                                                                                             >  2:26140/1‑151 (MQ=255)
                                      gCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACt                                                                                                          <  2:408756/151‑1 (MQ=255)
                                         gTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAg                                                                                                       <  2:955069/151‑1 (MQ=255)
                                         gTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑TGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGc                                                                                                      >  2:272525/1‑150 (MQ=255)
                                            gAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCa                                                                                                     >  1:695920/1‑150 (MQ=255)
                                            gAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGGACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCaa                                                                                                    <  1:585864/151‑1 (MQ=255)
                                              aaCTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCaa                                                                                                    >  1:544213/1‑149 (MQ=255)
                                               aCTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAgg                                                                                                 >  2:785049/1‑151 (MQ=255)
                                                  gctgAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGt                                                                                              <  2:698249/151‑1 (MQ=255)
                                                     gAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTtgt                                                                                           <  2:497799/151‑1 (MQ=255)
                                                         aaCTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCg                                                                                       >  1:916339/1‑151 (MQ=255)
                                                            tGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGc                                                                                      >  1:767847/1‑149 (MQ=255)
                                                              gATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCt                                                                                  <  2:561104/151‑1 (MQ=255)
                                                                 gCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGt                                                                               <  2:512586/151‑1 (MQ=255)
                                                                 gCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGt                                                                               >  1:977578/1‑151 (MQ=255)
                                                                  cTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCCCCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTg                                                                                <  2:451060/149‑1 (MQ=255)
                                                                    ccGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGtt                                                                              <  2:304978/149‑1 (MQ=255)
                                                                    ccGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTa                                                                             <  1:752819/150‑1 (MQ=255)
                                                                      gCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAag                                                                          <  1:1036399/151‑1 (MQ=255)
                                                                      gCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAag                                                                          <  1:491612/151‑1 (MQ=255)
                                                                       cTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTaa                                                                            >  2:707646/1‑148 (MQ=255)
                                                                         gATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAgaga                                                                       <  2:463515/151‑1 (MQ=255)
                                                                           tCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAag                                                                          >  2:759727/1‑146 (MQ=255)
                                                                            cGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCa                                                                    <  1:201341/151‑1 (MQ=255)
                                                                              tGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATc                                                                     <  1:228529/148‑1 (MQ=255)
                                                                               ggTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACt                                                                  <  1:727962/150‑1 (MQ=255)
                                                                                   gCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTAc                                                             >  1:829421/1‑151 (MQ=255)
                                                                                   gCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTAc                                                             <  2:492669/151‑1 (MQ=255)
                                                                                      aCTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtg                                                           <  1:627280/150‑1 (MQ=255)
                                                                                      aCTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACt                                                            <  1:627778/149‑1 (MQ=255)
                                                                                        tCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtgatg                                                        >  2:61794/1‑151 (MQ=255)
                                                                                           gggCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGAtt                                                     >  2:103951/1‑151 (MQ=255)
                                                                                               ggTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCt                                                  >  1:238462/1‑150 (MQ=255)
                                                                                                gTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTAc                                                >  1:1019672/1‑151 (MQ=255)
                                                                                                  aaaTCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTAcc                                               <  2:183156/150‑1 (MQ=255)
                                                                                                  aaaTCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCg                                              >  2:965546/1‑151 (MQ=255)
                                                                                                  aaaTCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCg                                              <  2:856881/151‑1 (MQ=255)
                                                                                                           cgcgCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTa                                      <  1:796153/150‑1 (MQ=255)
                                                                                                              gcAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGa                                  <  1:956181/151‑1 (MQ=255)
                                                                                                               cAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGaa                                 <  2:899726/151‑1 (MQ=255)
                                                                                                               cAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGa                                  >  2:922580/1‑150 (MQ=255)
                                                                                                                       aaaTACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGctct                         <  1:921992/151‑1 (MQ=255)
                                                                                                                       aaaTACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGctct                         >  2:214049/1‑151 (MQ=255)
                                                                                                                        aaTACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGACCTGGCTCTg                        <  1:48219/151‑1 (MQ=255)
                                                                                                                          tACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGc                       >  1:744267/1‑150 (MQ=255)
                                                                                                                             ttCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCag                     >  2:941671/1‑149 (MQ=255)
                                                                                                                             ttCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCag                     >  1:957260/1‑149 (MQ=255)
                                                                                                                               cccGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGc                  <  1:686513/150‑1 (MQ=255)
                                                                                                                                 cGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGt               <  2:1059971/151‑1 (MQ=255)
                                                                                                                                  ggATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGt               >  1:773742/1‑150 (MQ=255)
                                                                                                                                  ggATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTc              >  2:404281/1‑151 (MQ=255)
                                                                                                                                     tGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTgg           >  2:521502/1‑151 (MQ=255)
                                                                                                                                       ccaccaTCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGca         <  1:205221/151‑1 (MQ=255)
                                                                                                                                        caccaTCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcac        <  1:695029/151‑1 (MQ=255)
                                                                                                                                         accaTCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcaca       <  1:190813/151‑1 (MQ=255)
                                                                                                                                         accaTCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcac        >  2:38548/1‑150 (MQ=255)
                                                                                                                                         accaTCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcac        >  1:957997/1‑150 (MQ=255)
                                                                                                                                              ccTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAg  <  1:250461/151‑1 (MQ=255)
                                                                                                                                                      ||                                                                                                                                             
ACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑TGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAG  >  NC_000913/1756743‑1757033

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