Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,285,974 C→A 5.2% D241Y (GAC→TAC)  rsmI ← 16S rRNA C1402 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,285,9740CA 5.2% 1251.7 / 27.5 324D241Y (GAC→TAC) rsmI16S rRNA C1402 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (170/136);  new base A (11/7);  total (181/143)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 340 total aligned reads displayed.

TCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAACCAC  >  CP009273/3285825‑3286121
                                                                                                                                                     |                                                                                                                                                   
    cACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGtcttct                                                                                                                                                 >  1:3580731/1‑150 (MQ=255)
    cACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGtcttat                                                                                                                                                 <  2:3875458/150‑3 (MQ=255)
         ttACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCtg                                                                                                                                            <  2:3536599/150‑1 (MQ=255)
            ccccTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGc                                                                                                                                         >  1:1444707/1‑150 (MQ=255)
             cccTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCt                                                                                                                                        >  2:3596634/1‑150 (MQ=255)
             cccTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCt                                                                                                                                        <  1:3166399/150‑1 (MQ=255)
               ctgctgCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCttt                                                                                                                                      <  1:4314929/150‑1 (MQ=255)
                tgctgcTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTa                                                                                                                                     >  2:3690434/1‑150 (MQ=255)
                 gctgctCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTAt                                                                                                                                    <  1:3069817/150‑1 (MQ=255)
                 gctgctCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTATTCTTCCTGTGCTTTAt                                                                                                                                    >  1:2321182/1‑150 (MQ=255)
                  ctgctCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCCTCGGCGGGTAAGTCTTCTTCCTGTGCTTTATg                                                                                                                                   >  2:832969/1‑150 (MQ=255)
                  ctgctCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATg                                                                                                                                   >  2:1917273/1‑150 (MQ=255)
                    gctCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGAc                                                                                                                                 >  2:162158/1‑150 (MQ=255)
                          gcgcATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGa                                                                                                                           >  1:4905645/1‑150 (MQ=255)
                                  ttAGACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGc                                                                                                                   <  2:4278329/150‑1 (MQ=255)
                                      aCAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCa                                                                                                               <  2:2088138/150‑1 (MQ=255)
                                        aGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCAtt                                                                                                             >  1:3628988/1‑150 (MQ=255)
                                         gcgcATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCAttt                                                                                                            <  1:4856551/150‑1 (MQ=255)
                                          cgcATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTc                                                                                                           <  1:596288/150‑1 (MQ=255)
                                          cgcATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTATTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTc                                                                                                           <  2:2576073/150‑1 (MQ=255)
                                           gcATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCg                                                                                                          >  2:2942055/1‑150 (MQ=255)
                                           gcATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGGAAGTCTTCGTCCTGTGCTTTATGACCGTCGACAATCAGCACCATTTCg                                                                                                          >  1:3299277/1‑150 (MQ=255)
                                            cATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGc                                                                                                         >  2:1559369/1‑150 (MQ=255)
                                            cATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGc                                                                                                         >  1:1341130/1‑150 (MQ=255)
                                             aTTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGcc                                                                                                        >  2:4347732/1‑150 (MQ=255)
                                              tttttCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCt                                                                                                       <  2:1697479/150‑1 (MQ=255)
                                                 ttcttcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTg                                                                                                    <  2:175160/150‑1 (MQ=255)
                                                 ttcttcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTg                                                                                                    <  2:2239757/150‑1 (MQ=255)
                                                  tcttcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTATTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGc                                                                                                   >  1:1368938/1‑150 (MQ=255)
                                                   cttcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCg                                                                                                  <  2:4117054/150‑1 (MQ=255)
                                                      cACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACg                                                                                               <  2:1340946/150‑1 (MQ=255)
                                                      cACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACg                                                                                               <  2:2420902/150‑1 (MQ=255)
                                                       aCGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACgg                                                                                              >  1:2422650/1‑150 (MQ=255)
                                                        cGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTGTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGt                                                                                             >  1:3635003/1‑150 (MQ=255)
                                                          ccGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGttg                                                                                           >  1:2680186/1‑149 (MQ=255)
                                                           cGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGtttt                                                                                          >  2:2248563/1‑150 (MQ=255)
                                                           cGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGtttt                                                                                          <  2:4905645/150‑1 (MQ=255)
                                                           cGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGtttt                                                                                          >  2:4590199/1‑150 (MQ=255)
                                                            gTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTc                                                                                         <  1:977184/150‑1 (MQ=255)
                                                             tGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTGTGCGACGGTTTTCa                                                                                        >  2:2697247/1‑150 (MQ=255)
                                                               aaTTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGAGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTTCGACGGTTTTCATc                                                                                      >  1:4157640/1‑150 (MQ=255)
                                                               aaTTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGAGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATc                                                                                      >  1:3096713/1‑150 (MQ=255)
                                                                  ttCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTc                                                                                   >  2:2472246/1‑150 (MQ=255)
                                                                    cTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCt                                                                                 >  1:1853336/1‑150 (MQ=255)
                                                                     tGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCtt                                                                                >  2:2811297/1‑150 (MQ=255)
                                                                     tGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCtt                                                                                <  2:327458/150‑1 (MQ=255)
                                                                     tGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTATTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCtt                                                                                >  1:2106867/1‑150 (MQ=255)
                                                                      gCGGCCAGCGCCGCCGCTTTTTTCTGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGCCCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCtgt                                                                               <  1:1006711/150‑3 (MQ=255)
                                                                       cGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTa                                                                              <  1:4606675/150‑1 (MQ=255)
                                                                       cGGCCAGCGCCGCCGCTTTTTTAAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTa                                                                              >  1:150332/1‑150 (MQ=255)
                                                                         gCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCAATTCGCCTTTGCGACGGTTTTCATCTTCCTTTAcc                                                                            >  1:1388227/1‑150 (MQ=255)
                                                                           cAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCa                                                                          <  1:1053166/150‑1 (MQ=255)
                                                                               gccgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGcc                                                                      <  1:753554/150‑1 (MQ=255)
                                                                               gccgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGcc                                                                      >  2:57890/1‑150 (MQ=255)
                                                                               gccgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGcc                                                                      <  1:1196858/150‑1 (MQ=255)
                                                                               gccgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGcc                                                                      <  1:2833276/150‑1 (MQ=255)
                                                                               gccgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTATTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGcc                                                                      <  2:1770550/150‑1 (MQ=255)
                                                                                ccgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCca                                                                     <  1:4978782/150‑1 (MQ=255)
                                                                                 cgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCcag                                                                    <  2:3931153/150‑1 (MQ=255)
                                                                                   ccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCcagca                                                                  <  1:5035516/150‑1 (MQ=255)
                                                                                       ttttttCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTc                                                              <  1:3404414/150‑1 (MQ=255)
                                                                                       ttttttCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTc                                                              <  2:306427/150‑1 (MQ=255)
                                                                                           ttCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCa                                                          >  2:4015878/1‑150 (MQ=255)
                                                                                            tCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCaa                                                         >  1:2166277/1‑150 (MQ=255)
                                                                                             cAGCGGCAGTTCTGCCTGTAGCAGCTCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAAc                                                        <  2:4833495/150‑1 (MQ=255)
                                                                                                   cAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgcgc                                                  >  1:4718243/1‑150 (MQ=255)
                                                                                                    aGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgcgcg                                                 >  2:2268082/1‑150 (MQ=255)
                                                                                                     gTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgcgcgc                                                >  1:3252438/1‑150 (MQ=255)
                                                                                                         tGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTg                                            >  1:253607/1‑150 (MQ=255)
                                                                                                         tGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTg                                            <  1:2801960/150‑1 (MQ=255)
                                                                                                         tGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTg                                            <  1:1647052/150‑1 (MQ=255)
                                                                                                          gCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGa                                           >  2:2003661/1‑150 (MQ=255)
                                                                                                             tGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACGATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATg                                        >  1:818661/1‑150 (MQ=255)
                                                                                                             tGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATg                                        >  2:1023239/1‑150 (MQ=255)
                                                                                                              gTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATgg                                       >  2:1634311/1‑150 (MQ=255)
                                                                                                                agcagcGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGtt                                     <  2:1853336/150‑1 (MQ=255)
                                                                                                                  cagcGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTc                                   <  2:1932218/150‑1 (MQ=255)
                                                                                                                   agcGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTcc                                  >  1:3033670/1‑150 (MQ=255)
                                                                                                                   agcGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTcc                                  <  2:3764093/150‑1 (MQ=255)
                                                                                                                    gcgcCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTccc                                 <  2:2166277/150‑1 (MQ=255)
                                                                                                                     cgcCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCa                                <  2:5039114/150‑1 (MQ=255)
                                                                                                                     cgcCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCa                                <  1:2500269/150‑1 (MQ=255)
                                                                                                                      gcCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAg                               >  1:3390487/1‑150 (MQ=255)
                                                                                                                      gcCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAg                               <  2:3796043/150‑1 (MQ=255)
                                                                                                                       ccAGCGAGCGCAGGGCATCGGCGGGTAACAATTCTTCAAGTGCATTAAGACCTTCGACAATCAGCACCAAATCGCCTTTGAGACGGTTTTCATCTTCCAGTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGATTCCCAgg                              >  2:404312/1‑150 (MQ=255)
                                                                                                                        cAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGt                             >  1:3769087/1‑150 (MQ=255)
                                                                                                                           cGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGtttt                          >  1:2642573/1‑150 (MQ=255)
                                                                                                                             tGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTgg                        <  1:4237554/150‑1 (MQ=255)
                                                                                                                                  aGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGc                   >  1:2459133/1‑150 (MQ=255)
                                                                                                                                      cATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCAc               >  1:4118784/1‑150 (MQ=255)
                                                                                                                                       aTCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCAcg              <  1:3888280/150‑1 (MQ=255)
                                                                                                                                       aTCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCAcg              >  1:3486685/1‑150 (MQ=255)
                                                                                                                                         cggcggGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCAcgcg            <  1:3924418/150‑1 (MQ=255)
                                                                                                                                         cggcggGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCAcgcg            >  1:1043666/1‑150 (MQ=255)
                                                                                                                                          ggcggGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACgcgc           >  1:1921256/1‑150 (MQ=255)
                                                                                                                                            cggGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCa         >  2:3454281/1‑150 (MQ=255)
                                                                                                                                            cggGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCa         <  1:791451/150‑1 (MQ=255)
                                                                                                                                              ggTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGa       <  1:4442003/150‑1 (MQ=255)
                                                                                                                                                tAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAAc     >  2:672842/1‑150 (MQ=255)
                                                                                                                                                 aaGTATTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAAcc    >  2:724314/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAACCAC  >  CP009273/3285825‑3286121

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: