Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,285,974 C→A 5.4% 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.4% 1044.7 / 24.1 277D241Y (GAC→TAC) rsmI16S rRNA C1402 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (145/116);  new base A (3/12);  total (149/128)
Fisher's exact test for biased strand distribution p-value = 1.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 294 total aligned reads displayed.

CGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAACCACGT  >  CP009273/3285826‑3286123
                                                                                                                                                    |                                                                                                                                                     
cGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGtct                                                                                                                                                      >  1:2624770/1‑150 (MQ=255)
 ggCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGtctt                                                                                                                                                     >  2:4344668/1‑150 (MQ=255)
 ggCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGtctg                                                                                                                                                     >  2:1392255/1‑149 (MQ=255)
   cACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGtattct                                                                                                                                                   <  2:412800/150‑1 (MQ=255)
      tGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTcc                                                                                                                                                >  2:4455718/1‑150 (MQ=255)
             ccTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCtt                                                                                                                                         <  1:1100008/150‑1 (MQ=255)
              ctgctgCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCttt                                                                                                                                        <  1:692829/150‑1 (MQ=255)
                   gctCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGAc                                                                                                                                   <  2:1717967/150‑1 (MQ=255)
                   gctCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGAc                                                                                                                                   <  2:3098712/150‑1 (MQ=255)
                      ccAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCtt                                                                                                                                <  1:4308387/150‑1 (MQ=255)
                      ccAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCtt                                                                                                                                >  1:1007896/1‑150 (MQ=255)
                      ccAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCtt                                                                                                                                >  1:3735037/1‑150 (MQ=255)
                       cAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTc                                                                                                                               <  1:569850/150‑1 (MQ=255)
                        aGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCg                                                                                                                              >  2:3496746/1‑150 (MQ=255)
                           gcATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACa                                                                                                                           <  2:1218579/150‑1 (MQ=255)
                            cATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACaa                                                                                                                          >  2:3492095/1‑150 (MQ=255)
                            cATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACaa                                                                                                                          >  2:3491805/1‑150 (MQ=255)
                             aTACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAAt                                                                                                                         <  2:4046291/150‑1 (MQ=255)
                              tccTTATACAGCGCATTTTTCTTCACGCCGTGAATTTATGCGGCCAGCTGCGCCGCTTTTATCAGCGGAAGATCTGACAGTAGCAGAGCAAGCTAGCGCAGGGCATCGGCGGGTAAGTCTAAATCCAGTGCTTTATCAACTTCGACAATc                                                                                                                        >  2:2480702/3‑150 (MQ=255)
                              tACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATc                                                                                                                        >  1:4420094/1‑150 (MQ=255)
                               aCTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCa                                                                                                                       <  1:2893180/150‑1 (MQ=255)
                               aCTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCa                                                                                                                       <  2:3704977/150‑1 (MQ=255)
                               aCTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCa                                                                                                                       <  2:1821064/150‑1 (MQ=255)
                                cTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAg                                                                                                                      <  1:3268780/150‑1 (MQ=255)
                                cTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAg                                                                                                                      <  1:842588/150‑1 (MQ=255)
                                 ttATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGc                                                                                                                     >  2:3356249/1‑150 (MQ=255)
                                  tataCAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCa                                                                                                                    >  1:632544/1‑150 (MQ=255)
                                  tataCAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTATTCTTCCTGTGCTTTATGACCTTCGACCATCAGCa                                                                                                                    >  2:3648657/1‑150 (MQ=255)
                                   ataCAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCAc                                                                                                                   <  1:3839309/150‑1 (MQ=255)
                                      cAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCAt                                                                                                                >  1:3071403/1‑150 (MQ=255)
                                      cAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTATTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCAt                                                                                                                >  2:737158/1‑150 (MQ=255)
                                         cgcATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTc                                                                                                             >  1:1293908/1‑150 (MQ=255)
                                          gcATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCg                                                                                                            <  2:1282500/150‑1 (MQ=255)
                                            aTTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCTTCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGcc                                                                                                          <  2:4372203/150‑1 (MQ=255)
                                            aTTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCGGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGcc                                                                                                          >  1:514509/1‑150 (MQ=255)
                                              ttttCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCtt                                                                                                        <  2:667579/150‑1 (MQ=255)
                                                ttcttcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTg                                                                                                      <  1:2029117/150‑1 (MQ=255)
                                                  cttcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCg                                                                                                    >  2:1706969/1‑150 (MQ=255)
                                                   ttcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGa                                                                                                   >  1:3584574/1‑150 (MQ=255)
                                                   ttcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGa                                                                                                   >  1:2404780/1‑150 (MQ=255)
                                                    tcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGAc                                                                                                  >  2:2892431/1‑150 (MQ=255)
                                                    tcACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGAc                                                                                                  <  2:2589195/150‑1 (MQ=255)
                                                      aCGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACgg                                                                                                <  2:3514498/150‑1 (MQ=255)
                                                      aCGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACgg                                                                                                <  2:4416082/150‑1 (MQ=255)
                                                      aCGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACgg                                                                                                <  2:1242132/150‑1 (MQ=255)
                                                         ccGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGttt                                                                                             <  2:2758397/150‑1 (MQ=255)
                                                         ccGAGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGttt                                                                                             <  1:782754/150‑1 (MQ=255)
                                                           gTGAATTTCTGCGGCCTGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTc                                                                                           <  1:1542046/150‑1 (MQ=255)
                                                           gTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCTCCTTTGCGACGGTTTTc                                                                                           <  2:1640224/150‑1 (MQ=255)
                                                           gTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTc                                                                                           <  2:2426944/150‑1 (MQ=255)
                                                           gTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTc                                                                                           >  2:3289140/1‑150 (MQ=255)
                                                               aTTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCt                                                                                       >  1:1773919/1‑150 (MQ=255)
                                                               aTTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCt                                                                                       <  2:1098781/150‑1 (MQ=255)
                                                                tttCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCtt                                                                                      >  2:3916268/1‑150 (MQ=255)
                                                                tttCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCtt                                                                                      >  2:3812761/1‑150 (MQ=255)
                                                                tgtCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTATTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCtt                                                                                      <  2:2376811/148‑1 (MQ=255)
                                                                   cTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCt                                                                                   <  2:1696008/150‑1 (MQ=255)
                                                                          cAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCa                                                                            <  2:2867322/150‑1 (MQ=255)
                                                                              gccgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGcc                                                                        >  1:4082694/1‑150 (MQ=255)
                                                                               ccgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCca                                                                       >  2:2192242/1‑150 (MQ=255)
                                                                                cgccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCcag                                                                      >  2:2786166/1‑150 (MQ=255)
                                                                                 gccgcTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCAACATTTCGCCTTTGAGACGGTTTTCATCTTCCTTTACCCACGCcagc                                                                     >  1:1376170/1‑150 (MQ=255)
                                                                                     cTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTATTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCt                                                                 <  2:2211729/150‑1 (MQ=255)
                                                                                         tttCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGcc                                                             <  2:2860234/150‑1 (MQ=255)
                                                                                               cGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACggg                                                       >  1:4018682/1‑150 (MQ=255)
                                                                                                ggCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgc                                                      <  1:1060865/150‑1 (MQ=255)
                                                                                                ggCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgc                                                      >  1:1801353/1‑150 (MQ=255)
                                                                                                 gCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgcg                                                     >  2:3383336/1‑150 (MQ=255)
                                                                                                   aGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGgcgcg                                                   >  2:497464/1‑150 (MQ=255)
                                                                                                     ttCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGcc                                                 >  2:2730532/1‑150 (MQ=255)
                                                                                                      tCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCg                                                >  2:3870557/1‑150 (MQ=255)
                                                                                                         gCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGa                                             <  2:3157269/150‑1 (MQ=255)
                                                                                                         gCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTACTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGa                                             <  1:675047/150‑1 (MQ=255)
                                                                                                            tGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATg                                          <  1:1924601/150‑1 (MQ=255)
                                                                                                               agcagcGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGtt                                       >  1:2173451/1‑150 (MQ=255)
                                                                                                                    cgcCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCa                                  <  1:1706969/150‑1 (MQ=255)
                                                                                                                     gcCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAg                                 >  2:3770026/1‑150 (MQ=255)
                                                                                                                     gcCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAg                                 <  1:590613/150‑1 (MQ=255)
                                                                                                                      ccAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCGTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAgg                                <  2:129128/150‑1 (MQ=255)
                                                                                                                        aGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGtt                              >  2:2951688/1‑150 (MQ=255)
                                                                                                                        aGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGtt                              >  2:2951575/1‑150 (MQ=255)
                                                                                                                           gTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTg                           >  2:429779/1‑150 (MQ=255)
                                                                                                                             gcgcAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGt                         <  1:274074/150‑1 (MQ=255)
                                                                                                                                cAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAg                      <  1:3220450/150‑1 (MQ=255)
                                                                                                                                  gggCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCt                    >  2:432085/1‑150 (MQ=255)
                                                                                                                                  gggCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCt                    >  2:2017251/1‑150 (MQ=255)
                                                                                                                                  gggCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCt                    >  2:2017215/1‑150 (MQ=255)
                                                                                                                                   ggCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTc                   <  1:3339563/150‑1 (MQ=255)
                                                                                                                                   ggCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTc                   <  2:4435225/150‑1 (MQ=255)
                                                                                                                                   ggCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTc                   <  1:3338555/150‑1 (MQ=255)
                                                                                                                                    gCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCa                  <  1:1127517/150‑1 (MQ=255)
                                                                                                                                     cATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCAc                 >  2:3238937/1‑150 (MQ=255)
                                                                                                                                      aTCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCAcg                >  1:3508440/1‑150 (MQ=255)
                                                                                                                                           cggGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCa           >  2:1384891/1‑150 (MQ=255)
                                                                                                                                            gggTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAg          >  1:733951/1‑150 (MQ=255)
                                                                                                                                                 aGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAACCa     >  2:4217754/1‑150 (MQ=255)
                                                                                                                                                  gTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAACCAc    >  1:4351194/1‑150 (MQ=255)
                                                                                                                                                   tcttctTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAACCACg   <  2:1115580/150‑1 (MQ=255)
                                                                                                                                                    cttcttCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAACCACGt  >  1:125985/1‑150 (MQ=255)
                                                                                                                                                    cttcttCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAACCACGt  >  1:4049633/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTTCAGCGGCAGTTCTGCCTGTAGCAGCGCCAGCGTGCGCAGGGCATCGGCGGGTAAGTCTTCTTCCTGTGCTTTATGACCTTCGACAATCAGCACCATTTCGCCTTTGCGACGGTTTTCATCTTCCTTTACCCACGCCAGCAGCTCGCCAACGGGCGCGCCGTGAATGGTTTCCCAGGTTTTGGTCAGCTCACGCGCCAGAACCACGT  >  CP009273/3285826‑3286123

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: ‑