Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,526,796 A→T 8.8% Y207F (TAC→TTC)  pck → phosphoenolpyruvate carboxykinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,526,7960AT 8.8% 563.3 / 36.4 189Y207F (TAC→TTC) pckphosphoenolpyruvate carboxykinase
Reads supporting (aligned to +/- strand):  ref base A (90/82);  new base T (3/14);  total (93/96)
Fisher's exact test for biased strand distribution p-value = 9.30e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 190 total aligned reads displayed.

AAGAACTGGCAGGTTTCAAACCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCA  >  CP009273/3526649‑3526942
                                                                                                                                                   |                                                                                                                                                  
   aaCTGGCAGGTTTCAAACCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTAcgg                                                                                                                                                 >  2:2136628/1‑148 (MQ=255)
    aCTGGCAGGTTTCAAACCAGACTTTATCGTTATGAACGGCGCGCAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTAAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGATCGCATGCAGCTGATTGGCGGCACCTGGTTcggc                                                                                                                                                <  2:152709/148‑1 (MQ=255)
     cTGGCAGTTTTAAAACCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAATAACAGGGTCTAAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTAcggcg                                                                                                                                               <  2:1634886/148‑1 (MQ=255)
     cTGGCAGGTTTCAAACCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTAcggcg                                                                                                                                               <  1:1770797/148‑1 (MQ=255)
          aGGTTTCAAACCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGaa                                                                                                                                          <  1:3026009/148‑1 (MQ=255)
          aGGTTTCAAACCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGaa                                                                                                                                          <  1:673676/148‑1 (MQ=255)
              ttCAAACCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATga                                                                                                                                      <  1:1216134/148‑1 (MQ=255)
                cAAACCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATgaag                                                                                                                                    <  2:2246193/148‑1 (MQ=255)
                 aaaCCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATgaaga                                                                                                                                   <  1:1311723/148‑1 (MQ=255)
                 aaaCCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATgaaga                                                                                                                                   <  1:748024/148‑1 (MQ=255)
                 aaaCCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATgaaga                                                                                                                                   <  1:734138/148‑1 (MQ=255)
                  aaCCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATgaagaa                                                                                                                                  >  2:207094/1‑148 (MQ=255)
                   aCCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGaaa                                                                                                                                 >  2:894402/1‑148 (MQ=255)
                    ccAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAg                                                                                                                                <  1:1812973/148‑1 (MQ=255)
                     cAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAgg                                                                                                                               <  1:1485011/148‑1 (MQ=255)
                          tttATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGt                                                                                                                          <  2:883849/148‑1 (MQ=255)
                           ttATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGtt                                                                                                                         <  1:2111185/148‑1 (MQ=255)
                             aTCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTtct                                                                                                                       <  2:1645716/148‑1 (MQ=255)
                             aTCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTtct                                                                                                                       <  1:1214581/148‑1 (MQ=255)
                                gTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTTCGGCGGCGAAATGAAGAAAGGGATGTTCTCga                                                                                                                    <  2:906028/148‑1 (MQ=255)
                                gTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCga                                                                                                                    <  1:1380041/148‑1 (MQ=255)
                                  tATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTTCGGCGGCGAAATGAAGAAAGGGATGTTCTCgatg                                                                                                                  >  1:2544292/1‑148 (MQ=255)
                                  tATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCgatg                                                                                                                  <  2:2344059/148‑1 (MQ=255)
                                  tATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCgatg                                                                                                                  >  2:257318/1‑148 (MQ=255)
                                    tGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCgctgat                                                                                                                >  2:1336621/1‑148 (MQ=255)
                                     gAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGatgatg                                                                                                               >  1:2648528/1‑148 (MQ=255)
                                      aaCGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGAtgatga                                                                                                              >  1:1062799/1‑148 (MQ=255)
                                        cGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAAc                                                                                                            >  1:1075290/1‑148 (MQ=255)
                                        cGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAAc                                                                                                            >  1:1285076/1‑148 (MQ=255)
                                        cGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAAc                                                                                                            <  2:1507109/148‑1 (MQ=255)
                                         ggCGCGAAGTGCACTAACCCGCAGTTGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTTACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTTCGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACt                                                                                                           <  2:670356/148‑1 (MQ=255)
                                         ggCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACt                                                                                                           <  2:2925466/148‑1 (MQ=255)
                                            gcgAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTAcc                                                                                                        <  2:814524/148‑1 (MQ=255)
                                              gAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACctg                                                                                                      <  2:3183786/148‑1 (MQ=255)
                                               aaGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACctgc                                                                                                     <  1:1042837/148‑1 (MQ=255)
                                                       tAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTg                                                                                             >  2:3268425/1‑148 (MQ=255)
                                                       tAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTg                                                                                             >  1:2010015/1‑148 (MQ=255)
                                                         aCCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGaa                                                                                           >  2:3180204/1‑148 (MQ=255)
                                                          cccGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGaaa                                                                                          >  2:4939/1‑148 (MQ=255)
                                                            cGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAgg                                                                                        >  2:3145612/1‑148 (MQ=255)
                                                              cAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTa                                                                                      <  2:2408610/148‑1 (MQ=255)
                                                               aGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTAt                                                                                     <  2:1002363/148‑1 (MQ=255)
                                                               aGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTAt                                                                                     <  1:665699/148‑1 (MQ=255)
                                                                gTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATc                                                                                    >  1:449760/1‑148 (MQ=255)
                                                                   gAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCt                                                                                 >  2:2812875/1‑148 (MQ=255)
                                                                      aGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCCGCTGATTGGCGGCACCGGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTAGATGATGACCTCACTGCTGCCGCTGAAAGGTATCGCGTCt                                                                              >  1:3086995/1‑148 (MQ=255)
                                                                      aGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCt                                                                              >  1:1435052/1‑148 (MQ=255)
                                                                       gAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTTCGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTa                                                                             <  2:556997/148‑1 (MQ=255)
                                                                       gAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTTCGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTa                                                                             <  2:3203432/148‑1 (MQ=255)
                                                                        aaCAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGACGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTAt                                                                            >  2:1953255/1‑148 (MQ=255)
                                                                        aaCAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTAt                                                                            >  2:2859845/1‑148 (MQ=255)
                                                                        aaCAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTAt                                                                            >  2:369501/1‑148 (MQ=255)
                                                                        aaCAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTAt                                                                            >  2:431658/1‑148 (MQ=255)
                                                                         aCAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATg                                                                           >  1:1940874/1‑148 (MQ=255)
                                                                           aGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCa                                                                         <  1:716247/148‑1 (MQ=255)
                                                                           aGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCa                                                                         <  1:2565189/148‑1 (MQ=255)
                                                                            gggTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCAc                                                                        >  1:1235287/1‑148 (MQ=255)
                                                                                ctcAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTTCGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCt                                                                    <  1:1125246/148‑1 (MQ=255)
                                                                                ctcAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCt                                                                    >  1:938440/1‑148 (MQ=255)
                                                                                ctcAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCt                                                                    >  1:472903/1‑148 (MQ=255)
                                                                                   aaCTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCg                                                                 >  1:3097446/1‑148 (MQ=255)
                                                                                    aCTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGc                                                                >  2:1449834/1‑148 (MQ=255)
                                                                                    aCTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGc                                                                >  2:1476727/1‑148 (MQ=255)
                                                                                      tCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCa                                                              >  1:72424/1‑148 (MQ=255)
                                                                                      tCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCa                                                              >  1:2513995/1‑148 (MQ=255)
                                                                                      tCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCa                                                              >  2:59802/1‑148 (MQ=255)
                                                                                         gAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACg                                                           <  2:3175899/148‑1 (MQ=255)
                                                                                         gAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACg                                                           >  2:395509/1‑148 (MQ=255)
                                                                                             aCTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTgg                                                       <  1:3028379/148‑1 (MQ=255)
                                                                                              cTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGt                                                      <  1:894402/148‑1 (MQ=255)
                                                                                                    ggCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGaaa                                                <  2:353779/148‑1 (MQ=255)
                                                                                                           aaCCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATg                                         >  2:2428280/1‑148 (MQ=255)
                                                                                                            aCCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGt                                        <  2:1752609/148‑1 (MQ=255)
                                                                                                            aCCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGt                                        <  2:79012/148‑1 (MQ=255)
                                                                                                              cTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTg                                      <  2:1235287/148‑1 (MQ=255)
                                                                                                                gACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCg                                    <  1:2627895/148‑1 (MQ=255)
                                                                                                                    gAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGgtgt                                <  2:719878/148‑1 (MQ=255)
                                                                                                                     aGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGtt                               <  1:3179771/148‑1 (MQ=255)
                                                                                                                     aGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGtt                               <  1:2357934/148‑1 (MQ=255)
                                                                                                                        gcATGCAGCTGATTGGCGGCACCTGGTTCGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGttctt                            <  1:3082597/148‑1 (MQ=255)
                                                                                                                        gcATGCAGCTGATTGGCGGCACCTGGTTCGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGttctt                            <  1:329387/148‑1 (MQ=255)
                                                                                                                        gcATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGttctt                            <  2:673748/148‑1 (MQ=255)
                                                                                                                         cATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGttcttc                           <  2:116562/148‑1 (MQ=255)
                                                                                                                          aTGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCg                          >  1:157670/1‑148 (MQ=255)
                                                                                                                           tGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCgg                         >  1:2874995/1‑148 (MQ=255)
                                                                                                                            gCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGc                        >  1:949562/1‑148 (MQ=255)
                                                                                                                            gCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGc                        >  1:2603703/1‑148 (MQ=255)
                                                                                                                             cAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGcc                       <  2:449760/148‑1 (MQ=255)
                                                                                                                               gCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCtt                     >  1:485347/1‑148 (MQ=255)
                                                                                                                                cTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCttt                    <  1:2435622/148‑1 (MQ=255)
                                                                                                                                cTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCttt                    <  1:575028/148‑1 (MQ=255)
                                                                                                                                cTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGGATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCttt                    <  2:2780152/148‑1 (MQ=255)
                                                                                                                                   aTTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCg                 >  1:1465049/1‑148 (MQ=255)
                                                                                                                                   aTTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCg                 >  2:1059741/1‑148 (MQ=255)
                                                                                                                                       gcggcACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCAc             >  1:274275/1‑148 (MQ=255)
                                                                                                                                           cACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGt         >  2:2963566/1‑148 (MQ=255)
                                                                                                                                             ccTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTaa       >  1:1552769/1‑148 (MQ=255)
                                                                                                                                               tGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTaaaa     <  1:594543/148‑1 (MQ=255)
                                                                                                                                                 gTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAacc   >  1:85577/1‑148 (MQ=255)
                                                                                                                                                  tACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAacca  <  2:1154704/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAGAACTGGCAGGTTTCAAACCAGACTTTATCGTTATGAACGGCGCGAAGTGCACTAACCCGCAGTGGAAAGAACAGGGTCTCAACTCCGAAAACTTCGTGGCGTTTAACCTGACCGAGCGCATGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCA  >  CP009273/3526649‑3526942

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: