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.T93.9% 241.6 / 7.8 66coding (1195/1413 nt)pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base . (0/4);  new base T (34/28);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*NC_0009131,756,8922.G93.9% 256.4 / 8.7 66coding (1195/1413 nt)pykFpyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base . (0/4);  new base G (34/28);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑TGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAG  >  NC_000913/1756743‑1757033
                                                                                                                                                      ||                                                                                                                                             
aCCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑t                                                                                                                                              <  2:517494/151‑1 (MQ=255)
aCCGTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑t                                                                                                                                              <  1:399716/151‑1 (MQ=255)
   gTAAACTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑Tgtg                                                                                                                                           <  1:1046731/151‑3 (MQ=255)
     aaaCTGCGCATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑Tgtga                                                                                                                                          <  2:151387/150‑4 (MQ=255)
          gcgcATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑Tgtgaccacca                                                                                                                                    <  2:828761/151‑10 (MQ=255)
          gcgcATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCAC‑‑Tgtgaccacc                                                                                                                                     <  1:852242/150‑9 (MQ=255)
            gcATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAAc                                                                                                                                    <  1:101984/151‑1 (MQ=255)
            gcATTACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAAc                                                                                                                                    <  1:636092/151‑1 (MQ=255)
                tACCGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGaa                                                                                                                                 <  2:435845/150‑1 (MQ=255)
                   cGAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCAGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGaaaaaa                                                                                                                             <  2:1056504/151‑1 (MQ=255)
                    gAAGCGGTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAAc                                                                                                                            >  2:471463/1‑151 (MQ=255)
                         ggTATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCtc                                                                                                                       >  2:992730/1‑151 (MQ=255)
                           tATGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCtcat                                                                                                                     >  2:804507/1‑151 (MQ=255)
                            aTGCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCtcatc                                                                                                                    >  1:329779/1‑151 (MQ=255)
                              gCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCtcatca                                                                                                                   <  1:279292/150‑1 (MQ=255)
                              gCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCtcatca                                                                                                                   <  1:854700/150‑1 (MQ=255)
                              gCCGTGGTGCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCtcatc                                                                                                                    >  2:149057/1‑149 (MQ=255)
                                 gtggtgCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTg                                                                                                               <  1:706623/151‑1 (MQ=255)
                                     tgCCGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTAc                                                                                                           >  2:105431/1‑151 (MQ=255)
                                       ccGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTg                                                                                                         >  1:926576/1‑151 (MQ=255)
                                       ccGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTg                                                                                                         >  2:522040/1‑151 (MQ=255)
                                        cGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTg                                                                                                         >  1:495809/1‑150 (MQ=255)
                                        cGTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGa                                                                                                        <  2:264292/151‑1 (MQ=255)
                                         gTTGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAg                                                                                                       <  2:561627/151‑1 (MQ=255)
                                           tGAAACTGCTGAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGc                                                                                                      >  2:499518/1‑150 (MQ=255)
                                                  gctgAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGc                                                                                                <  2:186480/149‑1 (MQ=255)
                                                  gctgAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGc                                                                                                <  2:838103/149‑1 (MQ=255)
                                                  gctgAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGt                                                                                              >  2:155752/1‑151 (MQ=255)
                                                     gAAAAACTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTtgt                                                                                           >  1:138566/1‑151 (MQ=255)
                                                      aaaaaCTGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTtgtg                                                                                          <  2:200446/151‑1 (MQ=255)
                                                            tGGATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAGAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGc                                                                                      >  2:611069/1‑149 (MQ=255)
                                                             ggATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAg                                                                                    <  2:478261/150‑1 (MQ=255)
                                                             ggATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAg                                                                                    <  1:607995/150‑1 (MQ=255)
                                                             ggATGCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAg                                                                                    <  1:441185/150‑1 (MQ=255)
                                                                 gCTCCGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGt                                                                               >  2:61133/1‑151 (MQ=255)
                                                                    ccGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTa                                                                             >  2:563479/1‑150 (MQ=255)
                                                                    ccGCTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTa                                                                             <  2:824159/150‑1 (MQ=255)
                                                                       cTGATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTaa                                                                            >  1:980568/1‑148 (MQ=255)
                                                                         gATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAagag                                                                        <  1:485467/150‑1 (MQ=255)
                                                                         gATCGTGGTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAag                                                                          >  2:82194/1‑148 (MQ=255)
                                                                                gTTGCTACTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACtt                                                                 >  1:163724/1‑150 (MQ=255)
                                                                                      aCTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACt                                                            <  1:598632/149‑1 (MQ=255)
                                                                                       cTCAGGGCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACtgat                                                         >  2:431626/1‑151 (MQ=255)
                                                                                             gCGGTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTc                                                   <  1:70843/151‑1 (MQ=255)
                                                                                               ggTAAATCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTa                                                 <  1:599173/151‑1 (MQ=255)
                                                                                                  aaaTCTGCTCGCGCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCg                                              >  2:825461/1‑151 (MQ=255)
                                                                                                           cgcgCAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGt                                       <  2:841718/149‑1 (MQ=255)
                                                                                                             cgcAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTaa                                     >  2:156603/1‑149 (MQ=255)
                                                                                                              gcAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTaa                                     >  2:54473/1‑148 (MQ=255)
                                                                                                              gcAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTa                                      >  1:655134/1‑147 (MQ=255)
                                                                                                               cAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGaa                                 <  1:1763/151‑1 (MQ=255)
                                                                                                               cAGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGa                                  >  2:1039761/1‑150 (MQ=255)
                                                                                                                aGTACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAAc                                >  1:865187/1‑151 (MQ=255)
                                                                                                                  tACGTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTg                              >  1:668446/1‑151 (MQ=255)
                                                                                                                     gTAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGct                           >  2:731876/1‑151 (MQ=255)
                                                                                                                      tAAATACTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGctc                          >  1:621433/1‑151 (MQ=255)
                                                                                                                           aCTTCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGc                       >  1:134424/1‑149 (MQ=255)
                                                                                                                              tCCCGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGc                  >  1:428315/1‑151 (MQ=255)
                                                                                                                               cccGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAACACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCg                 >  2:1049880/1‑151 (MQ=255)
                                                                                                                               cccGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGc                  <  1:896385/150‑1 (MQ=255)
                                                                                                                                 cGGATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGt               >  2:68841/1‑151 (MQ=255)
                                                                                                                                   gATGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCt             <  2:891558/151‑1 (MQ=255)
                                                                                                                                    aTGCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTg            <  2:18351/151‑1 (MQ=255)
                                                                                                                                      gCCACCATCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGc          <  2:11625/151‑1 (MQ=255)
                                                                                                                                         accaTCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcac        >  1:352520/1‑150 (MQ=255)
                                                                                                                                         accaTCCTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGcac        >  2:26605/1‑150 (MQ=255)
                                                                                                                                              ccTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACaaa   <  2:803378/150‑1 (MQ=255)
                                                                                                                                              ccTGGCACTGTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAg  <  1:125990/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 < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: