Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,003,597 G→T 8.7% D100Y (GAC→TAC)  rlmL → fused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,003,5970GT 8.7% 941.8 / 38.6 252D100Y (GAC→TAC) rlmLfused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (129/100);  new base T (8/15);  total (137/115)
Fisher's exact test for biased strand distribution p-value = 7.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

GGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGCCAGT  >  CP009273/1003448‑1003746
                                                                                                                                                     |                                                                                                                                                     
ggAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATg                                                                                                                                                       <  2:2038609/150‑1 (MQ=255)
 gAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGa                                                                                                                                                      <  2:581198/150‑1 (MQ=255)
     cgccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACAcc                                                                                                                                                  <  1:1467049/150‑1 (MQ=255)
      gccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCa                                                                                                                                                 <  1:1870675/150‑1 (MQ=255)
        cTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATc                                                                                                                                               >  2:1626057/1‑150 (MQ=255)
         tGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATcc                                                                                                                                              >  1:109439/1‑150 (MQ=255)
            ccTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCa                                                                                                                                           >  1:845845/1‑150 (MQ=255)
             cTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCaa                                                                                                                                          <  1:1488184/150‑1 (MQ=255)
                 cgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAAcagt                                                                                                                                      >  1:1399907/1‑150 (MQ=255)
                 cgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAAcagt                                                                                                                                      >  2:3771096/1‑150 (MQ=255)
                  gTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAAcagtc                                                                                                                                     <  2:3163898/150‑1 (MQ=255)
                  gTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAAcagtc                                                                                                                                     >  2:2787961/1‑150 (MQ=255)
                    attatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAAcagtcag                                                                                                                                   >  1:3609968/1‑150 (MQ=255)
                      tatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTa                                                                                                                                 >  1:4505048/1‑150 (MQ=255)
                        tGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACg                                                                                                                               >  1:4297408/1‑150 (MQ=255)
                              cGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGa                                                                                                                         >  1:356163/1‑150 (MQ=255)
                                cTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGAtg                                                                                                                       >  1:2950300/1‑150 (MQ=255)
                                 tGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGAtga                                                                                                                      >  2:4550585/1‑150 (MQ=255)
                                   ggCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGAtgaaa                                                                                                                    >  2:1417465/1‑150 (MQ=255)
                                   ggCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGAtgaaa                                                                                                                    <  1:2185835/150‑1 (MQ=255)
                                      gAGTGTAAGTTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtg                                                                                                                 <  2:94153/150‑1 (MQ=255)
                                       aGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtga                                                                                                                >  2:3934773/1‑150 (MQ=255)
                                       aGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtga                                                                                                                >  2:2638468/1‑150 (MQ=255)
                                        gtgtAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaa                                                                                                               <  2:4419958/150‑1 (MQ=255)
                                          gtAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaag                                                                                                             <  1:719408/150‑1 (MQ=255)
                                          gtAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaag                                                                                                             >  2:1799642/1‑150 (MQ=255)
                                            aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAc                                                                                                           <  2:3458925/150‑1 (MQ=255)
                                            aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAc                                                                                                           <  1:3608983/150‑1 (MQ=255)
                                            aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAc                                                                                                           <  1:4119352/150‑1 (MQ=255)
                                            aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAc                                                                                                           <  2:3603003/150‑1 (MQ=255)
                                               gTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAcgcg                                                                                                        >  2:2990098/1‑150 (MQ=255)
                                                tttACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGa                                                                                                       <  1:2644949/150‑1 (MQ=255)
                                                 ttACAGCGAGTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAt                                                                                                      >  1:3824853/1‑150 (MQ=255)
                                                  tACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtc                                                                                                     <  1:3895426/150‑1 (MQ=255)
                                                  tACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtc                                                                                                     <  1:1692883/150‑1 (MQ=255)
                                                   aCAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtcg                                                                                                    <  1:2209024/150‑1 (MQ=255)
                                                     aGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAtcgtc                                                                                                  >  2:768430/1‑150 (MQ=255)
                                                        gATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGAt                                                                                               >  1:1373995/1‑150 (MQ=255)
                                                          tttAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGc                                                                                             >  2:2257372/1‑150 (MQ=255)
                                                            tAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCtt                                                                                           >  2:1137430/1‑150 (MQ=255)
                                                             aGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCttt                                                                                          >  2:1741160/1‑150 (MQ=255)
                                                             aGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCttt                                                                                          <  2:3666291/150‑1 (MQ=255)
                                                              gACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTc                                                                                         <  2:1187175/150‑1 (MQ=255)
                                                              gACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTc                                                                                         <  2:3505659/150‑1 (MQ=255)
                                                                ccTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAc                                                                                       <  1:2310527/150‑1 (MQ=255)
                                                                   ctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAcgcg                                                                                    >  2:3322827/1‑150 (MQ=255)
                                                                    tATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCgg                                                                                   <  1:2156914/150‑1 (MQ=255)
                                                                      tctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaa                                                                                 <  1:1585864/150‑1 (MQ=255)
                                                                       ctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaa                                                                                >  2:716945/1‑150 (MQ=255)
                                                                         cGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaaa                                                                              <  2:4381926/150‑1 (MQ=255)
                                                                          ggTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAAt                                                                             <  1:408104/150‑1 (MQ=255)
                                                                            tgtTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCt                                                                           <  1:2458169/150‑1 (MQ=255)
                                                                            tgtTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCt                                                                           >  2:1688440/1‑150 (MQ=255)
                                                                              ttCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGc                                                                         >  2:2620355/1‑150 (MQ=255)
                                                                                  ggCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcgcg                                                                     <  2:2421324/150‑1 (MQ=255)
                                                                                  ggCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcgcg                                                                     <  1:2029025/150‑1 (MQ=255)
                                                                                    cGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTc                                                                   <  2:1429105/150‑1 (MQ=255)
                                                                                     gATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTcc                                                                  >  2:1043146/1‑150 (MQ=255)
                                                                                     gATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTcc                                                                  >  1:3466947/1‑150 (MQ=255)
                                                                                           cTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGt                                                            >  2:4024428/1‑150 (MQ=255)
                                                                                            tGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGtt                                                           <  1:2508295/150‑1 (MQ=255)
                                                                                             ggACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTg                                                          >  1:557358/1‑150 (MQ=255)
                                                                                               aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGAt                                                        >  1:278684/1‑150 (MQ=255)
                                                                                               aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGAt                                                        >  1:3732584/1‑150 (MQ=255)
                                                                                               aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGAt                                                        >  1:2921006/1‑150 (MQ=255)
                                                                                                 agagATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATcg                                                      <  2:2527792/150‑1 (MQ=255)
                                                                                                 agagATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATcg                                                      >  1:1291494/1‑150 (MQ=255)
                                                                                                  gagaTGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATcgc                                                     >  1:3163551/1‑150 (MQ=255)
                                                                                                  gagaTGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATcgc                                                     >  1:3754307/1‑150 (MQ=255)
                                                                                                           aaTCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCg                                            >  2:4303561/1‑150 (MQ=255)
                                                                                                            aTCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCgg                                           >  2:3053547/1‑150 (MQ=255)
                                                                                                                     cgACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCcgcg                                  >  1:2947537/1‑150 (MQ=255)
                                                                                                                       aCCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGtt                                >  2:2396858/1‑150 (MQ=255)
                                                                                                                        ccTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTa                               <  1:4388160/150‑1 (MQ=255)
                                                                                                                           tCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACg                            <  1:3154792/150‑1 (MQ=255)
                                                                                                                              cTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCt                         <  1:3380009/150‑1 (MQ=255)
                                                                                                                              cTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCt                         <  2:669466/150‑1 (MQ=255)
                                                                                                                              cTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCt                         <  2:3286768/150‑1 (MQ=255)
                                                                                                                                gTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTgg                       >  1:3125667/1‑150 (MQ=255)
                                                                                                                                gTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTgg                       >  1:3344023/1‑150 (MQ=255)
                                                                                                                                gTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTgg                       <  1:798379/150‑1 (MQ=255)
                                                                                                                                    aCTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGc                   >  1:2686789/1‑150 (MQ=255)
                                                                                                                                      ttCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCAt                 <  2:4297408/150‑1 (MQ=255)
                                                                                                                                       tCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATa                <  1:3236100/150‑1 (MQ=255)
                                                                                                                                       tCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATa                >  1:3092413/1‑150 (MQ=255)
                                                                                                                                        cAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaa               >  2:366450/1‑150 (MQ=255)
                                                                                                                                         aGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaaa              <  1:132952/150‑1 (MQ=255)
                                                                                                                                          gTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAg             >  2:2090263/1‑150 (MQ=255)
                                                                                                                                            ggTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa           >  1:2759472/1‑150 (MQ=255)
                                                                                                                                            ggTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa           >  1:2007164/1‑150 (MQ=255)
                                                                                                                                             gTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCTCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaaa          >  2:1968979/1‑150 (MQ=255)
                                                                                                                                             gTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaaa          <  1:3588773/150‑1 (MQ=255)
                                                                                                                                               ttGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAAcc        <  1:3238642/150‑1 (MQ=255)
                                                                                                                                               ttGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAAcc        <  1:1044909/150‑1 (MQ=255)
                                                                                                                                                  aaTGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGcc     <  2:2127481/150‑1 (MQ=255)
                                                                                                                                                   aTGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGCCa    <  2:3584797/150‑1 (MQ=255)
                                                                                                                                                   aTGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGCCa    >  2:620636/1‑150 (MQ=255)
                                                                                                                                                   aTGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGCCa    >  1:2824694/1‑150 (MQ=255)
                                                                                                                                                    tGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGCCAg   <  1:628487/150‑1 (MQ=255)
                                                                                                                                                     gACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGCCAGt  >  1:2742757/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAGCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAACCGCCAGT  >  CP009273/1003448‑1003746

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: