Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,396,082 G→T 100% 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,0820GT99.3% 918.3 / ‑2.3 282R204S (CGT→AGT) csrDtargeting factor for csrBC sRNA degradation
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (144/134);  total (146/136)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 282 total aligned reads displayed.

CTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTA  >  CP009273/3395935‑3396228
                                                                                                                                                   |                                                                                                                                                  
ctactaTCTCCTTATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTTAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATAAACACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTaa                                                                                                                                                  >  2:1576411/1‑150 (MQ=255)
  aCTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAagc                                                                                                                                                <  1:180931/150‑1 (MQ=255)
   cTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAagca                                                                                                                                               <  2:1019286/150‑1 (MQ=255)
   cTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAagca                                                                                                                                               <  2:3671359/150‑1 (MQ=255)
    tttCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAagcag                                                                                                                                              >  1:3985039/1‑150 (MQ=255)
       ctcCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGGGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGt                                                                                                                                           >  1:4992056/1‑150 (MQ=255)
        tcCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTa                                                                                                                                          <  1:829329/150‑1 (MQ=255)
        tcCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTa                                                                                                                                          <  1:1959793/150‑1 (MQ=255)
         ccTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTAt                                                                                                                                         <  2:1265185/150‑1 (MQ=255)
         ccTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGGGCTGAATTTCACTAAGCAGCGTAt                                                                                                                                         >  1:2769949/1‑150 (MQ=255)
           tGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAAATTCACTAAGCAGCGTATcc                                                                                                                                       >  1:2431893/1‑150 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCa                                                                                                                                      <  2:571063/150‑1 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACAGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCa                                                                                                                                      <  1:318405/150‑1 (MQ=255)
               cTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAgcg                                                                                                                                   <  1:4680494/150‑1 (MQ=255)
                ttCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAgcgc                                                                                                                                  >  1:4045575/1‑150 (MQ=255)
                 tCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCa                                                                                                                                 >  2:2730492/1‑150 (MQ=255)
                   cAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGctct                                                                                                                               <  2:2822058/150‑1 (MQ=255)
                    aGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActg                                                                                                                              >  1:1560754/1‑150 (MQ=255)
                      tAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgct                                                                                                                            <  2:4323352/150‑1 (MQ=255)
                       aaCGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActgctg                                                                                                                           <  1:3678881/150‑1 (MQ=255)
                          gTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGtt                                                                                                                        <  1:539727/150‑1 (MQ=255)
                            tGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCt                                                                                                                      <  1:3778168/150‑1 (MQ=255)
                            tGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCt                                                                                                                      <  2:2297202/150‑1 (MQ=255)
                            tGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCt                                                                                                                      >  1:2543109/1‑150 (MQ=255)
                            tGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCt                                                                                                                      <  1:2190945/150‑1 (MQ=255)
                            tGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCt                                                                                                                      <  1:931651/150‑1 (MQ=255)
                             gCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTAACAAAGCAGCGTATCCAGCGCACTGCTGGTTCTg                                                                                                                     >  1:3026486/1‑150 (MQ=255)
                               tAACTGATTATCGAAAAAGAGTCGGGTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCCTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATATCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTggg                                                                                                                   >  2:5038182/1‑150 (MQ=255)
                                aaCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggc                                                                                                                  >  1:5044477/1‑150 (MQ=255)
                                 aCTGATTATCGAAAAAGAGTCGGTTATTTTGGCCGGTTTTAACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTAGCGTGCGGTCTGAATTTCAATAAGCAGCGTATCCAGCGCACAGCAGGTTCTGggcg                                                                                                                 >  1:1042863/1‑150 (MQ=255)
                                 aCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggcg                                                                                                                 <  1:1254251/150‑1 (MQ=255)
                                  cTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggcgg                                                                                                                >  2:1376263/1‑150 (MQ=255)
                                   tGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCAATGCTGGTTATGggcggc                                                                                                               >  1:2918002/1‑150 (MQ=255)
                                    gATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGcc                                                                                                              >  1:586944/1‑150 (MQ=255)
                                      ttATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAt                                                                                                            >  2:3281480/1‑150 (MQ=255)
                                      ttATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAt                                                                                                            >  2:4705522/1‑150 (MQ=255)
                                       tATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAtt                                                                                                           >  2:3568846/1‑150 (MQ=255)
                                          cGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCAt                                                                                                        >  2:2329256/1‑150 (MQ=255)
                                           gAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATa                                                                                                       <  2:2361001/150‑1 (MQ=255)
                                           gAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATa                                                                                                       >  2:789452/1‑150 (MQ=255)
                                            aaaaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAg                                                                                                      >  2:1387296/1‑150 (MQ=255)
                                              aaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAATAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATACATCGCACTGCTGGTTCTGGGCGGCCATTCATAGAt                                                                                                    >  1:344888/1‑150 (MQ=255)
                                              aaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGAt                                                                                                    >  2:4483297/1‑150 (MQ=255)
                                               aaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATg                                                                                                   >  1:4426536/1‑150 (MQ=255)
                                                  agTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGtt                                                                                                <  1:4221236/150‑1 (MQ=255)
                                                  agTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGtt                                                                                                >  2:4273787/1‑150 (MQ=255)
                                                    tCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTcc                                                                                              >  2:3830923/1‑150 (MQ=255)
                                                    tCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTcc                                                                                              >  1:4304621/1‑150 (MQ=255)
                                                     cGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTccc                                                                                             >  1:1105546/1‑150 (MQ=255)
                                                        ttATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAac                                                                                          >  1:1873312/1‑150 (MQ=255)
                                                           ttGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAAcaca                                                                                       >  1:1509652/1‑150 (MQ=255)
                                                              aGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTa                                                                                    >  2:475383/1‑150 (MQ=255)
                                                              aGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTa                                                                                    <  2:2012238/150‑1 (MQ=255)
                                                                   ggTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGcc                                                                               <  1:516969/150‑1 (MQ=255)
                                                                   ggTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGcc                                                                               >  2:1933991/1‑150 (MQ=255)
                                                                    gTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCATATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCa                                                                              <  1:1046275/150‑1 (MQ=255)
                                                                     ttttCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCAc                                                                             <  1:3493458/150‑1 (MQ=255)
                                                                      tttCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACg                                                                            <  1:3896750/150‑1 (MQ=255)
                                                                       ttCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGc                                                                           >  2:3934206/1‑150 (MQ=255)
                                                                        tCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCt                                                                          <  1:3842633/150‑1 (MQ=255)
                                                                        tCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCt                                                                          <  2:3144090/150‑1 (MQ=255)
                                                                         cACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTc                                                                         >  1:4824533/1‑150 (MQ=255)
                                                                           cGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCAc                                                                       <  1:4096610/150‑1 (MQ=255)
                                                                            gTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCAcc                                                                      >  1:1706409/1‑150 (MQ=255)
                                                                            gTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCAcc                                                                      <  2:2428959/150‑1 (MQ=255)
                                                                                tGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGttt                                                                  <  2:4423173/150‑1 (MQ=255)
                                                                                tGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGttt                                                                  >  1:4950203/1‑150 (MQ=255)
                                                                                  ggcggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTaa                                                                <  1:4492690/150‑1 (MQ=255)
                                                                                  ggcggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTaa                                                                <  2:3304628/150‑1 (MQ=255)
                                                                                     ggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGAt                                                             <  2:2627302/150‑1 (MQ=255)
                                                                                          acgcgCGGATCAGCGTATCAAGACGGCTGTGTTGGTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAg                                                        <  2:4399105/146‑1 (MQ=255)
                                                                                          aaGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAg                                                        <  1:553889/150‑1 (MQ=255)
                                                                                                 gATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCg                                                 <  1:3569199/150‑1 (MQ=255)
                                                                                                    cAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGt                                              <  2:4045575/150‑1 (MQ=255)
                                                                                                      gCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGttt                                            <  2:178779/150‑1 (MQ=255)
                                                                                                       cGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTc                                           >  2:3406033/1‑150 (MQ=255)
                                                                                                       cGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTc                                           <  2:1529181/150‑1 (MQ=255)
                                                                                                        gTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTcc                                          <  2:3689232/150‑1 (MQ=255)
                                                                                                         tATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCa                                         <  1:2376603/150‑1 (MQ=255)
                                                                                                         tATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCa                                         >  1:3340016/1‑150 (MQ=255)
                                                                                                          aTCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAg                                        >  2:2015325/1‑150 (MQ=255)
                                                                                                           tCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGc                                       <  1:4636942/150‑1 (MQ=255)
                                                                                                              aGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAAt                                    >  1:4355125/1‑150 (MQ=255)
                                                                                                              aGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAAt                                    >  2:3097460/1‑150 (MQ=255)
                                                                                                               gACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAAtt                                   >  1:3643201/1‑150 (MQ=255)
                                                                                                                  ggCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCtt                                <  1:3288607/150‑1 (MQ=255)
                                                                                                                     tgtgTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTTTCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCTATTCTTGcc                             <  1:2609054/150‑1 (MQ=255)
                                                                                                                     tgtgTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGcc                             >  2:4572198/1‑150 (MQ=255)
                                                                                                                      gtgtTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGccc                            <  1:1454212/150‑1 (MQ=255)
                                                                                                                              cgcgTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGtt                    >  2:617316/1‑150 (MQ=255)
                                                                                                                              cgcgTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGtt                    <  2:4977627/150‑1 (MQ=255)
                                                                                                                               gcgtgcgtTCTGAATTTCACTCAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTg                   <  2:3453288/150‑1 (MQ=255)
                                                                                                                                cgtgcgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGc                  <  1:1412187/150‑1 (MQ=255)
                                                                                                                                     gtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTg             <  1:4518482/150‑1 (MQ=255)
                                                                                                                                         tGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCACGGCAAGTTGCCGTTGTAAc         >  2:2424175/1‑150 (MQ=255)
                                                                                                                                           aaTTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCa       >  1:2994526/1‑150 (MQ=255)
                                                                                                                                           aaTTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCa       <  2:3272414/150‑1 (MQ=255)
                                                                                                                                           aaTTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCa       <  2:2543109/150‑1 (MQ=255)
                                                                                                                                             tttCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGc     <  2:2872950/150‑1 (MQ=255)
                                                                                                                                              ttCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCg    <  2:979410/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTA  >  CP009273/3395935‑3396228

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: