Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,933,166 C→A 5.0% R421S (CGT→AGT)  pykA → pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,933,1660CA 5.0% 991.0 / 15.8 255R421S (CGT→AGT) pykApyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base C (128/114);  new base A (9/4);  total (137/118)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 255 total aligned reads displayed.

CAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGT  >  CP009273/1933020‑1933312
                                                                                                                                                  |                                                                                                                                                  
cAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCg                                                                                                                                                   <  1:1555744/148‑1 (MQ=255)
  gCTAACCTCCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTg                                                                                                                                                 >  2:2581942/1‑148 (MQ=255)
  gCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTg                                                                                                                                                 >  2:1385028/1‑148 (MQ=255)
  gCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTg                                                                                                                                                 >  1:930431/1‑148 (MQ=255)
     aaCCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCTTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCg                                                                                                                                              <  1:332833/148‑1 (MQ=255)
     aaCCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCg                                                                                                                                              <  2:2000533/148‑1 (MQ=255)
       ccaccTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGtt                                                                                                                                            >  2:536421/1‑148 (MQ=255)
         accTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTAc                                                                                                                                          >  2:1466202/1‑148 (MQ=255)
             gAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCg                                                                                                                                      <  2:1674319/148‑1 (MQ=255)
              aaaGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCgg                                                                                                                                     >  1:613242/1‑148 (MQ=255)
              aaaGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCgg                                                                                                                                     <  2:430276/148‑1 (MQ=255)
                aGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGAGTGCTCTCTATCGTGGCGTTACGCCGGTg                                                                                                                                   <  2:3253750/148‑1 (MQ=255)
                aGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTg                                                                                                                                   <  2:3253478/148‑1 (MQ=255)
                 ggCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGc                                                                                                                                  >  1:831875/1‑148 (MQ=255)
                  gCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCa                                                                                                                                 >  1:586876/1‑148 (MQ=255)
                     ttACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACtt                                                                                                                              <  1:681694/148‑1 (MQ=255)
                      tACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACttt                                                                                                                             <  1:2558241/148‑1 (MQ=255)
                      tACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACttt                                                                                                                             <  1:1914552/148‑1 (MQ=255)
                      tACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACttt                                                                                                                             >  2:3156649/1‑148 (MQ=255)
                             atcatcACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAgcg                                                                                                                      <  1:300699/148‑1 (MQ=255)
                                atcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTa                                                                                                                   <  1:2301566/148‑1 (MQ=255)
                                 tcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTaa                                                                                                                  >  1:1780212/1‑148 (MQ=255)
                                 tcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTaa                                                                                                                  <  2:1076179/148‑1 (MQ=255)
                                 tcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTaa                                                                                                                  <  2:3067181/148‑1 (MQ=255)
                                 tcaCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTaa                                                                                                                  <  2:2201570/148‑1 (MQ=255)
                                   aCCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATg                                                                                                                <  2:3416274/148‑1 (MQ=255)
                                    ccATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGa                                                                                                               >  2:1201867/1‑148 (MQ=255)
                                     cATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGAc                                                                                                              >  1:1663708/1‑148 (MQ=255)
                                       tGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACgg                                                                                                            <  1:1037261/148‑1 (MQ=255)
                                       tGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACgg                                                                                                            <  2:2749418/148‑1 (MQ=255)
                                       tGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACgg                                                                                                            <  1:2486973/148‑1 (MQ=255)
                                         aCCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCg                                                                                                          >  1:1221642/1‑148 (MQ=255)
                                         aCCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCg                                                                                                          <  1:629959/148‑1 (MQ=255)
                                         aCCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCg                                                                                                          >  1:22808/1‑148 (MQ=255)
                                           cGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTa                                                                                                        >  2:1236254/1‑148 (MQ=255)
                                             aaTCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTagc                                                                                                      <  1:3490507/148‑1 (MQ=255)
                                             aaTCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTagc                                                                                                      >  2:607228/1‑148 (MQ=255)
                                                  ggTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTg                                                                                                 <  2:2391002/148‑1 (MQ=255)
                                                  ggTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTGTGATAGCGCTAATGACGGCGTAGCAGCTg                                                                                                 >  2:183281/1‑148 (MQ=255)
                                                    tCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGcc                                                                                               >  1:426861/1‑148 (MQ=255)
                                                        aCCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCg                                                                                           <  1:3497050/148‑1 (MQ=255)
                                                        aCCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTAGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCTTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTCGCAGCTGCCAgag                                                                                           >  2:1770596/1‑146 (MQ=255)
                                                         ccGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGa                                                                                          <  1:944810/148‑1 (MQ=255)
                                                          cgcgCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGaa                                                                                         <  1:1839450/148‑1 (MQ=255)
                                                           gcgcTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAg                                                                                        <  1:1719811/148‑1 (MQ=255)
                                                             gcTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCg                                                                                      >  1:2870495/1‑148 (MQ=255)
                                                               tgatgaCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCAATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGt                                                                                    >  2:918013/1‑148 (MQ=255)
                                                                gatgaCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGtt                                                                                   <  2:248691/148‑1 (MQ=255)
                                                                gatgaCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGtt                                                                                   >  2:366581/1‑148 (MQ=255)
                                                                  tgaCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTaa                                                                                 <  2:279937/148‑1 (MQ=255)
                                                                    aCCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATc                                                                               <  2:2554321/148‑1 (MQ=255)
                                                                     ccTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATct                                                                              >  2:1425802/1‑148 (MQ=255)
                                                                        cccGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATctgct                                                                           >  2:1028577/1‑148 (MQ=255)
                                                                         ccGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATctgctg                                                                          <  1:19233/148‑1 (MQ=255)
                                                                           gTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTgcg                                                                        <  1:1378980/148‑1 (MQ=255)
                                                                            tATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGGGGTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTgcgc                                                                       >  2:2090610/1‑148 (MQ=255)
                                                                             aTCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGcgcg                                                                      <  2:2469276/148‑1 (MQ=255)
                                                                             aTCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGcgcg                                                                      >  2:730776/1‑148 (MQ=255)
                                                                             aTCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGcgcg                                                                      >  2:3436842/1‑148 (MQ=255)
                                                                                   tctGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAg                                                                <  2:709521/148‑1 (MQ=255)
                                                                                     tGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTAGCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGt                                                              <  2:3313732/148‑1 (MQ=255)
                                                                                        tCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTAc                                                           >  2:1527040/1‑148 (MQ=255)
                                                                                         cTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACt                                                          >  1:2141478/1‑148 (MQ=255)
                                                                                           gCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTCCTTg                                                        >  1:871454/1‑148 (MQ=255)
                                                                                           gCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTg                                                        >  1:3473003/1‑148 (MQ=255)
                                                                                             cAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGAt                                                      >  2:1686390/1‑148 (MQ=255)
                                                                                             cAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGAt                                                      >  2:3447800/1‑148 (MQ=255)
                                                                                               aTTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGt                                                    <  1:2699707/148‑1 (MQ=255)
                                                                                               aTTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGt                                                    <  1:2699102/148‑1 (MQ=255)
                                                                                                  ttCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTg                                                 >  1:146466/1‑148 (MQ=255)
                                                                                                    cGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGt                                               >  1:616903/1‑148 (MQ=255)
                                                                                                     gCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTg                                              >  2:3482665/1‑148 (MQ=255)
                                                                                                            cgcgCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGt                                       >  2:2157917/1‑148 (MQ=255)
                                                                                                              cgcCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGa                                     <  2:656786/148‑1 (MQ=255)
                                                                                                              cgcCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGa                                     >  1:1430352/1‑148 (MQ=255)
                                                                                                               gcCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGAt                                    <  2:1797962/148‑1 (MQ=255)
                                                                                                                ccATGAACGTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGAtt                                   >  2:1592251/1‑148 (MQ=255)
                                                                                                                    gAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCa                               <  1:2833254/148‑1 (MQ=255)
                                                                                                                    gAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCa                               <  2:3366897/148‑1 (MQ=255)
                                                                                                                       cGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCAccc                            >  2:1083130/1‑148 (MQ=255)
                                                                                                                       cGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCAccc                            >  2:376827/1‑148 (MQ=255)
                                                                                                                        gTACGCTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCa                           <  2:2238369/148‑1 (MQ=255)
                                                                                                                          aCGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAgg                         >  1:1186209/1‑148 (MQ=255)
                                                                                                                           cGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAggg                        <  2:818026/148‑1 (MQ=255)
                                                                                                                             cTGAACCTGACTGCTCTCTATAGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCg                      >  1:2806184/1‑148 (MQ=255)
                                                                                                                              tGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGa                     >  1:1071656/1‑148 (MQ=255)
                                                                                                                                 aCCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGt                  <  2:1627767/148‑1 (MQ=255)
                                                                                                                                 aCCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGt                  <  2:586876/148‑1 (MQ=255)
                                                                                                                                   cTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGtga                <  2:831875/148‑1 (MQ=255)
                                                                                                                                    tGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGtgat               >  2:2311284/1‑148 (MQ=255)
                                                                                                                                    tGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGtgat               <  2:1254617/148‑1 (MQ=255)
                                                                                                                                    tGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGtgat               >  1:1395014/1‑148 (MQ=255)
                                                                                                                                         gCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTa          >  1:180310/1‑148 (MQ=255)
                                                                                                                                          ctctctATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTAc         <  2:1780212/148‑1 (MQ=255)
                                                                                                                                             tctATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGt      >  1:303241/1‑148 (MQ=255)
                                                                                                                                              ctATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTg     <  2:2182519/148‑1 (MQ=255)
                                                                                                                                                aTCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTggg   <  2:3018782/148‑1 (MQ=255)
                                                                                                                                                aTCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTggg   >  1:278939/1‑148 (MQ=255)
                                                                                                                                                aTCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTggg   <  1:1249316/148‑1 (MQ=255)
                                                                                                                                                 tCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGt  >  1:774659/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGCTAACCACCTGAAAGGCGTTACGGCGATCATCACCATGACCGAATCGGGTCGTACCGCGCTGATGACCTCCCGTATCAGCTCTGGTCTGCCAATTTTCGCCATGTCGCGCCATGAACGTACGCTGAACCTGACTGCTCTCTATCGTGGCGTTACGCCGGTGCACTTTGATAGCGCTAATGACGGCGTAGCAGCTGCCAGCGAAGCGGTTAATCTGCTGCGCGATAAAGGTTACTTGATGTCTGGTGACCTGGTGATTGTCACCCAGGGCGACGTGATGAGTACCGTGGGT  >  CP009273/1933020‑1933312

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