Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,396,082 G→T 49.8% 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,0820GT49.8% 71.7 / inf 286R204S (CGT→AGT) csrDtargeting factor for csrBC sRNA degradation
Reads supporting (aligned to +/- strand):  ref base G (77/66);  new base T (68/75);  total (145/141)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

CTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTACC  >  CP009273/3395935‑3396230
                                                                                                                                                   |                                                                                                                                                    
cTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTaa                                                                                                                                                    <  2:3084281/150‑1 (MQ=255)
cTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTCGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTaa                                                                                                                                                    <  2:4926038/150‑1 (MQ=255)
cTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTAATTGAGGCCGGTTTTCACGTCCTGGGCGGCCTAAGAGAGGATCAGCGTATCAACAAGGCTGTGTTTTTCGCGTTCGTTCTGAATTTCACaaa                                                                                                                                                    >  2:4487044/1‑150 (MQ=255)
 tACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAag                                                                                                                                                   <  1:1431547/150‑1 (MQ=255)
 tACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAag                                                                                                                                                   <  2:4315292/150‑1 (MQ=255)
  aCTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCCTAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTAACGAagc                                                                                                                                                  >  1:1807699/1‑150 (MQ=255)
  aCTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAagc                                                                                                                                                  <  1:1524909/150‑1 (MQ=255)
  aCTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTAACGTCCTGGGCGGCATAAGAGCGGATCAGCGGATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAagc                                                                                                                                                  >  1:2250537/1‑150 (MQ=255)
        tcCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTa                                                                                                                                            <  1:1032767/150‑1 (MQ=255)
        tcCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTa                                                                                                                                            <  1:3651554/150‑1 (MQ=255)
        tcCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGGGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCgga                                                                                                                                            >  1:3114732/1‑148 (MQ=255)
          cTGATCTTCCAGTAACGTGGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATc                                                                                                                                          <  2:930829/150‑1 (MQ=255)
           tGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATcc                                                                                                                                         <  2:3469391/150‑1 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCa                                                                                                                                        <  1:489069/150‑1 (MQ=255)
             aTCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAg                                                                                                                                       >  2:4375990/1‑150 (MQ=255)
               cTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAgcg                                                                                                                                     <  2:2417748/150‑1 (MQ=255)
                  ccAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCAc                                                                                                                                  >  1:3487765/1‑150 (MQ=255)
                   cAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCAct                                                                                                                                 <  2:255323/150‑1 (MQ=255)
                      tAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgct                                                                                                                              <  1:3018919/150‑1 (MQ=255)
                      tAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCActgct                                                                                                                              >  2:5209961/1‑150 (MQ=255)
                          gTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGtt                                                                                                                          >  2:3593308/1‑150 (MQ=255)
                           ttGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAATAGCGGATCAGCGTATAAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTc                                                                                                                         >  2:4008129/1‑150 (MQ=255)
                           ttGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTc                                                                                                                         <  1:1942322/150‑1 (MQ=255)
                           ttGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTc                                                                                                                         <  2:3015076/150‑1 (MQ=255)
                            tGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCt                                                                                                                        >  1:1236061/1‑150 (MQ=255)
                             gCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTg                                                                                                                       >  2:1526597/1‑150 (MQ=255)
                               tAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTggg                                                                                                                     >  2:2490895/1‑150 (MQ=255)
                                aaCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggc                                                                                                                    <  1:1672248/150‑1 (MQ=255)
                                 aCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggcg                                                                                                                   >  1:4812254/1‑150 (MQ=255)
                                 aCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggcg                                                                                                                   <  1:4571249/150‑1 (MQ=255)
                                 aCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggcg                                                                                                                   >  2:1063980/1‑150 (MQ=255)
                                    gATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGcc                                                                                                                <  1:4952758/150‑1 (MQ=255)
                                    gATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCCGGTTCTGGGCGGcc                                                                                                                >  1:3733354/1‑150 (MQ=255)
                                     aTTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCa                                                                                                               >  1:5144627/1‑150 (MQ=255)
                                              aaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGAt                                                                                                      <  1:4121660/150‑1 (MQ=255)
                                                  agTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGtt                                                                                                  <  2:604725/150‑1 (MQ=255)
                                                  agTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGtt                                                                                                  >  1:671847/1‑150 (MQ=255)
                                                   gTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTc                                                                                                 <  1:3855862/150‑1 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTACCa                                                                                              >  2:4654699/1‑150 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCa                                                                                              <  2:5039220/150‑1 (MQ=255)
                                                       gTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCaa                                                                                             >  2:2707789/1‑150 (MQ=255)
                                                         tATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAaca                                                                                           >  1:3976591/1‑150 (MQ=255)
                                                          aTTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAacac                                                                                          <  1:4017094/150‑1 (MQ=255)
                                                           ttGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAAcaca                                                                                         >  1:426186/1‑150 (MQ=255)
                                                           ttGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGGATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAAcaca                                                                                         <  2:4499049/150‑1 (MQ=255)
                                                           tagAGGCCGGTCTGCACGTCCCGGGCGGCATTAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAAcaca                                                                                         <  2:3241801/148‑1 (MQ=255)
                                                            tGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACAt                                                                                        >  1:3927350/1‑150 (MQ=255)
                                                               ggCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTag                                                                                     >  1:750545/1‑150 (MQ=255)
                                                               ggCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTag                                                                                     <  2:4247/150‑1 (MQ=255)
                                                                     ttttCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCAc                                                                               <  2:865905/150‑1 (MQ=255)
                                                                        tCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCt                                                                            <  2:529077/150‑1 (MQ=255)
                                                                         cACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTc                                                                           <  2:1399184/150‑1 (MQ=255)
                                                                         cACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTc                                                                           >  1:3786331/1‑150 (MQ=255)
                                                                           cGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGATCCACGCTctc                                                                         >  2:4298203/1‑148 (MQ=255)
                                                                           cGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCAc                                                                         >  2:2066449/1‑150 (MQ=255)
                                                                           cGTCCTGGGCGGCATAAGAGCGGATCAGCGACACAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCAc                                                                         >  2:2610657/1‑150 (MQ=255)
                                                                            gTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCAcc                                                                        <  2:4951976/150‑1 (MQ=255)
                                                                                    cggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGa                                                                >  2:724832/1‑150 (MQ=255)
                                                                                     ggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGAt                                                               <  2:4552945/150‑1 (MQ=255)
                                                                                      gcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATa                                                              <  2:286159/150‑1 (MQ=255)
                                                                                      gcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATa                                                              <  2:4847195/150‑1 (MQ=255)
                                                                                        aTAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACg                                                            >  2:4137447/1‑150 (MQ=255)
                                                                                        aTAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACg                                                            >  2:4433345/1‑150 (MQ=255)
                                                                                           agagCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGt                                                         <  1:940064/150‑1 (MQ=255)
                                                                                            gagCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTa                                                        <  1:572152/150‑1 (MQ=255)
                                                                                            gagCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTa                                                        <  1:777620/150‑1 (MQ=255)
                                                                                             agCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAg                                                       >  1:3308376/1‑150 (MQ=255)
                                                                                                ggATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGccc                                                    >  1:2694068/1‑150 (MQ=255)
                                                                                                    cAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGt                                                >  1:2914096/1‑150 (MQ=255)
                                                                                                     aGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGtt                                               <  1:2137146/150‑1 (MQ=255)
                                                                                                     aGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGtt                                               <  2:5176503/150‑1 (MQ=255)
                                                                                                      gCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGttt                                              <  1:2779970/150‑1 (MQ=255)
                                                                                                           tCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGc                                         >  2:1189932/1‑150 (MQ=255)
                                                                                                            cAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCa                                        <  2:1946815/150‑1 (MQ=255)
                                                                                                             aaGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCaa                                       >  1:2451822/1‑150 (MQ=255)
                                                                                                              aGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAAt                                      <  1:4487031/150‑1 (MQ=255)
                                                                                                                  ggCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCtt                                  >  2:2198823/1‑150 (MQ=255)
                                                                                                                  ggCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCtt                                  >  2:1636515/1‑150 (MQ=255)
                                                                                                                    cTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGc                                <  2:1820923/150‑1 (MQ=255)
                                                                                                                       tgttgtTCGCGTGCGTTCTGAATTTCACTATGCAGCGTTTACTGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACCTTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCg                             <  2:2443504/150‑1 (MQ=255)
                                                                                                                       tgttgtTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCg                             <  2:2981914/150‑1 (MQ=255)
                                                                                                                       tgttgtTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCg                             <  1:4876603/150‑1 (MQ=255)
                                                                                                                        gttgttCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCgg                            >  1:4776930/1‑150 (MQ=255)
                                                                                                                                 gtgcgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGcc                   >  1:2538642/1‑150 (MQ=255)
                                                                                                                                 gtgcgtTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTAAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGcc                   >  2:1604927/1‑150 (MQ=255)
                                                                                                                                  tgcgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCg                  <  1:1703814/150‑1 (MQ=255)
                                                                                                                                    cgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGtt                <  1:3145341/150‑1 (MQ=255)
                                                                                                                                    cgtTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGtt                <  1:2490895/150‑1 (MQ=255)
                                                                                                                                     gtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTg               >  1:1381828/1‑150 (MQ=255)
                                                                                                                                     gtTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTg               >  1:3333331/1‑150 (MQ=255)
                                                                                                                                     gtTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTg               <  1:5123430/150‑1 (MQ=255)
                                                                                                                                        cTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTaa            >  2:2162235/1‑150 (MQ=255)
                                                                                                                                         tGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAAc           <  2:2939306/150‑1 (MQ=255)
                                                                                                                                         tGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAAc           <  2:1343737/150‑1 (MQ=255)
                                                                                                                                         tGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAAc           >  2:3355350/1‑150 (MQ=255)
                                                                                                                                          gAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAAcc          <  2:3846021/150‑1 (MQ=255)
                                                                                                                                            aTTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAg        <  1:2381718/150‑1 (MQ=255)
                                                                                                                                              ttCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCg      >  1:694599/1‑150 (MQ=255)
                                                                                                                                              ttCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCg      <  1:3371576/150‑1 (MQ=255)
                                                                                                                                                  cGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTAcc  >  2:5132397/1‑150 (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: