Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,396,082 G→T 87.4% 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,0820GT87.4% 657.2 / 89.2 247R204S (CGT→AGT) csrDtargeting factor for csrBC sRNA degradation
Reads supporting (aligned to +/- strand):  ref base G (14/17);  new base T (127/87);  total (143/104)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 247 total aligned reads displayed.

CTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTAC  >  CP009273/3395935‑3396229
                                                                                                                                                   |                                                                                                                                                   
 tACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAag                                                                                                                                                  >  1:3288214/1‑150 (MQ=255)
  aCTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAagc                                                                                                                                                 >  1:2238290/1‑150 (MQ=255)
    tttCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAagcag                                                                                                                                               >  1:4901940/1‑150 (MQ=255)
    tttCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAagcag                                                                                                                                               <  2:418997/150‑1 (MQ=255)
      tctcCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCg                                                                                                                                             >  1:3887418/1‑150 (MQ=255)
         ccTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTAt                                                                                                                                          <  2:259227/150‑1 (MQ=255)
         ccTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTAt                                                                                                                                          <  1:2879423/150‑1 (MQ=255)
         ccTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTAt                                                                                                                                          <  1:2553005/150‑1 (MQ=255)
          cTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATc                                                                                                                                         >  1:291807/1‑150 (MQ=255)
          cTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATc                                                                                                                                         >  1:2905610/1‑150 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCa                                                                                                                                       >  1:778865/1‑150 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCa                                                                                                                                       <  2:236646/150‑1 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGGTCTGAATTTCACTAAGCAGCGTATCCa                                                                                                                                       >  1:1111388/1‑150 (MQ=255)
             aTCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTNTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAg                                                                                                                                      <  1:4791836/150‑1 (MQ=255)
              tCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAgc                                                                                                                                     >  1:2974063/1‑150 (MQ=255)
                ttCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTAACTAAGCAGCGTATCCAgcgc                                                                                                                                   >  1:305544/1‑150 (MQ=255)
                  ccAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCAc                                                                                                                                 <  2:2139883/150‑1 (MQ=255)
                    aGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActg                                                                                                                               <  2:1648454/150‑1 (MQ=255)
                    aGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActg                                                                                                                               <  1:1784753/150‑1 (MQ=255)
                      tAACGTTGCTAACTGATGATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgct                                                                                                                             <  2:4721/150‑1 (MQ=255)
                       aacggagCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgctg                                                                                                                            <  1:1657882/144‑1 (MQ=255)
                       aaCGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgctg                                                                                                                            >  2:2888195/1‑150 (MQ=255)
                        aCGTTGCTAACTGATTATCGAAAAAGAGTCGGGTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTgg                                                                                                                           >  2:3517095/1‑150 (MQ=255)
                          gTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGtt                                                                                                                         >  1:2466242/1‑150 (MQ=255)
                          gTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGtt                                                                                                                         >  2:4765918/1‑150 (MQ=255)
                             gCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTg                                                                                                                      <  2:3827573/150‑1 (MQ=255)
                                aaCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGTTCAGCGTATCAAGAAGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAATCAGCGTATCCAGCGCACTGCTGGTTCTGggc                                                                                                                   >  2:4558427/1‑150 (MQ=255)
                                aaCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggc                                                                                                                   <  1:405188/150‑1 (MQ=255)
                                   tGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCAATTCCAGGTCGGCAAAAGAGCGGAACAGCTTATCAACACCGCTGTGTTGTTCGCGTGAGTTATGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggcggc                                                                                                                >  2:3219983/1‑150 (MQ=255)
                                     aTTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAgc                                                                                                                                     <  1:4682430/127‑1 (MQ=255)
                                     aTTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCa                                                                                                              >  2:2675230/1‑150 (MQ=255)
                                      ttATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGTGCGGATCAGCGTATCAAGACGGCTGTGGTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAt                                                                                                             >  2:3632259/1‑150 (MQ=255)
                                       tATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAtt                                                                                                            >  2:2320492/1‑150 (MQ=255)
                                        aTCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTc                                                                                                           <  2:4267551/150‑1 (MQ=255)
                                          cGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCAGAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACAGCTGGTTCTGGGCGGCCATTCAt                                                                                                         >  2:987388/1‑150 (MQ=255)
                                            aaaaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAg                                                                                                       <  1:4156388/150‑1 (MQ=255)
                                            aaaaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAg                                                                                                       <  2:4641990/150‑1 (MQ=255)
                                            aaaaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAg                                                                                                       <  2:4642688/150‑1 (MQ=255)
                                              aaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGAt                                                                                                     >  2:3797966/1‑150 (MQ=255)
                                              aaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGAt                                                                                                     <  1:1011134/150‑1 (MQ=255)
                                              aaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGCTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGAt                                                                                                     <  1:972263/150‑1 (MQ=255)
                                               aaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATg                                                                                                    >  2:769211/1‑150 (MQ=255)
                                                 gagTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGt                                                                                                  >  1:3106286/1‑150 (MQ=255)
                                                  agTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGtt                                                                                                 >  1:2614086/1‑150 (MQ=255)
                                                   gTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTc                                                                                                <  1:691606/150‑1 (MQ=255)
                                                   gTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTc                                                                                                >  2:2274415/1‑150 (MQ=255)
                                                     cGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTccc                                                                                              <  2:3485473/150‑1 (MQ=255)
                                                     cGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTccc                                                                                              <  1:4660681/150‑1 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCa                                                                                             <  2:1882990/150‑1 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCa                                                                                             >  2:1446198/1‑150 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCa                                                                                             >  2:4916804/1‑150 (MQ=255)
                                                        ttATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAac                                                                                           >  2:1005224/1‑150 (MQ=255)
                                                        ttATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAac                                                                                           >  2:218410/1‑150 (MQ=255)
                                                          aTTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAacac                                                                                         >  1:4696290/1‑150 (MQ=255)
                                                            tGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACAt                                                                                       <  1:4765918/150‑1 (MQ=255)
                                                              aGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTa                                                                                     <  1:1944148/150‑1 (MQ=255)
                                                               ggCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTag                                                                                    <  1:2109681/150‑1 (MQ=255)
                                                                 ccGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTagag                                                                                  <  2:2569610/150‑1 (MQ=255)
                                                                  cGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGc                                                                                 <  1:3762305/150‑1 (MQ=255)
                                                                     ttttCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCAc                                                                              <  1:4495832/150‑1 (MQ=255)
                                                                      tttCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACg                                                                             >  2:664601/1‑150 (MQ=255)
                                                                      tttCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACg                                                                             <  1:4410785/150‑1 (MQ=255)
                                                                       ttCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGc                                                                            >  1:4342273/1‑150 (MQ=255)
                                                                          aCGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCa                                                                         <  2:3797860/150‑1 (MQ=255)
                                                                            gTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCAcc                                                                       >  2:1727423/1‑150 (MQ=255)
                                                                               cTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGtt                                                                    >  1:4207410/1‑150 (MQ=255)
                                                                                 gggCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTa                                                                  <  1:2696119/150‑1 (MQ=255)
                                                                                     ggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGAt                                                              >  2:1148978/1‑150 (MQ=255)
                                                                                     ggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGAt                                                              >  2:2741857/1‑150 (MQ=255)
                                                                                      gcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATa                                                             >  2:4364222/1‑150 (MQ=255)
                                                                                      gcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATa                                                             >  1:2215876/1‑150 (MQ=255)
                                                                                        aTAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACg                                                           >  1:3400342/1‑150 (MQ=255)
                                                                                          aaGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAg                                                         >  2:326812/1‑150 (MQ=255)
                                                                                          aaGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAg                                                         <  1:1944590/150‑1 (MQ=255)
                                                                                           agagCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGt                                                        <  2:4134919/150‑1 (MQ=255)
                                                                                            gagCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTa                                                       >  2:5052967/1‑150 (MQ=255)
                                                                                             agCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAg                                                      >  1:1225272/1‑150 (MQ=255)
                                                                                             agCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAg                                                      <  1:3172955/150‑1 (MQ=255)
                                                                                                 gATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCg                                                  >  1:1223870/1‑150 (MQ=255)
                                                                                                 gATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCg                                                  >  2:680984/1‑150 (MQ=255)
                                                                                                    cAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGt                                               >  2:1359543/1‑150 (MQ=255)
                                                                                                        gTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTcc                                           <  1:1045430/150‑1 (MQ=255)
                                                                                                         tATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCa                                          >  1:2072482/1‑150 (MQ=255)
                                                                                                          aTCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAg                                         >  2:5143684/1‑150 (MQ=255)
                                                                                                            cAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCa                                       <  1:2392473/150‑1 (MQ=255)
                                                                                                             aaGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCaa                                      <  1:4513872/150‑1 (MQ=255)
                                                                                                                         ttgttCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGc                          >  1:3434693/1‑150 (MQ=255)
                                                                                                                             tCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGt                      >  2:3331744/1‑150 (MQ=255)
                                                                                                                                 gtgcgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGcc                  <  1:4241391/150‑1 (MQ=255)
                                                                                                                                 gtgcgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGcc                  >  2:1120632/1‑150 (MQ=255)
                                                                                                                                 gtgcgtTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGcc                  >  1:3795050/1‑150 (MQ=255)
                                                                                                                                     gtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTg              >  2:5168150/1‑150 (MQ=255)
                                                                                                                                        cTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTaa           <  1:1760139/150‑1 (MQ=255)
                                                                                                                                           aaTTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCa        <  1:553287/150‑1 (MQ=255)
                                                                                                                                            aTTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAg       <  1:1437965/150‑1 (MQ=255)
                                                                                                                                             tttCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGc      <  1:2553790/150‑1 (MQ=255)
                                                                                                                                               tCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGt    >  2:1773073/1‑150 (MQ=255)
                                                                                                                                                cACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTa   <  2:2816209/150‑1 (MQ=255)
                                                                                                                                                cACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTa   >  1:231039/1‑150 (MQ=255)
                                                                                                                                                 aCGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTAc  <  2:3417054/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTAC  >  CP009273/3395935‑3396229

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: