Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,396,082 G→T 37.1% R204S (CGT→AGT)  csrD ← targeting factor for csrBC sRNA degradation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,396,0820GT37.1% 437.9 / inf 338R204S (CGT→AGT) csrDtargeting factor for csrBC sRNA degradation
Reads supporting (aligned to +/- strand):  ref base G (113/98);  new base T (58/67);  total (173/165)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

CTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTACC  >  CP009273/3395935‑3396230
                                                                                                                                                   |                                                                                                                                                    
cTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGa                                                                                                                                                     <  2:2417870/149‑1 (MQ=255)
cTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGa                                                                                                                                                     <  2:3778510/149‑1 (MQ=255)
 tACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGTGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAATACGGGTGTGTTGTTCGCGTGCGTTCTGAATTTAACGaa                                                                                                                                                    >  2:4517348/1‑149 (MQ=255)
 tACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTaa                                                                                                                                                    <  2:3643981/149‑1 (MQ=255)
 tACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATAAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTaa                                                                                                                                                    >  1:2026488/1‑149 (MQ=255)
  aCTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAag                                                                                                                                                   <  1:2001477/149‑1 (MQ=255)
      tctcCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAagcagc                                                                                                                                               <  1:2078511/149‑1 (MQ=255)
         ccTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTa                                                                                                                                            <  2:1466345/149‑1 (MQ=255)
          cTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTAGCGTGCGTTCTGAATTTCACGAAGCAGCGTAt                                                                                                                                           >  1:1067850/1‑149 (MQ=255)
           tGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATc                                                                                                                                          <  1:4848912/149‑1 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATcc                                                                                                                                         <  1:771032/149‑1 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATcc                                                                                                                                         >  1:1820567/1‑149 (MQ=255)
             aTCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCCTAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCa                                                                                                                                        >  1:4141021/1‑149 (MQ=255)
             aTCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCa                                                                                                                                        >  1:4590686/1‑149 (MQ=255)
              tCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAg                                                                                                                                       >  2:132632/1‑149 (MQ=255)
                ttCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATAAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTAAAGAAGCAGCGTATACAgcg                                                                                                                                     >  1:3095690/1‑149 (MQ=255)
                     gTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActg                                                                                                                                >  2:3877038/1‑149 (MQ=255)
                     gTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCActg                                                                                                                                >  2:1993272/1‑149 (MQ=255)
                     gTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCActg                                                                                                                                <  2:1788074/149‑1 (MQ=255)
                      tAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgc                                                                                                                               <  2:2673298/149‑1 (MQ=255)
                      tAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCActgc                                                                                                                               >  2:2892960/1‑149 (MQ=255)
                       aaCGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgct                                                                                                                              <  2:3057038/149‑1 (MQ=255)
                       aaCGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgct                                                                                                                              >  2:1296371/1‑149 (MQ=255)
                       aaCGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgct                                                                                                                              <  2:1889177/149‑1 (MQ=255)
                          gTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGt                                                                                                                           >  2:4995948/1‑149 (MQ=255)
                           ttGCTAACTGATTATCGAAAAAGAGTCGGTTATTTAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACAGCTGGtt                                                                                                                          >  1:790051/1‑149 (MQ=255)
                           ttGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGtt                                                                                                                          <  1:1240504/149‑1 (MQ=255)
                               tAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTgg                                                                                                                      >  2:4257745/1‑149 (MQ=255)
                                aaCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTggg                                                                                                                     <  1:1163328/149‑1 (MQ=255)
                                       tATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAt                                                                                                              >  2:2957652/1‑149 (MQ=255)
                                        aTCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAtt                                                                                                             <  2:1670080/149‑1 (MQ=255)
                                         tCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATAAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATACAGCGCACTGCTGGTTCTGGGCGGCCATTc                                                                                                            >  1:3667334/1‑149 (MQ=255)
                                          cGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCa                                                                                                           <  1:682237/149‑1 (MQ=255)
                                            aaaaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATa                                                                                                         >  2:3569016/1‑149 (MQ=255)
                                             aaaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAg                                                                                                        <  2:3167357/149‑1 (MQ=255)
                                              aaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGa                                                                                                       <  2:2289606/149‑1 (MQ=255)
                                               aaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGAt                                                                                                      <  1:833822/149‑1 (MQ=255)
                                                agatTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATg                                                                                                     >  1:4948795/1‑149 (MQ=255)
                                                agagTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATg                                                                                                     >  2:434751/1‑149 (MQ=255)
                                                   gTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGtt                                                                                                  <  2:4855894/149‑1 (MQ=255)
                                                    tCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTc                                                                                                 >  1:4102625/1‑149 (MQ=255)
                                                     cGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTcc                                                                                                <  1:2688521/149‑1 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTccc                                                                                               <  2:7385/149‑1 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTccc                                                                                               <  1:2193464/149‑1 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTccc                                                                                               >  2:3737372/1‑149 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTccc                                                                                               >  1:5175943/1‑149 (MQ=255)
                                                        ttATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCaa                                                                                             >  2:3915995/1‑149 (MQ=255)
                                                         tATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAac                                                                                            >  1:1746011/1‑149 (MQ=255)
                                                             gAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACAt                                                                                        <  1:2425608/149‑1 (MQ=255)
                                                              aGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACAtt                                                                                       <  1:4008013/149‑1 (MQ=255)
                                                               ggCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCATGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTa                                                                                      >  1:2817005/1‑149 (MQ=255)
                                                               ggCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTa                                                                                      >  1:298033/1‑149 (MQ=255)
                                                                gCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTag                                                                                     >  2:4358873/1‑149 (MQ=255)
                                                                 ccGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTaga                                                                                    >  2:4472705/1‑149 (MQ=255)
                                                                  cGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTagag                                                                                   <  1:3202372/149‑1 (MQ=255)
                                                                      tttCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCAc                                                                               <  1:3547295/149‑1 (MQ=255)
                                                                        tCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGc                                                                             <  1:792628/149‑1 (MQ=255)
                                                                        tCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGc                                                                             <  2:4797118/149‑1 (MQ=255)
                                                                             tCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCAcc                                                                        >  1:2482068/1‑149 (MQ=255)
                                                                                tGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGtt                                                                     <  1:3306027/149‑1 (MQ=255)
                                                                                tGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGtt                                                                     <  2:3669944/149‑1 (MQ=255)
                                                                                 gggCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGttt                                                                    <  1:3400972/149‑1 (MQ=255)
                                                                                 gggCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGttt                                                                    <  1:4837170/149‑1 (MQ=255)
                                                                                     ggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGa                                                                >  1:441727/1‑149 (MQ=255)
                                                                                      gcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGAt                                                               <  2:4614629/149‑1 (MQ=255)
                                                                                      gcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGAt                                                               <  2:2771648/149‑1 (MQ=255)
                                                                                       cATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATa                                                              <  1:1322355/149‑1 (MQ=255)
                                                                                       cATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATa                                                              <  2:324797/149‑1 (MQ=255)
                                                                                         tAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACg                                                            >  1:2790208/1‑149 (MQ=255)
                                                                                             agCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTa                                                        >  1:289567/1‑149 (MQ=255)
                                                                                             agCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTa                                                        <  1:2924099/149‑1 (MQ=255)
                                                                                              gCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAg                                                       <  1:4564508/149‑1 (MQ=255)
                                                                                                  aTCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCg                                                   <  2:3214093/149‑1 (MQ=255)
                                                                                                   tCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCgg                                                  <  2:1549339/149‑1 (MQ=255)
                                                                                                         tATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTcc                                            <  1:4455732/149‑1 (MQ=255)
                                                                                                             aaGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCa                                        >  1:1447706/1‑149 (MQ=255)
                                                                                                               gACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAAt                                      <  2:401383/149‑1 (MQ=255)
                                                                                                               gACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAAt                                      <  2:1051081/149‑1 (MQ=255)
                                                                                                                 cGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTc                                    <  2:595086/149‑1 (MQ=255)
                                                                                                                  ggCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCt                                   >  1:101230/1‑149 (MQ=255)
                                                                                                                   gCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCtt                                  >  2:669971/1‑149 (MQ=255)
                                                                                                                       tgttgtTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGccc                              <  2:2643096/149‑1 (MQ=255)
                                                                                                                       tgttgtTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGccc                              <  2:4733530/149‑1 (MQ=255)
                                                                                                                           gttCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCa                          >  1:944910/1‑149 (MQ=255)
                                                                                                                            ttCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCaa                         <  1:172045/149‑1 (MQ=255)
                                                                                                                                 gtgcgtTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGc                    <  1:529131/149‑1 (MQ=255)
                                                                                                                                 gtgcgtTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGc                    >  1:1669941/1‑149 (MQ=255)
                                                                                                                                    cgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGt                 <  1:1986632/149‑1 (MQ=255)
                                                                                                                                     gtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGtt                <  2:2186139/149‑1 (MQ=255)
                                                                                                                                      ttCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTg               >  1:3717189/1‑149 (MQ=255)
                                                                                                                                          gAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAAc           <  2:3493701/149‑1 (MQ=255)
                                                                                                                                          gAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCTTTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAAc           <  2:2542269/149‑1 (MQ=255)
                                                                                                                                           aaTTTCACGAAGCAGCGTATCTAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAAcc          <  2:2809512/149‑1 (MQ=255)
                                                                                                                                            aTTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCa         >  2:3344754/1‑149 (MQ=255)
                                                                                                                                             tttCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAg        >  1:2265445/1‑149 (MQ=255)
                                                                                                                                              ttCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGc       <  1:669937/149‑1 (MQ=255)
                                                                                                                                                 aCGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTa    <  1:3237369/149‑1 (MQ=255)
                                                                                                                                                 aCGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTa    >  2:5203035/1‑149 (MQ=255)
                                                                                                                                                  cGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTAc   >  2:56480/1‑149 (MQ=255)
                                                                                                                                                   gAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTAcc  >  1:4740339/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTACC  >  CP009273/3395935‑3396230

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: