Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,269,251 A→T 18.1% Q282L (CAG→CTG)  yjcD → inner membrane putative guanine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,269,2510AT18.1% 456.2 / 91.0 192Q282L (CAG→CTG) yjcDinner membrane putative guanine permease
Reads supporting (aligned to +/- strand):  ref base A (68/89);  new base T (21/14);  total (89/103)
Fisher's exact test for biased strand distribution p-value = 9.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 194 total aligned reads displayed.

CTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCG  >  CP009273/4269104‑4269397
                                                                                                                                                   |                                                                                                                                                  
cTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCa                                                                                                                                                    <  1:823656/148‑1 (MQ=255)
 ttCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAg                                                                                                                                                   >  1:2620310/1‑148 (MQ=255)
      tGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGaa                                                                                                                                              <  1:354132/148‑1 (MQ=255)
      tGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGaa                                                                                                                                              >  1:654097/1‑148 (MQ=255)
       gAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAAc                                                                                                                                             <  2:530292/148‑1 (MQ=255)
            aTGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACctgct                                                                                                                                        <  1:1609553/148‑1 (MQ=255)
              gAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTgg                                                                                                                                      >  1:2225709/1‑148 (MQ=255)
              gAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTgg                                                                                                                                      >  2:3766297/1‑148 (MQ=255)
              gAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTgg                                                                                                                                      <  2:3585969/148‑1 (MQ=255)
              gAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTgg                                                                                                                                      >  2:2247028/1‑148 (MQ=255)
                aaaCGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGAt                                                                                                                                    >  2:357379/1‑148 (MQ=255)
                  aCGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATaa                                                                                                                                  >  1:713810/1‑148 (MQ=255)
                      cAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGAt                                                                                                                              <  2:3671174/148‑1 (MQ=255)
                       aaTTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATg                                                                                                                             <  2:2378506/148‑1 (MQ=255)
                        aTTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATgg                                                                                                                            >  1:3083120/1‑148 (MQ=255)
                         ttCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGACCCTGCTGGATAAAGATggg                                                                                                                           >  2:2296324/1‑148 (MQ=255)
                             cTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGa                                                                                                                       >  1:1757484/1‑148 (MQ=255)
                              tGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGat                                                                                                                      <  1:335598/148‑1 (MQ=255)
                               gATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGatc                                                                                                                     >  1:2983433/1‑148 (MQ=255)
                               gATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGatc                                                                                                                     >  1:3114824/1‑148 (MQ=255)
                               gATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGatc                                                                                                                     <  2:2034511/148‑1 (MQ=255)
                                  tGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGatcatc                                                                                                                  >  2:928919/1‑148 (MQ=255)
                                   ggCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGTTGGATAAAGATGGGCAGATCATCg                                                                                                                 <  2:924005/148‑1 (MQ=255)
                                    gCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGa                                                                                                                <  2:1927292/148‑1 (MQ=255)
                                         cTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTTTTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTg                                                                                                           <  1:200776/148‑1 (MQ=255)
                                           ggACATTATGGGCGCGCTTAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGTATAAAGATGGGCAGATCATCGACGGTggg                                                                                                         <  1:3013660/148‑1 (MQ=255)
                                           ggACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTggg                                                                                                         <  2:1356745/148‑1 (MQ=255)
                                            gACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGa                                                                                                        >  1:912727/1‑148 (MQ=255)
                                             aCATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGaa                                                                                                       <  2:1125632/148‑1 (MQ=255)
                                             aCATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGaa                                                                                                       <  2:811040/148‑1 (MQ=255)
                                                 tATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGca                                                                                                   >  2:119186/1‑148 (MQ=255)
                                                 tATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGca                                                                                                   >  2:3751251/1‑148 (MQ=255)
                                                 tATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGca                                                                                                   >  1:1397061/1‑148 (MQ=255)
                                                 tATGGGCGCGCGGAATCCTGTAGGCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGca                                                                                                   <  2:968266/148‑1 (MQ=255)
                                                  aTGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcac                                                                                                  <  2:2176572/148‑1 (MQ=255)
                                                    gggCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTGTGGCGCTGTTGATGACTGCTTTATTTGATGCAACCGGAACTATCCGTGCCGTGGCCGGCCTGGCGAACCTGCTGGATAAAGAGGGGCTGATCATCGACGGTGGGAAAGcactg                                                                                                <  1:2867730/148‑1 (MQ=255)
                                                        gcgcTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcactgacca                                                                                            <  2:1191097/148‑1 (MQ=255)
                                                         cgcTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcactgaccac                                                                                           <  1:2463631/148‑1 (MQ=255)
                                                            tGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcactgaccactga                                                                                        >  1:255843/1‑148 (MQ=255)
                                                             gAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcactgaccactgac                                                                                       >  1:3487349/1‑148 (MQ=255)
                                                             gAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcactgaccactgac                                                                                       >  2:3463716/1‑148 (MQ=255)
                                                              aaTCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACt                                                                                      >  2:2159996/1‑148 (MQ=255)
                                                              aaTCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACt                                                                                      >  2:3339647/1‑148 (MQ=255)
                                                              aaTCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACt                                                                                      >  2:3339687/1‑148 (MQ=255)
                                                                 ccTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCa                                                                                   >  1:1165665/1‑148 (MQ=255)
                                                                  cTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCAt                                                                                  <  1:596564/148‑1 (MQ=255)
                                                                    gTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGa                                                                                >  1:3634225/1‑148 (MQ=255)
                                                                     tAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGag                                                                               >  2:2375622/1‑148 (MQ=255)
                                                                     tAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGag                                                                               >  2:2376280/1‑148 (MQ=255)
                                                                       gTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGagca                                                                             >  1:3458591/1‑148 (MQ=255)
                                                                       gTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGagca                                                                             >  1:3459105/1‑148 (MQ=255)
                                                                        tCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGagcag                                                                            <  1:2404586/148‑1 (MQ=255)
                                                                           tGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGagcagctt                                                                         <  1:529673/148‑3 (MQ=255)
                                                                            gCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGtt                                                                        <  1:1522230/148‑1 (MQ=255)
                                                                              cAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGtttt                                                                      <  2:1639632/148‑1 (MQ=255)
                                                                              cAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGtttt                                                                      <  1:1870372/148‑1 (MQ=255)
                                                                                aGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTtct                                                                    <  2:654097/148‑1 (MQ=255)
                                                                                aGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTtct                                                                    <  2:2569429/148‑1 (MQ=255)
                                                                                 gCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCTTGCCGTCGCCGGCCAGGCTAACCTGCTGGATAAAGATGGGCATATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTtctc                                                                   <  1:3784935/148‑1 (MQ=255)
                                                                                 gCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTtctc                                                                   <  2:364212/148‑1 (MQ=255)
                                                                                 gCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGTATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTtctc                                                                   <  1:1633236/148‑1 (MQ=255)
                                                                                 gCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTtctc                                                                   <  2:3086896/148‑1 (MQ=255)
                                                                                    ttCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTgg                                                                >  1:2296133/1‑148 (MQ=255)
                                                                                    ttCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTgg                                                                >  2:255154/1‑148 (MQ=255)
                                                                                    ttCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTgg                                                                >  2:254667/1‑148 (MQ=255)
                                                                                     tCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGca                                                                                                   >  1:2112323/1‑112 (MQ=255)
                                                                                     tCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGca                                                                                                   <  2:2112323/112‑1 (MQ=255)
                                                                                     tCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGc                                                               >  1:1040987/1‑148 (MQ=255)
                                                                                     tCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGc                                                               >  1:1041570/1‑148 (MQ=255)
                                                                                     tCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGc                                                               >  2:1443651/1‑148 (MQ=255)
                                                                                             tGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTggtggg                                                       <  2:3439518/148‑1 (MQ=255)
                                                                                              ggTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTggtgggt                                                      <  2:824117/148‑1 (MQ=255)
                                                                                              ggTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTggtgggt                                                      >  2:343560/1‑148 (MQ=255)
                                                                                               gTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTggtgggtg                                                     <  2:2416825/148‑1 (MQ=255)
                                                                                                  atgaCGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCgg                                                  >  2:3055752/1‑148 (MQ=255)
                                                                                                               ttGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGt                                     <  2:1255740/148‑1 (MQ=255)
                                                                                                               ttGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGt                                     <  2:3502611/148‑1 (MQ=255)
                                                                                                                tGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGta                                    <  2:3674729/148‑1 (MQ=255)
                                                                                                                  aTGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGtata                                  <  1:1313792/148‑1 (MQ=255)
                                                                                                                   tGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATAc                                 <  1:2217955/148‑1 (MQ=255)
                                                                                                                     ccACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACAt                               <  2:1848966/148‑1 (MQ=255)
                                                                                                                     ccACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACAt                               >  1:3181742/1‑148 (MQ=255)
                                                                                                                        ccGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGa                            >  2:319912/1‑148 (MQ=255)
                                                                                                                          ggAAATAGCCGTGCCGTGGCCGGACAGGCGATCATGGTGGATAAAGTTGGGGTGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGt                          <  1:3214217/148‑1 (MQ=255)
                                                                                                                             aCTATCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTg                       <  1:1197249/148‑1 (MQ=255)
                                                                                                                             aCTAGCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTg                       <  2:2011795/148‑1 (MQ=255)
                                                                                                                                aTCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTgcgg                    >  1:3718644/1‑148 (MQ=255)
                                                                                                                                 tCCGTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTgcggc                   >  2:3320531/1‑148 (MQ=255)
                                                                                                                                 tCCGTGCCGTCGCAGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTgcggc                   >  1:851205/1‑148 (MQ=255)
                                                                                                                                    gTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCggg                <  1:2194510/148‑1 (MQ=255)
                                                                                                                                    gTGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCggg                <  1:2159996/148‑1 (MQ=255)
                                                                                                                                     tGCCGTCGCCGGCCTGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGt               >  1:985528/1‑148 (MQ=255)
                                                                                                                                     tGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGt               <  1:3504051/148‑1 (MQ=255)
                                                                                                                                       ccGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGtcc             >  2:1521480/1‑146 (MQ=255)
                                                                                                                                       ccGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTGTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAc             >  2:3647759/1‑148 (MQ=255)
                                                                                                                                           cGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAcggcg         >  2:459841/1‑148 (MQ=255)
                                                                                                                                             ccGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAcggcggc       >  1:1577294/1‑148 (MQ=255)
                                                                                                                                             ccGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAcggcggc       >  1:1577257/1‑148 (MQ=255)
                                                                                                                                                gCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCggcggg    <  1:181684/148‑1 (MQ=255)
                                                                                                                                                  cAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCggcgggcg  >  2:3624788/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGACATTATGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGTCGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTGACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCG  >  CP009273/4269104‑4269397

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: