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,0820GT95.6% 838.3 / 23.7 262R204S (CGT→AGT) csrDtargeting factor for csrBC sRNA degradation
Reads supporting (aligned to +/- strand):  ref base G (5/7);  new base T (142/107);  total (148/114)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 262 total aligned reads displayed.

CTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGCGTA  >  CP009273/3395935‑3396228
                                                                                                                                                   |                                                                                                                                                  
cTACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGATCGGATCAGCGTATCAAGACGGCTGTGGTGTTCGCGTGCGTTCTGAATTTCACTaa                                                                                                                                                  >  1:4053330/1‑150 (MQ=255)
 tACTTTCTCCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAag                                                                                                                                                 >  1:88058/1‑150 (MQ=255)
        tcCTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTa                                                                                                                                          >  1:3803598/1‑150 (MQ=255)
         ccTGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTAt                                                                                                                                         <  2:1226621/150‑1 (MQ=255)
           tGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATcc                                                                                                                                       >  1:4055232/1‑150 (MQ=255)
           tGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATcc                                                                                                                                       <  1:1958312/150‑1 (MQ=255)
           tGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATcc                                                                                                                                       >  1:1753223/1‑150 (MQ=255)
           tGATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGGTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATcc                                                                                                                                       >  1:2473877/1‑150 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCa                                                                                                                                      <  2:906599/150‑1 (MQ=255)
            gATCTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCa                                                                                                                                      <  1:3353947/150‑1 (MQ=255)
               cTTCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAgcg                                                                                                                                   >  1:2844382/1‑150 (MQ=255)
                ttCCAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAgcgc                                                                                                                                  <  1:3070520/150‑1 (MQ=255)
                   cAGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCAct                                                                                                                               >  2:3179605/1‑150 (MQ=255)
                    aGTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCActg                                                                                                                              <  2:3861204/150‑1 (MQ=255)
                     gTAACGTTGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACTTCCATGGCGGCATAAGAGCAGAACAGCGATTCAAGACGGCACTGTAGATCGCGTCCGTTAAGAAATTAACTAACCAGCCTATCAAGCGCAcagc                                                                                                                             >  2:4080462/1‑150 (MQ=255)
                            tGCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCt                                                                                                                      <  1:4560671/150‑1 (MQ=255)
                             gCTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTg                                                                                                                     <  2:3052744/150‑1 (MQ=255)
                              cTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTgg                                                                                                                    >  1:1924013/1‑150 (MQ=255)
                              cTAACTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTAGCGTGCGTTCTGAAAGTCACTAAGCAGCGTATACAGCGCACTGCTGGTTCTgg                                                                                                                    >  2:4219880/1‑150 (MQ=255)
                              cTAACTGATTATCGAAAAAGAGTCGGGTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCCGCGTATCCAGCGCACTGCTGGTTCTgg                                                                                                                    >  2:4219800/1‑150 (MQ=255)
                                aaCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggc                                                                                                                  <  2:132035/150‑1 (MQ=255)
                                 aCTGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggcg                                                                                                                 <  1:1167935/150‑1 (MQ=255)
                                   tGATTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGggcggc                                                                                                               >  2:1720479/1‑150 (MQ=255)
                                     aTTATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCa                                                                                                             >  2:3521268/1‑150 (MQ=255)
                                      ttATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAt                                                                                                            <  1:2435260/150‑1 (MQ=255)
                                      ttATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAt                                                                                                            <  2:2791025/150‑1 (MQ=255)
                                       tATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAATAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAtt                                                                                                           >  1:4442524/1‑150 (MQ=255)
                                       tATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAtt                                                                                                           <  1:2167549/150‑1 (MQ=255)
                                       tATCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCAtt                                                                                                           <  2:3658520/150‑1 (MQ=255)
                                        aTCGAAAAAGAGTCGGTTATTGATGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCATCGTATCAAGAAGGCTGTTTTTTTAGCGTGCGTTCTGAATTTCACTAAGCAGCGTATACAGCGCACTGCTGGTTCTGGGCGGCCATTc                                                                                                          >  2:3511532/1‑150 (MQ=255)
                                        aTCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTc                                                                                                          >  1:2560546/1‑150 (MQ=255)
                                        aTCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTc                                                                                                          >  1:583445/1‑150 (MQ=255)
                                         tCGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCa                                                                                                         <  1:1819501/150‑1 (MQ=255)
                                          cGAAAAAGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCAt                                                                                                        >  2:605944/1‑150 (MQ=255)
                                            aaaaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAg                                                                                                      >  1:3188439/1‑150 (MQ=255)
                                             aaaaGAGTCGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGa                                                                                                     >  2:1575970/1‑150 (MQ=255)
                                                     cGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTccc                                                                                             >  2:490069/1‑150 (MQ=255)
                                                     cGGTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTccc                                                                                             <  2:4165855/150‑1 (MQ=255)
                                                      ggTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCa                                                                                            <  2:1484147/150‑1 (MQ=255)
                                                       gTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCaa                                                                                           >  1:826025/1‑150 (MQ=255)
                                                       gTTATTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCaa                                                                                           <  1:1586541/150‑1 (MQ=255)
                                                          aTTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAacac                                                                                        <  2:3979361/150‑1 (MQ=255)
                                                          aTTGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAacac                                                                                        >  1:2364334/1‑150 (MQ=255)
                                                           ttGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAAcaca                                                                                       >  2:1800454/1‑150 (MQ=255)
                                                            tGAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACAt                                                                                      >  1:2968572/1‑150 (MQ=255)
                                                             gAGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCTGCCATTCATAGATGGTTCCCAACACAtt                                                                                     >  2:816481/1‑150 (MQ=255)
                                                              aGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTa                                                                                    <  1:1403401/150‑1 (MQ=255)
                                                              aGGCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTa                                                                                    >  2:1931801/1‑150 (MQ=255)
                                                               ggCCGGTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTag                                                                                   >  1:1536679/1‑150 (MQ=255)
                                                                    gTTTTCACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCa                                                                              <  1:4596925/150‑1 (MQ=255)
                                                                         cACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTc                                                                         >  1:1528210/1‑150 (MQ=255)
                                                                         cACGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTc                                                                         >  1:3238353/1‑150 (MQ=255)
                                                                          aCGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCa                                                                        <  1:3015179/150‑1 (MQ=255)
                                                                          aCGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCa                                                                        <  1:1001544/150‑1 (MQ=255)
                                                                          aCGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCa                                                                        <  2:2119941/150‑1 (MQ=255)
                                                                          aCGTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCa                                                                        >  2:1374293/1‑150 (MQ=255)
                                                                            gTCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCAcc                                                                      <  1:116642/150‑1 (MQ=255)
                                                                             tCCTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCg                                                                     <  1:3116456/150‑1 (MQ=255)
                                                                              ccTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGt                                                                    >  2:2484155/1‑150 (MQ=255)
                                                                               cTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGtt                                                                   <  2:1580466/150‑1 (MQ=255)
                                                                               cTGGGCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGtt                                                                   >  2:741125/1‑150 (MQ=255)
                                                                                 gggCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTa                                                                 >  1:3322508/1‑150 (MQ=255)
                                                                                 gggCGGCATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTa                                                                 >  1:805399/1‑150 (MQ=255)
                                                                                  ggcggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTaa                                                                <  1:3930017/150‑1 (MQ=255)
                                                                                  ggcggcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTaa                                                                >  2:2349111/1‑150 (MQ=255)
                                                                                      gcATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATa                                                            <  1:3531526/150‑1 (MQ=255)
                                                                                       cATAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATAc                                                           >  1:136084/1‑150 (MQ=255)
                                                                                         tAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGa                                                         <  1:4487426/150‑1 (MQ=255)
                                                                                         tAAGAGCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGa                                                         <  2:1030545/150‑1 (MQ=255)
                                                                                           agagCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGt                                                       >  2:4375170/1‑150 (MQ=255)
                                                                                            gagCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTa                                                      <  2:4421922/150‑1 (MQ=255)
                                                                                            gagCGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTa                                                      >  2:3251063/1‑150 (MQ=255)
                                                                                               cGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGcc                                                   >  1:62971/1‑150 (MQ=255)
                                                                                               cGGATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGcc                                                   <  2:3745608/150‑1 (MQ=255)
                                                                                                ggATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGccc                                                  >  2:4070995/1‑150 (MQ=255)
                                                                                                ggATCAGCGTATCAAGACGGCTGTGTTGGTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGccc                                                  >  2:1667080/1‑150 (MQ=255)
                                                                                                 gATCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCg                                                 <  1:3304255/150‑1 (MQ=255)
                                                                                                   tCAGCGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCggg                                               >  1:610235/1‑150 (MQ=255)
                                                                                                      gCGTATCAAGACTGCTGTGTTGTTAGCGTGCGTTCTGAATTTCACTAAGCATCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAAAACATTAGAGCCACGCTCACCGATTAAGATACGAGTAGCCCGGGttt                                            <  2:2196533/150‑1 (MQ=255)
                                                                                                       cGTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTc                                           <  2:1165584/150‑1 (MQ=255)
                                                                                                        gTATCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTcc                                          <  1:2619951/150‑1 (MQ=255)
                                                                                                          aTCAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAg                                        <  2:3591346/150‑1 (MQ=255)
                                                                                                            cAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCa                                      >  2:2438071/1‑150 (MQ=255)
                                                                                                            cAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCa                                      >  2:127210/1‑150 (MQ=255)
                                                                                                            cAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCa                                      >  1:3556816/1‑150 (MQ=255)
                                                                                                            cAAGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCa                                      >  2:1678907/1‑150 (MQ=255)
                                                                                                             aaGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGATTAAGATACGATTAGCCCGGGTTTCCAGCaa                                     >  2:3347699/1‑150 (MQ=255)
                                                                                                              aGACGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAAt                                    >  2:4600176/1‑150 (MQ=255)
                                                                                                                 cGGCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCt                                 >  2:1880102/1‑150 (MQ=255)
                                                                                                                   gCTGTGTTGTTCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTg                               <  1:1971696/150‑1 (MQ=255)
                                                                                                                          tgttCGCGTGCGTTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCa                        >  1:3986422/1‑150 (MQ=255)
                                                                                                                               gcgtgcgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTg                   <  1:2726810/150‑1 (MQ=255)
                                                                                                                                cgtgcgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGc                  >  1:3108255/1‑150 (MQ=255)
                                                                                                                                cgtgcgtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGc                  >  2:1574875/1‑150 (MQ=255)
                                                                                                                                     gtTCTGAATTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTg             >  2:4521995/1‑150 (MQ=255)
                                                                                                                                         tGAATTTCACGAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAAc         <  1:4303526/150‑1 (MQ=255)
                                                                                                                                           aaTTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCa       >  1:1284479/1‑150 (MQ=255)
                                                                                                                                           aaTTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCa       <  1:338466/150‑1 (MQ=255)
                                                                                                                                           aaTTTCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCa       <  1:1049595/150‑1 (MQ=255)
                                                                                                                                             tttCACTAAGCAGCGTATCCAGCGCACTGCTGGTTCTGGGCGGCCATTCATAGATGGTTCCCAACACATTAGAGCCACGCTCACCGTTTAAGATACGAGTAGCCCGGGTTTCCAGCAATTCTTGCCCGGCAAGTTGCCGTTGTAACCAGc     >  2:2915015/1‑150 (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: ‑