Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,003,597 G→T 100% 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,5970GT98.8% 781.6 / 0.5 240D100Y (GAC→TAC) rlmLfused 23S rRNA m(2)G2445 and m(7)G2069 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (114/123);  total (115/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 240 total aligned reads displayed.

GCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAA  >  CP009273/1003451‑1003738
                                                                                                                                                  |                                                                                                                                             
gccgccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTACGCTGTCCACTTCAGTGGTTTGAATTac                                                                                                                                             >  2:3261590/1‑149 (MQ=255)
    ccTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCa                                                                                                                                         >  1:1517470/1‑149 (MQ=255)
          cTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCa                                                                                                                                   >  2:513420/1‑149 (MQ=255)
           tCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCaa                                                                                                                                  >  2:2087439/1‑149 (MQ=255)
           tCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCaa                                                                                                                                  >  2:1897998/1‑149 (MQ=255)
            cgcgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAc                                                                                                                                 <  1:512760/149‑1 (MQ=255)
              cgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcag                                                                                                                               <  2:1786535/149‑1 (MQ=255)
              cgTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGAGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcag                                                                                                                               <  1:1094658/149‑1 (MQ=255)
                tattatGTTGCCGGTGGTCGAGTGTAAGGTTTACATCGATTTAGACCTCTATCTGGGTGTTCAGTCGATCAACTGGACATAGATGTTTATTCCTGGCGCGTCCTTGGCTGTCCACTTCAGTGGTTTGAATTACACCATCGGCAAcagtc                                                                                                                             <  1:3606851/149‑1 (MQ=255)
                tattatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtc                                                                                                                             <  1:844465/149‑1 (MQ=255)
                tattatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtc                                                                                                                             <  2:910959/149‑1 (MQ=255)
                  ttatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtcag                                                                                                                           >  1:3732195/1‑149 (MQ=255)
                   tatGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAAcagtcagt                                                                                                                          >  1:2068177/1‑149 (MQ=255)
                     tGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTAc                                                                                                                        <  1:1762657/149‑1 (MQ=255)
                        tGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGgc                                                                                                                     <  1:203995/149‑1 (MQ=255)
                               gggCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtga                                                                                                              <  2:3142772/149‑1 (MQ=255)
                               gggCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtga                                                                                                              >  2:1012532/1‑149 (MQ=255)
                               gggCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtga                                                                                                              >  1:3879271/1‑149 (MQ=255)
                               gggCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtga                                                                                                              >  2:1012292/1‑149 (MQ=255)
                                ggCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaa                                                                                                             >  2:2668185/1‑149 (MQ=255)
                                   gAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCAATTCAGTTGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagt                                                                                                          <  1:3786702/149‑1 (MQ=255)
                                    aGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtg                                                                                                         <  2:3832700/149‑1 (MQ=255)
                                    aGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtg                                                                                                         <  2:724672/149‑1 (MQ=255)
                                       gtAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGAtgaaagtgaaa                                                                                                      <  2:3339116/149‑1 (MQ=255)
                                         aaGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGa                                                                                                    <  1:3582627/149‑1 (MQ=255)
                                          aGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAc                                                                                                   >  2:3843846/1‑149 (MQ=255)
                                          aGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAc                                                                                                   <  1:256323/149‑1 (MQ=255)
                                           ggTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAcg                                                                                                  >  1:2650866/1‑149 (MQ=255)
                                            gTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGAcgc                                                                                                 >  1:1511556/1‑149 (MQ=255)
                                              ttACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGa                                                                                               <  1:964025/149‑1 (MQ=255)
                                               tACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGAt                                                                                              <  2:467759/149‑1 (MQ=255)
                                                     gATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGa                                                                                        <  1:3263349/149‑1 (MQ=255)
                                                      aTTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGAt                                                                                       >  1:2959017/1‑149 (MQ=255)
                                                       tttAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATg                                                                                      >  1:1053070/1‑149 (MQ=255)
                                                        ttAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGc                                                                                     <  2:2497165/149‑1 (MQ=255)
                                                           gACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCttt                                                                                  >  1:560508/1‑149 (MQ=255)
                                                           gACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCttt                                                                                  >  1:560746/1‑149 (MQ=255)
                                                              ctctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAc                                                                               >  1:2055806/1‑149 (MQ=255)
                                                              ctctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAc                                                                               <  1:425190/149‑1 (MQ=255)
                                                              ctctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAc                                                                               <  1:2667211/149‑1 (MQ=255)
                                                               tctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAcg                                                                              >  2:4094787/1‑149 (MQ=255)
                                                               tctATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCAcg                                                                              >  2:3574857/1‑149 (MQ=255)
                                                                  aTCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCgg                                                                           >  2:431871/1‑149 (MQ=255)
                                                                   tctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGa                                                                          >  2:2368925/1‑149 (MQ=255)
                                                                   tctcGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGa                                                                          >  1:3634472/1‑149 (MQ=255)
                                                                      cGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaa                                                                       <  2:1186798/149‑1 (MQ=255)
                                                                       ggTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaaa                                                                      <  1:3461252/149‑1 (MQ=255)
                                                                       ggTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGaaaaa                                                                      >  2:3362097/1‑149 (MQ=255)
                                                                         tgtTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATc                                                                    >  2:1715523/1‑149 (MQ=255)
                                                                           ttCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTg                                                                  <  2:1877894/149‑1 (MQ=255)
                                                                             cAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAGAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGcc                                                                >  2:1232910/1‑149 (MQ=255)
                                                                             cAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGcc                                                                >  2:1437944/1‑149 (MQ=255)
                                                                             cAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGcc                                                                >  2:2025517/1‑149 (MQ=255)
                                                                               ggCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCcgc                                                              <  2:3826333/149‑1 (MQ=255)
                                                                                   aTCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTcc                                                          >  2:3935418/1‑149 (MQ=255)
                                                                                      aaCTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCaaa                                                       <  2:2433642/149‑1 (MQ=255)
                                                                                         tGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGt                                                    >  1:634311/1‑149 (MQ=255)
                                                                                          ggACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGtt                                                   <  1:158263/149‑1 (MQ=255)
                                                                                          ggACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGtt                                                   <  1:2313910/149‑1 (MQ=255)
                                                                                          ggACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGtt                                                   <  1:3517167/149‑1 (MQ=255)
                                                                                            aCAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGa                                                 <  2:625522/149‑1 (MQ=255)
                                                                                             cAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGAt                                                <  2:2651507/149‑1 (MQ=255)
                                                                                                agaTGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATcgc                                             <  2:3630556/149‑1 (MQ=255)
                                                                                                    gTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATg                                         >  1:713458/1‑149 (MQ=255)
                                                                                                          tCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCgg                                   >  1:1447216/1‑149 (MQ=255)
                                                                                                            cTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAACAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGat                                 >  2:2854548/1‑149 (MQ=255)
                                                                                                             tGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGata                                >  2:2930112/1‑149 (MQ=255)
                                                                                                             tGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGata                                <  2:449952/149‑1 (MQ=255)
                                                                                                              ggCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGatat                               <  2:3704257/149‑1 (MQ=255)
                                                                                                              ggCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGatat                               <  2:1523059/149‑1 (MQ=255)
                                                                                                              ggCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGatat                               >  1:1875833/1‑149 (MQ=255)
                                                                                                              ggCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGatat                               <  1:391981/149‑1 (MQ=255)
                                                                                                               gcgcGACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATc                              <  1:622517/149‑1 (MQ=255)
                                                                                                               gcgcGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATc                              <  1:806960/149‑1 (MQ=255)
                                                                                                                 gcgACCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCcg                            <  1:519915/149‑1 (MQ=255)
                                                                                                                    aCCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGt                         <  2:1396929/149‑1 (MQ=255)
                                                                                                                    aCCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGt                         >  1:1437287/1‑149 (MQ=255)
                                                                                                                    aCCTTCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGt                         >  1:1438058/1‑149 (MQ=255)
                                                                                                                        tCGCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAAc                     <  2:1338994/149‑1 (MQ=255)
                                                                                                                          gCTGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGt                   >  1:1004057/1‑149 (MQ=255)
                                                                                                                            tGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCt                 <  1:3515494/149‑1 (MQ=255)
                                                                                                                            tGTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCt                 >  2:395419/1‑149 (MQ=255)
                                                                                                                             gTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTg                >  1:961987/1‑149 (MQ=255)
                                                                                                                             gTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTg                >  1:3018584/1‑149 (MQ=255)
                                                                                                                             gTCCACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTg                >  2:821702/1‑149 (MQ=255)
                                                                                                                                cACTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCt             <  1:3210395/149‑1 (MQ=255)
                                                                                                                                 aCTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTg            >  1:2306307/1‑149 (MQ=255)
                                                                                                                                 aCTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTg            >  1:2480597/1‑149 (MQ=255)
                                                                                                                                 aCTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTg            >  1:1664460/1‑149 (MQ=255)
                                                                                                                                 aCTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTg            <  2:145595/149‑1 (MQ=255)
                                                                                                                                  cTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGc           >  1:2129679/1‑149 (MQ=255)
                                                                                                                                  cTTCAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGc           >  1:1211869/1‑149 (MQ=255)
                                                                                                                                     cAGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATa        >  2:61281/1‑149 (MQ=255)
                                                                                                                                      aGTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaa       <  2:164331/149‑1 (MQ=255)
                                                                                                                                       gTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGACAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTGAACGTCTGGCTGCATata      >  2:2325311/1‑147 (MQ=255)
                                                                                                                                       gTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaaa      <  1:1240324/149‑1 (MQ=255)
                                                                                                                                       gTGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATaaa      >  2:3206610/1‑149 (MQ=255)
                                                                                                                                        tGGTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAg     <  2:2104407/149‑1 (MQ=255)
                                                                                                                                          gTTTGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaa   <  2:1159480/149‑1 (MQ=255)
                                                                                                                                           tttGAATTACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGaaa  >  2:260170/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GCCGCCTGGCCTCGCGTATTATGTTGCCGCTGGGCGAGTGTAAGGTTTACAGCGATTTAGACCTCTATCTCGGTGTTCAGGCGATCAACTGGACAGAGATGTTTAATCCTGGCGCGACCTTCGCTGTCCACTTCAGTGGTTTGAATGACACCATCCGCAACAGTCAGTACGGCGCGATGAAAGTGAAAGACGCGATCGTCGATGCTTTCACGCGGAAAAATCTGCCGCGTCCAAATGTTGATCGCGATGCGCCGGATATCCGCGTTAACGTCTGGCTGCATAAAGAAA  >  CP009273/1003451‑1003738

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: