Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 361,673 A→T 8.5% V24E (GTG→GAG)  lacI ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273361,6730AT 8.5% 24608.1 / inf 8163V24E (GTG→GAG) lacIDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (3407/4058);  new base T (306/387);  total (3716/4447)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 8505 total aligned reads displayed.

AGGGCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGT  >  CP009273/361526‑361820
                                                                                                                                                   |                                                                                                                                                    
  ggCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTcacc                                                                                                                                                    >  1:356878/1‑148 (MQ=255)
  ggCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTcacc                                                                                                                                                    <  1:6189/148‑1 (MQ=255)
  ggCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTcacc                                                                                                                                                    >  2:1636486/1‑148 (MQ=255)
   gCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTTATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTcacca                                                                                                                                                   <  1:1202979/148‑1 (MQ=255)
      aGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCTCCAcgc                                                                                                                                                >  2:1284514/1‑148 (MQ=255)
        aCTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCTCCACGCgg                                                                                                                                              <  1:1142133/148‑1 (MQ=255)
          tGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGa                                                                                                                                           >  2:2366359/1‑148 (MQ=255)
             aGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAc                                                                                                                                        <  2:170331/148‑1 (MQ=255)
                tGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTCATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGt                                                                                                                                     <  2:1761048/148‑1 (MQ=255)
                tGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGt                                                                                                                                     <  2:3418782/148‑1 (MQ=255)
                 ggCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTc                                                                                                                                    <  1:3147492/148‑1 (MQ=255)
                    aaCGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGa                                                                                                                                 <  2:3795251/148‑1 (MQ=255)
                    aaCGCCAATCAGCAACGACTGTTTGCCCCCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAACGGTCTg                                                                                                                                  <  2:3490005/148‑1 (MQ=255)
                     aCGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAACGGTCTGa                                                                                                                                 <  1:1845432/148‑1 (MQ=255)
                     aCGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTGTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCTCCACGCGGGAAACGGTCTGAt                                                                                                                                <  2:2831859/148‑1 (MQ=255)
                      cGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCGGGTTCCCCACGCGGGAAACGGTCTGATa                                                                                                                               >  2:732653/1‑148 (MQ=255)
                       gCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAACGGTCTGATa                                                                                                                               >  1:3488637/1‑148 (MQ=255)
                         cAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAaga                                                                                                                            <  1:149838/148‑1 (MQ=255)
                             cAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGCAAACGGTCGGATAAGAGac                                                                                                                         >  1:1000396/1‑148 (MQ=255)
                             cAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGaca                                                                                                                        <  2:3023708/148‑1 (MQ=255)
                             cAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGaca                                                                                                                        <  1:2597016/148‑1 (MQ=255)
                               gCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACAcc                                                                                                                      <  1:2342410/148‑1 (MQ=255)
                                     aCTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCCTa                                                                                                                >  2:2696276/1‑148 (MQ=255)
                                     aCTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATa                                                                                                                >  1:1976254/1‑148 (MQ=255)
                                        gTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCGTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAACGGTCTGATAAGAGACACCGGCATAct                                                                                                              <  2:2284439/148‑1 (MQ=255)
                                         tgtGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATActct                                                                                                            <  2:2886381/146‑1 (MQ=255)
                                               cGCCAGTTGTTGTGCCACGCGGTTGGGAATGTTATTCAGCTCCGCCATCGCAGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACa                                                                                                      >  2:2765286/1‑148 (MQ=255)
                                               cGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGTAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACa                                                                                                      >  1:3761537/1‑148 (MQ=255)
                                                gCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACAt                                                                                                     <  1:2412590/148‑1 (MQ=255)
                                                   aGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGt                                                                                                  <  1:2740411/148‑1 (MQ=255)
                                                    gttgttGTGCCACGCGGTTGGGAATGTCATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTTGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGta                                                                                                 <  2:2429801/148‑1 (MQ=255)
                                                      tgttgtGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCCCTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCGTATAAGAGACACCGGCATACTCTGCGACCTCGgata                                                                                               >  2:3560506/1‑148 (MQ=255)
                                                       gtagtgCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATaa                                                                                              <  1:3070517/145‑1 (MQ=255)
                                                        ttgtGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAAc                                                                                             >  2:793174/1‑148 (MQ=255)
                                                          gtgCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACg                                                                                            <  2:1572271/148‑1 (MQ=255)
                                                          gtgCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGt                                                                                           <  2:3736019/148‑1 (MQ=255)
                                                             ccACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCAGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAAGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTAc                                                                                        <  1:2801505/148‑1 (MQ=255)
                                                               aCGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCTCCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTg                                                                                      <  2:201433/148‑1 (MQ=255)
                                                                cgcgGTTGGGAATGTAAGTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTgg                                                                                     <  1:2266815/148‑1 (MQ=255)
                                                                 gtggTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGGCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGt                                                                                    >  1:763852/3‑148 (MQ=255)
                                                                 gcgGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACTTTACTGGt                                                                                    <  2:610862/148‑1 (MQ=255)
                                                                     ttGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTca                                                                                >  1:997992/1‑148 (MQ=255)
                                                                     ttGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTca                                                                                >  1:541224/1‑148 (MQ=255)
                                                                     ttGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTca                                                                                >  1:3369337/1‑148 (MQ=255)
                                                                      tGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCTCGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTcac                                                                               <  2:5464/148‑1 (MQ=255)
                                                                      tGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTca                                                                                <  1:3228572/148‑1 (MQ=255)
                                                                         gAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCTCGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACAtt                                                                            >  2:256660/1‑148 (MQ=255)
                                                                         gAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACAtt                                                                            >  1:286689/1‑148 (MQ=255)
                                                                         gAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAAGGTCTGATAAGAGACACCGGCATACTCTGCGCCATCGTCTAACGTTACTGGTTTCACAtt                                                                            >  1:1915954/1‑148 (MQ=255)
                                                                            tGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGtttcactttcac                                                                         <  2:3093409/148‑1 (MQ=255)
                                                                            tGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTcac                                                                         >  1:395195/1‑148 (MQ=255)
                                                                             gTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTATCACATTcacc                                                                        >  1:2195334/1‑148 (MQ=255)
                                                                              tAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTcacca                                                                       <  2:2141995/148‑1 (MQ=255)
                                                                                aTTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCaccacc                                                                     <  1:845105/148‑1 (MQ=255)
                                                                                 ttCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCaccacc                                                                     <  2:3692923/148‑1 (MQ=255)
                                                                                  tCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCAccc                                                                    >  1:3006834/1‑148 (MQ=255)
                                                                                  tCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCt                                                                   >  1:2931957/1‑148 (MQ=255)
                                                                                   cAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAAGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTg                                                                  >  1:1284534/1‑148 (MQ=255)
                                                                                    aGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAAAGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGa                                                                 <  1:2678695/148‑1 (MQ=255)
                                                                                     gCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGaa                                                                >  1:3303309/1‑148 (MQ=255)
                                                                                        ccGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTg                                                             >  2:2627463/1‑148 (MQ=255)
                                                                                            cATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGActc                                                         >  2:1157511/1‑148 (MQ=255)
                                                                                             aTCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCAGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCAAACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGActct                                                        >  1:177637/1‑148 (MQ=255)
                                                                                              tCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGActctc                                                       <  1:1071441/148‑1 (MQ=255)
                                                                                               cgccgcTTCCACATTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGCAAACGGTCTGATAAGAGACACCCGCATACTCTGCGACATCGGATAACGTTACTGGTTTCACATTCACCACCCTGAATTTActctc                                                       >  2:664387/1‑148 (MQ=255)
                                                                                               cgccgcTTCCACATTTTCCCGCGTTTTCGCAGAAACGGGGCTGGCCAGGTTCACCACGCCGGAAACCGTCTTATAAGAGACACCGGCATACACTGCGACATCCTAAAAAGTAAATGGTTTCAAATTCACCACCCAGAATTGActctca                                                      >  2:2496328/1‑147 (MQ=255)
                                                                                                gccgcTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCTCGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCtt                                                     <  2:3328306/148‑1 (MQ=255)
                                                                                                 ccgcTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCTCCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTc                                                    <  2:1884648/148‑1 (MQ=255)
                                                                                                    cTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCgg                                                 >  1:2050413/1‑148 (MQ=255)
                                                                                                     ttCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCTCCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCggg                                                >  2:422456/1‑148 (MQ=255)
                                                                                                     ttCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCggg                                                >  1:3428233/1‑148 (MQ=255)
                                                                                                     ttCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCggg                                                >  1:694168/1‑148 (MQ=255)
                                                                                                        cACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGgcgc                                             <  2:2077612/148‑1 (MQ=255)
                                                                                                          cTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTa                                           <  1:1558708/148‑1 (MQ=255)
                                                                                                          cTTTGTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTa                                           <  2:869189/148‑1 (MQ=255)
                                                                                                           tttttCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTAt                                          >  2:122473/1‑148 (MQ=255)
                                                                                                            ttttCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGTTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATc                                         >  1:2480049/1‑148 (MQ=255)
                                                                                                               tCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATg                                      <  2:107012/148‑1 (MQ=255)
                                                                                                                cccGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGc                                     >  2:941788/1‑148 (MQ=255)
                                                                                                                cccGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGc                                     <  2:3260993/148‑1 (MQ=255)
                                                                                                                    cgTTTTCGCAGAAACGTGGCTGGCCTGGTTCTCCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATa                                 >  1:1359765/1‑148 (MQ=255)
                                                                                                                    cgTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATa                                 >  2:484757/1‑148 (MQ=255)
                                                                                                                      ttttCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATAcc                               >  1:2535257/1‑148 (MQ=255)
                                                                                                                              gAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCGATCATGCCATACCGCGAAAgg                       <  1:1715833/148‑1 (MQ=255)
                                                                                                                               aaaCGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGt                      >  2:3179833/1‑148 (MQ=255)
                                                                                                                                aaCGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGtt                     <  1:470749/148‑1 (MQ=255)
                                                                                                                                aaCGTGGCTGGCCTGGTTCACCACGCGGGAAAAGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGtt                     >  2:1436822/1‑148 (MQ=255)
                                                                                                                                 aCGTGGCTGGCCTGGTTCACCACGCGGTAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGttt                    >  1:162882/1‑148 (MQ=255)
                                                                                                                                  cGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGtttt                   >  1:3697503/1‑148 (MQ=255)
                                                                                                                                  cGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTATCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGtttt                   >  2:3633755/1‑148 (MQ=255)
                                                                                                                                    tggctggcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTgc                 >  1:3628373/1‑148 (MQ=255)
                                                                                                                                         ggcCTGGTTCTCCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAt            >  1:915758/1‑148 (MQ=255)
                                                                                                                                         ggcCTGGGTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAt            <  2:1962203/148‑1 (MQ=255)
                                                                                                                                          gcCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACATCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCAtt           >  2:2223091/1‑148 (MQ=255)
                                                                                                                                             tGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCCTTCGa        >  1:1427223/1‑148 (MQ=255)
                                                                                                                                             tGGTTCACCACGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGa        >  2:1500101/1‑148 (MQ=255)
                                                                                                                                               gTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATg      >  1:2116411/1‑148 (MQ=255)
                                                                                                                                                  cacctcGCGGGAAACGGTCTGATAAGAGACACCGGCAGACTCGGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGgtg   <  2:1830596/148‑1 (MQ=255)
                                                                                                                                                   accacGCGGGAAACGGTCTTATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGgtgt  <  2:136713/148‑1 (MQ=255)
                                                                                                                                                   accacGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTAACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCACCATTCGATGgtgt  <  2:296239/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGGCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGT  >  CP009273/361526‑361820

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