Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,894,334 C→A 5.0% G128G (GGC→GGA yqcE → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,894,3340CA 5.0% 997.5 / 17.3 258G128G (GGC→GGAyqcEMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base C (118/127);  new base A (9/4);  total (127/131)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 263 total aligned reads displayed.

TTCCGCGATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTAT  >  CP009273/2894187‑2894481
                                                                                                                                                   |                                                                                                                                                   
ttCCGCGATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGc                                                                                                                                                     >  1:3310571/1‑148 (MQ=255)
   cgcgATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGAt                                                                                                                                                  >  1:521514/1‑148 (MQ=255)
    gcgATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATc                                                                                                                                                 <  2:238433/148‑1 (MQ=255)
        tgatCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTTCATCGTTGCTTGGCGATCATAg                                                                                                                                             >  1:3402824/1‑148 (MQ=255)
        tgatCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAg                                                                                                                                             >  1:2392811/1‑148 (MQ=255)
        tgatCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAg                                                                                                                                             <  2:3244702/148‑1 (MQ=255)
          atCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCg                                                                                                                                           <  2:3007227/148‑1 (MQ=255)
          atCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCg                                                                                                                                           <  1:1754443/148‑1 (MQ=255)
           tCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGa                                                                                                                                          >  1:2609759/1‑148 (MQ=255)
              ttACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCa                                                                                                                                       >  2:2073072/1‑148 (MQ=255)
              ttACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCa                                                                                                                                       <  1:2306875/148‑1 (MQ=255)
                cccGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCGGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAg                                                                                                                                     >  2:723146/2‑148 (MQ=255)
                aCCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAg                                                                                                                                     <  1:7704/148‑1 (MQ=255)
                 ccGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAgg                                                                                                                                    >  2:542948/1‑148 (MQ=255)
                  cGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAggg                                                                                                                                   <  2:3271411/148‑1 (MQ=255)
                   ggATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGa                                                                                                                                  <  1:2990065/148‑1 (MQ=255)
                     aTTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGaaa                                                                                                                                >  1:2973565/1‑148 (MQ=255)
                     aTTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGAGATCATAGCGAGCAAGGGaaa                                                                                                                                >  1:2704248/1‑148 (MQ=255)
                      ttACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGaaaa                                                                                                                               <  2:1884368/148‑1 (MQ=255)
                        aCTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAAtt                                                                                                                             <  2:12470/148‑1 (MQ=255)
                        aCTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAAtt                                                                                                                             >  2:794391/1‑148 (MQ=255)
                        aCTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAAtt                                                                                                                             <  2:2164985/148‑1 (MQ=255)
                           gggTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATg                                                                                                                          <  2:1983811/148‑1 (MQ=255)
                            ggTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAAGTATgg                                                                                                                         >  2:2253467/1‑148 (MQ=255)
                                 gTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCtgg                                                                                                                    >  2:3792622/1‑148 (MQ=255)
                                  ttAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCtgga                                                                                                                   <  1:726752/148‑1 (MQ=255)
                                   tAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCtggat                                                                                                                  >  1:2625197/1‑148 (MQ=255)
                                     aTGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCtggatgg                                                                                                                <  1:3799395/148‑1 (MQ=255)
                                      tGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCtggatgga                                                                                                               <  1:2438193/148‑1 (MQ=255)
                                       ggCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGaa                                                                                                              <  2:1242524/148‑1 (MQ=255)
                                        gCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAg                                                                                                             >  1:3708631/1‑148 (MQ=255)
                                           aCGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGc                                                                                                          >  1:888864/1‑148 (MQ=255)
                                            cGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCt                                                                                                         <  1:2220840/148‑1 (MQ=255)
                                                         gtgGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGgt                                                                                            <  1:2400468/148‑1 (MQ=255)
                                                         gtgGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGgt                                                                                            <  2:2539458/148‑1 (MQ=255)
                                                          tgGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGgtg                                                                                           <  2:2706559/148‑1 (MQ=255)
                                                          tgGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGgtg                                                                                           <  1:1517466/148‑1 (MQ=255)
                                                             gTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTaa                                                                                        <  2:1883586/148‑1 (MQ=255)
                                                             gTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTaa                                                                                        >  1:3348409/1‑148 (MQ=255)
                                                             gTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTaa                                                                                        <  1:3634806/148‑1 (MQ=255)
                                                             gTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTaa                                                                                        <  1:1560710/148‑1 (MQ=255)
                                                             gTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTaa                                                                                        <  2:650099/148‑1 (MQ=255)
                                                                aTGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGt                                                                                     >  2:13646/1‑148 (MQ=255)
                                                                 tGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTc                                                                                    >  1:1379589/1‑148 (MQ=255)
                                                                 tGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTc                                                                                    >  2:2772877/1‑148 (MQ=255)
                                                                  gCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGAGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCg                                                                                   >  2:3791258/1‑148 (MQ=255)
                                                                   ctctGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGAGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGc                                                                                  >  1:2107591/1‑148 (MQ=255)
                                                                     ctGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTg                                                                                >  2:144224/1‑148 (MQ=255)
                                                                         aTTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCgg                                                                            >  1:1228890/1‑148 (MQ=255)
                                                                           tCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGgtg                                                                          >  2:3133073/1‑148 (MQ=255)
                                                                           tCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCATCGGTGTAATGTCGCTGGCGgtg                                                                          >  2:2563556/1‑148 (MQ=255)
                                                                              gATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGttt                                                                       >  2:1152844/1‑148 (MQ=255)
                                                                                 cGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTAcc                                                                    >  2:3329946/1‑148 (MQ=255)
                                                                                 cGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTAcc                                                                    <  1:2032470/148‑1 (MQ=255)
                                                                                  gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCa                                                                   <  1:518478/148‑1 (MQ=255)
                                                                                  gCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCa                                                                   <  1:3792608/148‑1 (MQ=255)
                                                                                   ccTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCAt                                                                  >  2:1791061/1‑148 (MQ=255)
                                                                                    cTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCAtg                                                                 >  2:2562744/1‑148 (MQ=255)
                                                                                      ttGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCAtgtg                                                               <  1:2403843/148‑1 (MQ=255)
                                                                                      ttGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCAtgtg                                                               >  2:3128338/1‑148 (MQ=255)
                                                                                          gATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTc                                                           <  2:3501771/148‑1 (MQ=255)
                                                                                          gATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTc                                                           <  2:312061/148‑1 (MQ=255)
                                                                                           aTAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCt                                                          <  2:384620/148‑1 (MQ=255)
                                                                                           aTAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCt                                                          <  2:922520/148‑1 (MQ=255)
                                                                                           aTAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCt                                                          <  2:84956/148‑1 (MQ=255)
                                                                                            tAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCtt                                                         <  1:304751/148‑1 (MQ=255)
                                                                                             aaCGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCttt                                                        <  1:1095109/148‑1 (MQ=255)
                                                                                             aaCGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCttt                                                        >  1:785612/1‑148 (MQ=255)
                                                                                              acgacgATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGAGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCtttt                                                       <  2:2201582/148‑1 (MQ=255)
                                                                                                    aTTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCt                                                 >  2:435993/1‑148 (MQ=255)
                                                                                                     ttttAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCTCtt                                                >  2:3659160/1‑148 (MQ=255)
                                                                                                      tttAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCttt                                               >  2:1994352/1‑148 (MQ=255)
                                                                                                      tttAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGAGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCttt                                               >  1:2849271/1‑148 (MQ=255)
                                                                                                             cTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCgg                                        <  2:1508874/148‑1 (MQ=255)
                                                                                                              tgtgGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGa                                       <  1:1636845/148‑1 (MQ=255)
                                                                                                                tgGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATg                                     <  2:1428260/148‑1 (MQ=255)
                                                                                                                tgGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATg                                     <  2:2875209/148‑1 (MQ=255)
                                                                                                                  gTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGAc                                   <  1:1891975/148‑1 (MQ=255)
                                                                                                                   tCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACa                                  <  2:2187375/148‑1 (MQ=255)
                                                                                                                          cGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAg                           >  2:3728238/1‑148 (MQ=255)
                                                                                                                           gATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGc                          <  2:3059385/148‑1 (MQ=255)
                                                                                                                           gATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGc                          <  1:178981/148‑1 (MQ=255)
                                                                                                                              tAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTg                       >  1:3657622/1‑148 (MQ=255)
                                                                                                                                  gCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGaaaa                   >  2:734718/1‑148 (MQ=255)
                                                                                                                                   cTGCATCGTTGCTTGGCGATCATAGCGAGCACGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAAc                  >  1:3088983/1‑148 (MQ=255)
                                                                                                                                    tGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACg                 >  1:2181619/1‑148 (MQ=255)
                                                                                                                                    tGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACg                 >  1:252402/1‑148 (MQ=255)
                                                                                                                                     gCATCGTTGCTTGGCGATCATAGCGAGCACGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACgg                >  1:2542391/1‑148 (MQ=255)
                                                                                                                                     gCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACgg                >  1:1845537/1‑148 (MQ=255)
                                                                                                                                     gCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACgg                <  1:1047113/148‑1 (MQ=255)
                                                                                                                                     gCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACgg                <  2:2168085/148‑1 (MQ=255)
                                                                                                                                       aTCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTc              <  1:3545842/148‑1 (MQ=255)
                                                                                                                                       aTCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCGGGATGGAAGGGCTGCGCGGCGTCGGTGTACTGTCGCTGGCGGTGTGTACCATGGGGGTGTTTGCTCGCTTTGCACCGGATGACAGCACCAGCCTGAACACGGTc              >  2:1897846/1‑148 (MQ=255)
                                                                                                                                           ttGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTAGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTa          >  2:3388159/1‑148 (MQ=255)
                                                                                                                                            tGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTat         >  2:3616840/1‑148 (MQ=255)
                                                                                                                                             gCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTatc        <  2:2346390/148‑1 (MQ=255)
                                                                                                                                               ttGGAGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTatcat      <  2:178366/148‑1 (MQ=255)
                                                                                                                                                tGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTatcatc     >  1:2739269/1‑148 (MQ=255)
                                                                                                                                                 ggCGATCATAGCGAGCTAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCt    <  2:1019466/148‑1 (MQ=255)
                                                                                                                                                   cGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCtat  <  2:3324999/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTCCGCGATGATCATTACCGGATTACTGGGTCTGTTAATGGCAACGTATCCACCGCTGTGGGTAATGCTCTGTATTCAGATCGCCTTTGCGATAACGACGATTTTAATGCTGTGGTCGGTGTCGATTAAAGCTGCATCGTTGCTTGGCGATCATAGCGAGCAAGGGAAAATTATGGGCTGGATGGAAGGGCTGCGCGGCGTCGGTGTAATGTCGCTGGCGGTGTTTACCATGTGGGTCTTTTCTCGCTTTGCACCGGATGACAGCACCAGCCTGAAAACGGTCATTATCATCTAT  >  CP009273/2894187‑2894481

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: