Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,308,595 C→A 17.6% E277* (GAG→TAG)  phnJ ← carbon‑phosphorus lyase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,308,5950CA17.6% 1003.5 / 157.5 333E277* (GAG→TAG) phnJcarbon‑phosphorus lyase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (143/131);  new base A (24/35);  total (167/166)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 333 total aligned reads displayed.

GTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTC  >  CP009273/4308448‑4308742
                                                                                                                                                   |                                                                                                                                                   
 tCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGct                                                                                                                                                  <  1:386460/150‑1 (MQ=255)
 tCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGct                                                                                                                                                  <  1:387932/150‑1 (MQ=255)
 tCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGct                                                                                                                                                  <  2:4692991/150‑1 (MQ=255)
 tCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGct                                                                                                                                                  <  1:4325850/150‑1 (MQ=255)
  cTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAACAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGctc                                                                                                                                                 <  1:4630173/150‑1 (MQ=255)
   ttCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGctct                                                                                                                                                <  2:2548760/150‑1 (MQ=255)
     cccGGAGCCGGATTCCCCGACACTGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGt                                                                                                                                              <  1:2250913/150‑1 (MQ=255)
     cccGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGt                                                                                                                                              <  2:1789652/150‑1 (MQ=255)
         gAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTggc                                                                                                                                          >  1:4208426/1‑150 (MQ=255)
            ccGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTACAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTggcggc                                                                                                                                       >  2:2614152/1‑150 (MQ=255)
             cGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCa                                                                                                                                      <  2:2552083/150‑1 (MQ=255)
             cGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCa                                                                                                                                      <  1:2045707/150‑1 (MQ=255)
              ggATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCaa                                                                                                                                     <  2:2348442/150‑1 (MQ=255)
                aTTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCaata                                                                                                                                   <  1:3455454/150‑1 (MQ=255)
                  tCCCCGACAATGCCCAGCACGTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCaataat                                                                                                                                 <  1:2454093/150‑1 (MQ=255)
                     ccGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTGTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCgg                                                                                                                              <  2:1290196/150‑1 (MQ=255)
                       gACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTa                                                                                                                            >  2:2828313/1‑150 (MQ=255)
                         cAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATc                                                                                                                          <  2:677969/150‑1 (MQ=255)
                         cAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATc                                                                                                                          >  2:2127026/1‑150 (MQ=255)
                          aaTGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAACGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCg                                                                                                                         <  2:1632709/150‑1 (MQ=255)
                               ccAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCa                                                                                                                    >  2:1691844/1‑150 (MQ=255)
                                 aGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGa                                                                                                                  >  1:487226/1‑150 (MQ=255)
                                   cACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATCATCGGTATCGGAGCAGACa                                                                                                                >  2:929897/1‑150 (MQ=255)
                                   cACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACa                                                                                                                >  1:1695865/1‑150 (MQ=255)
                                    aCTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACaa                                                                                                               >  1:938322/1‑150 (MQ=255)
                                          cccGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGc                                                                                                         >  2:2647298/1‑150 (MQ=255)
                                           ccGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCg                                                                                                        <  2:312240/150‑1 (MQ=255)
                                           ccGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCg                                                                                                        <  2:2793798/150‑1 (MQ=255)
                                           ccGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCg                                                                                                        <  1:1265350/150‑1 (MQ=255)
                                           ccGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCg                                                                                                        <  2:1128246/150‑1 (MQ=255)
                                                  tAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCg                                                                                                 >  2:2341272/1‑150 (MQ=255)
                                                  tAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCg                                                                                                 <  1:1755266/150‑1 (MQ=255)
                                                  tAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCg                                                                                                 >  2:2557569/1‑150 (MQ=255)
                                                   aaaTCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCgg                                                                                                >  2:329593/1‑150 (MQ=255)
                                                   aaaTCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCgg                                                                                                >  1:2918470/1‑150 (MQ=255)
                                                       cAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTc                                                                                            >  1:1191715/1‑150 (MQ=255)
                                                       cAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTc                                                                                            >  1:2820454/1‑150 (MQ=255)
                                                        aaaaGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCa                                                                                           >  1:3071963/1‑150 (MQ=255)
                                                           agagACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATcc                                                                                        >  2:4324118/1‑150 (MQ=255)
                                                             agaCATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAg                                                                                      >  2:3397389/1‑150 (MQ=255)
                                                              gaCATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGc                                                                                     <  1:1796956/150‑1 (MQ=255)
                                                              gaCATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGc                                                                                     <  1:3685131/150‑1 (MQ=255)
                                                                 aTCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGcacc                                                                                  >  2:2198242/1‑150 (MQ=255)
                                                                    gCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACt                                                                               <  2:3601826/150‑1 (MQ=255)
                                                                    gCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACt                                                                               <  2:1427051/150‑1 (MQ=255)
                                                                    gCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACt                                                                               <  1:4120061/150‑1 (MQ=255)
                                                                       aaaGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtca                                                                            >  2:1899428/1‑150 (MQ=255)
                                                                        aaGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtcat                                                                           >  1:3553647/1‑150 (MQ=255)
                                                                         aGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtcatc                                                                          >  2:3860229/1‑150 (MQ=255)
                                                                         aGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtcatc                                                                          >  1:2595576/1‑150 (MQ=255)
                                                                          gCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtcatca                                                                         >  2:629247/1‑150 (MQ=255)
                                                                          gCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtcatca                                                                         >  2:4604445/1‑150 (MQ=255)
                                                                          gCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtcatca                                                                         >  2:1734603/1‑150 (MQ=255)
                                                                          gCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtcatca                                                                         <  2:2417683/150‑1 (MQ=255)
                                                                           ccTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCaa                                                                        >  2:1658496/1‑150 (MQ=255)
                                                                             tttGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATAGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGa                                                                      <  2:1577002/150‑1 (MQ=255)
                                                                                 cccGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGc                                                                  <  1:1534486/150‑1 (MQ=255)
                                                                                   cGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCtg                                                                <  1:702211/150‑1 (MQ=255)
                                                                                    ggCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCtgt                                                               <  2:767855/150‑1 (MQ=255)
                                                                                     gcgcGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCtgtg                                                              <  1:2112693/150‑1 (MQ=255)
                                                                                      cgcgTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTgg                                                             >  1:995647/1‑150 (MQ=255)
                                                                                      cgcgTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTgg                                                             <  1:715037/150‑1 (MQ=255)
                                                                                       gcgTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTggg                                                            <  1:1092436/150‑1 (MQ=255)
                                                                                        cgTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGt                                                           <  2:3635019/150‑1 (MQ=255)
                                                                                        cgTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGt                                                           >  1:1662215/1‑150 (MQ=255)
                                                                                        cgTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGt                                                           >  1:249545/1‑150 (MQ=255)
                                                                                          tAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCg                                                         >  2:4141764/1‑150 (MQ=255)
                                                                                            aaGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGAt                                                       >  1:573119/1‑150 (MQ=255)
                                                                                             aGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATc                                                      >  2:4031443/1‑150 (MQ=255)
                                                                                              ggTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATcc                                                     <  2:3901940/150‑1 (MQ=255)
                                                                                              ggTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATcc                                                     <  2:1937210/150‑1 (MQ=255)
                                                                                              ggTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATcc                                                     <  2:3901899/150‑1 (MQ=255)
                                                                                              ggTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATcc                                                     <  1:1018876/150‑1 (MQ=255)
                                                                                               gTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCg                                                    <  1:4652912/150‑1 (MQ=255)
                                                                                                   gTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGa                                                <  2:1246734/150‑1 (MQ=255)
                                                                                                          tATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCa                                         <  1:1358039/150‑1 (MQ=255)
                                                                                                            ttGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATg                                       <  1:2705307/150‑1 (MQ=255)
                                                                                                               aCCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCt                                    <  2:4428205/150‑1 (MQ=255)
                                                                                                                ccGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTc                                   <  1:4505476/150‑1 (MQ=255)
                                                                                                                ccGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTc                                   <  1:640155/150‑1 (MQ=255)
                                                                                                                  gAAAGTAACGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGgc                                            <  1:1727741/139‑1 (MQ=255)
                                                                                                                      gTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCa                             >  1:4106708/1‑150 (MQ=255)
                                                                                                                         aCGGTTGATTCATTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActg                          <  1:4009373/150‑1 (MQ=255)
                                                                                                                             ttGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtga                      <  2:844332/150‑1 (MQ=255)
                                                                                                                             ttGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtga                      <  2:2319275/150‑1 (MQ=255)
                                                                                                                               gATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaac                    <  1:616007/150‑1 (MQ=255)
                                                                                                                                   cATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtg                >  1:160729/1‑150 (MQ=255)
                                                                                                                                   cATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtg                <  1:3150677/150‑1 (MQ=255)
                                                                                                                                    aTTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtga               >  1:3295886/1‑150 (MQ=255)
                                                                                                                                    aTTGGTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtga               >  2:550955/1‑150 (MQ=255)
                                                                                                                                      tGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaac             >  1:2777225/1‑150 (MQ=255)
                                                                                                                                      tGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaac             >  1:3277679/1‑150 (MQ=255)
                                                                                                                                       ggTTTTTTGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaacg            <  1:1828692/150‑1 (MQ=255)
                                                                                                                                           ttttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGTTGAACGGGTg        <  2:1098402/150‑1 (MQ=255)
                                                                                                                                             ttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAt      >  2:1952378/1‑150 (MQ=255)
                                                                                                                                             ttGCCTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAt      <  1:3430471/150‑1 (MQ=255)
                                                                                                                                              tGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtc     <  2:2063760/150‑1 (MQ=255)
                                                                                                                                                ccTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcgt   >  1:3023935/1‑150 (MQ=255)
                                                                                                                                                ccTAGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcgt   >  2:24402/1‑150 (MQ=255)
                                                                                                                                                 cTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcgtc  <  2:3173879/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCCGGCGCGTAAAGGTGGGTCAGGTTATTGACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTC  >  CP009273/4308448‑4308742

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: