Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,890,151 G→T 38.8% A6A (GCG→GCT yieH → phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,890,1510GT38.8% 301.3 / 245.1 247A6A (GCG→GCTyieHphosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (58/93);  new base T (57/39);  total (115/132)
Fisher's exact test for biased strand distribution p-value = 1.65e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 251 total aligned reads displayed.

GCTATCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGAAATTATCG  >  CP009273/3890003‑3890299
                                                                                                                                                    |                                                                                                                                                    
gCTATCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCg                                                                                                                                                      >  1:1938978/1‑149 (MQ=255)
 cTATCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCgg                                                                                                                                                     <  2:2469872/149‑1 (MQ=255)
  tATCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAATCGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGt                                                                                                                                                    >  2:3567432/1‑149 (MQ=255)
   aTCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTa                                                                                                                                                   <  1:1210201/149‑1 (MQ=255)
    tCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTAt                                                                                                                                                  >  2:1773429/1‑149 (MQ=255)
    tCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTAt                                                                                                                                                  >  2:3815176/1‑149 (MQ=255)
      gTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTAttt                                                                                                                                                >  2:406955/1‑149 (MQ=255)
      gTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTAttt                                                                                                                                                <  1:1984970/149‑1 (MQ=255)
      gTATACCAAAATGCGACAATAGGCGCGGTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTAttt                                                                                                                                                >  2:3742406/1‑149 (MQ=255)
        ataCCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTAttttt                                                                                                                                              <  2:3952755/149‑1 (MQ=255)
         taCCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTc                                                                                                                                             >  2:804293/1‑149 (MQ=255)
           ccAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGa                                                                                                                                           <  2:2831744/149‑1 (MQ=255)
            cAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGAc                                                                                                                                          <  2:1203728/149‑1 (MQ=255)
            cAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGAc                                                                                                                                          <  1:2749594/149‑1 (MQ=255)
            cAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGAc                                                                                                                                          <  2:1894849/149‑1 (MQ=255)
            cAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGAc                                                                                                                                          <  1:1466911/149‑1 (MQ=255)
                aTGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCg                                                                                                                                      >  2:2890517/1‑149 (MQ=255)
                      cAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTa                                                                                                                                <  2:919817/149‑1 (MQ=255)
                        aTAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATAAGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGTAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACg                                                                                                                              >  2:2643625/1‑149 (MQ=255)
                         tAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGc                                                                                                                             <  1:1321837/149‑1 (MQ=255)
                          aGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCt                                                                                                                            >  2:4017112/1‑149 (MQ=255)
                           ggCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTg                                                                                                                           >  2:3233783/1‑149 (MQ=255)
                             cgcgTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGt                                                                                                                         <  1:442740/149‑1 (MQ=255)
                              gcgTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTc                                                                                                                        <  1:2444831/149‑1 (MQ=255)
                                 tttGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGAc                                                                                                                     >  1:368847/1‑149 (MQ=255)
                                 tttGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGAc                                                                                                                     >  1:3810391/1‑149 (MQ=255)
                                  ttGTGATAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAACTCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACa                                                                                                                    <  2:1017900/149‑1 (MQ=255)
                                  ttGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACa                                                                                                                    >  1:3880552/1‑149 (MQ=255)
                                  ttGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACa                                                                                                                    <  2:1957653/149‑1 (MQ=255)
                                      gagagaCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGa                                                                                                                >  2:3887307/1‑149 (MQ=255)
                                       agagaCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGaa                                                                                                               <  1:1269535/149‑1 (MQ=255)
                                        gagaCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAg                                                                                                              >  1:645134/1‑149 (MQ=255)
                                          gaCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTc                                                                                                            >  1:3977017/1‑149 (MQ=255)
                                           aCTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCa                                                                                                           >  1:985248/1‑149 (MQ=255)
                                            cTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCAt                                                                                                          <  2:3014717/149‑1 (MQ=255)
                                             ttttttATTGGATTTACTTATACTTCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCAtt                                                                                                         <  1:3603182/149‑1 (MQ=255)
                                             ttttttATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCAtt                                                                                                         >  1:2354048/1‑149 (MQ=255)
                                                   aTTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGctct                                                                                                   <  1:3548444/149‑1 (MQ=255)
                                                   aTTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGctct                                                                                                   <  1:1677814/149‑1 (MQ=255)
                                                    ttGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCtctc                                                                                                  >  1:2189687/1‑149 (MQ=255)
                                                     tGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTcg                                                                                                 >  2:348048/1‑149 (MQ=255)
                                                       gATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTcgcg                                                                                               >  1:1484609/1‑149 (MQ=255)
                                                         tttACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCa                                                                                             >  1:4085599/1‑149 (MQ=255)
                                                          ttACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCat                                                                                            <  2:2823159/149‑1 (MQ=255)
                                                           tACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCata                                                                                           <  2:2053449/149‑1 (MQ=255)
                                                            aCTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCatat                                                                                          >  2:2160337/1‑149 (MQ=255)
                                                                ataCGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTaa                                                                                      >  1:1314788/1‑149 (MQ=255)
                                                                  aCGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACg                                                                                    <  2:11836/149‑1 (MQ=255)
                                                                    ggCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGAt                                                                                  <  1:2514810/149‑1 (MQ=255)
                                                                    ggCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGAt                                                                                  <  1:3374517/149‑1 (MQ=255)
                                                                     gCAACGATTGCGTTGCGCAAATCGGCTCTACGTAAACGTTACTCTGATTAAAAAGGAAAGGCGTGTCCCGGCTAGACGCTGTCTTTTTCGACTGCGCCGGTACGCTGGTCGAGAGTGTAGTCATGTGCTCTCGCGCATCTGTAACGATg                                                                                 >  2:2376975/1‑149 (MQ=255)
                                                                     gCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCCTATGTAACGATg                                                                                 >  2:2377036/1‑149 (MQ=255)
                                                                     gCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATg                                                                                 <  2:1758928/149‑1 (MQ=255)
                                                                     gCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATg                                                                                 >  2:2800089/1‑149 (MQ=255)
                                                                      cAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGt                                                                                <  1:591611/149‑1 (MQ=255)
                                                                      cAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGt                                                                                <  1:592359/149‑1 (MQ=255)
                                                                         cGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTc                                                                             <  1:471766/149‑1 (MQ=255)
                                                                            ttgcgttgcgCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAgg                                                                          <  2:1350267/149‑1 (MQ=255)
                                                                               cgttgcgCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAAt                                                                       <  1:4093846/149‑1 (MQ=255)
                                                                                gttgcgCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAAtt                                                                      <  1:2001275/149‑1 (MQ=255)
                                                                                 ttgcgCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAAttt                                                                     >  2:2944906/1‑149 (MQ=255)
                                                                                     gcAACTCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTa                                                                 <  2:3333707/149‑1 (MQ=255)
                                                                                       aaaTCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTAtt                                                               <  1:2227171/149‑1 (MQ=255)
                                                                                        aaTCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACCGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTa                                                              <  2:3166938/149‑1 (MQ=255)
                                                                                        aaTCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTa                                                              <  2:3630205/149‑1 (MQ=255)
                                                                                         aTCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTAc                                                             <  2:3145391/149‑1 (MQ=255)
                                                                                                ataagAAAAAGTTAAAGTGATTAAAAAGGAAAGGCATGTACCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGAt                                                      <  1:764729/145‑1 (MQ=255)
                                                                                                 tACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATc                                                     >  1:1586840/1‑149 (MQ=255)
                                                                                                   cGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCt                                                   >  1:1883933/1‑149 (MQ=255)
                                                                                                   cGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCt                                                   >  1:3655266/1‑149 (MQ=255)
                                                                                                    gAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTg                                                  <  1:1834146/149‑1 (MQ=255)
                                                                                                       aaCGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAag                                               >  1:3410179/1‑149 (MQ=255)
                                                                                                        aCGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAaga                                              <  1:2377245/149‑1 (MQ=255)
                                                                                                         cGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAagag                                             >  2:2270532/1‑149 (MQ=255)
                                                                                                           ttAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGt                                           <  1:3899983/149‑1 (MQ=255)
                                                                                                           ttAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGt                                           <  2:1982733/149‑1 (MQ=255)
                                                                                                                  gATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCaaa                                    >  1:1629000/1‑149 (MQ=255)
                                                                                                                   aTTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAAc                                   <  1:3291431/149‑1 (MQ=255)
                                                                                                                   aTTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAAc                                   >  1:1428928/1‑149 (MQ=255)
                                                                                                                   aTTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAAc                                   <  2:1938699/149‑1 (MQ=255)
                                                                                                                    ttAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACg                                  >  2:782540/1‑149 (MQ=255)
                                                                                                                          aggaaaggCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCaa                            >  1:3189062/1‑149 (MQ=255)
                                                                                                                          aggaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCaa                            >  2:2427924/1‑149 (MQ=255)
                                                                                                                           ggaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCaaa                           >  1:1143628/1‑149 (MQ=255)
                                                                                                                           ggaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCaaa                           <  2:1446230/149‑1 (MQ=255)
                                                                                                                            gaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAg                          <  1:3196751/149‑1 (MQ=255)
                                                                                                                            gaaaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAg                          <  1:1048601/149‑1 (MQ=255)
                                                                                                                              aaggCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGgt                        <  1:927379/149‑1 (MQ=255)
                                                                                                                               aggCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGgtg                       <  1:3335562/149‑1 (MQ=255)
                                                                                                                               aggCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGgtg                       >  2:44144/1‑149 (MQ=255)
                                                                                                                               aggCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGgtg                       >  2:44421/1‑149 (MQ=255)
                                                                                                                                ggCATGTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGgtgt                      >  1:4010849/1‑149 (MQ=255)
                                                                                                                                    tGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTaaaa                  <  2:2032577/149‑1 (MQ=255)
                                                                                                                                     gTCCCGGATAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAAc                 <  1:2914689/149‑1 (MQ=255)
                                                                                                                                         cGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGCTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTa             <  1:3847263/149‑1 (MQ=255)
                                                                                                                                            aTAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGa          <  2:1605168/149‑1 (MQ=255)
                                                                                                                                            aTAGAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGa          <  2:1605062/149‑1 (MQ=255)
                                                                                                                                               gAAGCTGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGAAAt       >  1:3213206/1‑149 (MQ=255)
                                                                                                                                               gAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGAAAt       <  2:592980/149‑1 (MQ=255)
                                                                                                                                                    ggTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGAAATTATCg  >  2:471943/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTATCGTATACCAAAATGCGACAATAGGCGCGTTTGTGAGAGACTTTTTTATTGGATTTACTTATACGGCAACGATTGCGTTGCGCAAATCGGCTCTACGAAAACGTTAAACTGATTAAAAAGGAAAGGCATGTCCCGGATAGAAGCGGTATTTTTCGACTGCGACGGTACGCTGGTCGACAGTGAAGTCATTTGCTCTCGCGCATATGTAACGATGTTTCAGGAATTTGGTATTACGCTCGATCCTGAAGAGGTATTCAAACGTTTCAAAGGTGTAAAACTGTACGAAATTATCG  >  CP009273/3890003‑3890299

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: