Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,279,042 C→G 11.3% intergenic (+67/+39) rlmG → / ← PP_4649 ribosomal RNA large subunit methyltransferase G/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,279,0420CG11.3% 612.4 / 15.3 205intergenic (+67/+39)rlmG/PP_4649ribosomal RNA large subunit methyltransferase G/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (111/67);  new base G (2/25);  total (113/92)
Fisher's exact test for biased strand distribution p-value = 5.75e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

ACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCT  >  NC_002947/5278893‑5279190
                                                                                                                                                       |                                                                                                                                                      
acaTCAAGGTGGCGGGGTTGTTCCGGGGGGTGTAGCAGGTGGCGTAAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAgcgc                                                                                                                                                        <  1:287444/150‑1 (MQ=39)
acaGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAgcgc                                                                                                                                                        >  2:2753452/1‑150 (MQ=255)
acaGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTAGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAgcgc                                                                                                                                                        <  2:4678649/150‑1 (MQ=255)
  aGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGt                                                                                                                                                      <  1:4621909/150‑1 (MQ=255)
  aGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGt                                                                                                                                                      >  1:2612966/1‑150 (MQ=255)
  aGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGt                                                                                                                                                      >  1:1662258/1‑150 (MQ=255)
  aGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCACGt                                                                                                                                                      >  1:6112083/1‑150 (MQ=255)
   gCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGtt                                                                                                                                                     <  2:1531870/150‑1 (MQ=255)
    cAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCGAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTg                                                                                                                                                    >  2:2884383/1‑150 (MQ=255)
    cAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTg                                                                                                                                                    >  1:3697292/1‑150 (MQ=255)
     aaGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGc                                                                                                                                                   <  2:3022791/150‑1 (MQ=255)
      atctgggcgatgagcgtcgcagcgctgggtgtggcggcctgtgtgatgcggctgttcatgggggcggcgggcctgaccacacactgataaagagcggcAGCGGGCTTGCCCACGAAGCAGGCAACGCCTGTGCCATCTAGCGCGTTGCCt                                                                                                                                                 >  1:2139045‑M1/99‑150 (MQ=255)
      aGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGcc                                                                                                                                                  <  2:6855489/150‑1 (MQ=255)
         tGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGc                                                                                                                                              <  2:6218571/150‑1 (MQ=255)
         tGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGc                                                                                                                                              >  2:5043097/1‑150 (MQ=255)
           gcgcGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCtt                                                                                                                                            >  1:3408962/1‑150 (MQ=255)
             gcgGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCg                                                                                                                                          >  2:2026195/1‑150 (MQ=255)
              cgGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTgtgttcg                                                                                                                                          <  2:1692375/150‑7 (MQ=255)
               ggTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTg                                                                                                                                        <  1:3375083/150‑1 (MQ=255)
               ggTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTg                                                                                                                                        <  1:5054131/150‑1 (MQ=255)
               ggTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTg                                                                                                                                        <  1:6393849/150‑1 (MQ=255)
               ggTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTg                                                                                                                                        >  2:6502402/1‑150 (MQ=255)
               atgagcgtcgcagcgctgggtgtggcggcctgtgtgatgcggctgttcatgggggcggcgggcctgaccacacactgataaagagcggcAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTg                                                                                                                                        >  1:3178147‑M1/90‑150 (MQ=255)
               atgagcgtcgcagcgctgggtgtggcggcctgtgtgatgcggctgttcatgggggcggcgggcctgaccacacactgataaagagcggcAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTg                                                                                                                                        >  1:40327‑M1/90‑150 (MQ=255)
                tgagcgtcgcagcgctgggtgtggcggcctgtgtgatgcggctgttcatgggggcggcgggcctgaccacacactgataaagagcggcAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTgg                                                                                                                                       >  1:3039239‑M1/89‑150 (MQ=255)
                gttgttCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTgg                                                                                                                                       <  1:5901820/150‑1 (MQ=255)
                gttgttCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTgg                                                                                                                                       <  1:1313868/150‑1 (MQ=255)
                   gttCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCa                                                                                                                                    >  1:5254745/1‑150 (MQ=255)
                        ggggggTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGccc                                                                                                                              <  1:1116772/150‑1 (MQ=255)
                           gggTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCt                                                                                                                           <  1:999567/150‑1 (MQ=255)
                              tGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGcc                                                                                                                        >  2:2124264/1‑150 (MQ=255)
                                    aGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCtt                                                                                                                  >  1:4141516/1‑150 (MQ=255)
                                       tGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTAt                                                                                                               >  1:2476700/1‑150 (MQ=255)
                                        ggcggcAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATc                                                                                                              >  2:3957921/1‑150 (MQ=255)
                                         gcggcACCGCCCATGGTGGTCTTTCTGAACGCCCGAATGTCGGCGCGCCACCGCCCCGGTAGGTGGGGGCTTGCCCGCGAACACCGGCAACGCCTGGGCCATCAAGCGGGTTGCCTGCTTGGTGGGCAAGCCCGCTGCCGCTCTTTATCa                                                                                                             <  2:6077970/150‑1 (MQ=255)
                                         gcgcaaaCGCCCACTTTGGTATTTCTCAAAGCCCGATCTTTGGCGCGCCACCCCCCCGGTTGGTGGGGGCTGCCCCGCTAACCCGGGCACGGCCTGGCCCATCTCGCGGGTGGCGTGCTGCGGGGGAACGCCCGCTGCCGCTCTTTATCa                                                                                                             <  2:3107004/145‑1 (MQ=255)
                                           ggcACCGCCCTTGTGCGTCTTTCGGAAAGCCCGCAAGTTGGCGCGCCACCGCCCCCGTCGGTGGGGGTTTGCCCGCTAAAACCGGAACCGCCTGTGCAATATAGCGCTTGGCGTGCTTCTTGGGCAAGCCCGCTGCCGCTCTTTATCAgt                                                                                                           <  2:3788833/150‑1 (MQ=255)
                                           ggaaaCGCCCACGTTCTTCATTCTCAAAGCCCGCATGTCGGCGCGCCACCGCCCCCTTAGGCGCGGGCTGGCCCGCCAACACCGCCAACGCCTGTGCCATCTAGCGGGTTGCCTGCTTGGGGGCAACGCCCGCTGCCGCTCTTTATCAgt                                                                                                           <  2:3788625/147‑1 (MQ=255)
                                            gcAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAgtg                                                                                                          >  1:2055140/1‑150 (MQ=255)
                                                   ccATGTTCTTCATTCTGAAAGCCCGCACGTAGGCGCGCCACCGCCCCGGTCGGTGCGGGCTTGCCCGCGAACACCGGCACCGCCTGTGCCATCTCGCGGGTTGCGTGGTTGGTGGGAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTc                                                                                                   <  2:3191076/150‑1 (MQ=255)
                                                   ccAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTc                                                                                                   >  1:1676919/1‑150 (MQ=255)
                                                   ccAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTc                                                                                                   >  2:6395256/1‑150 (MQ=255)
                                                    cAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTGGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCAACTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCCCTGCCGCTCGTTATAAGTGTGTGGTCa                                                                                                  >  1:1758138/1‑150 (MQ=38)
                                                     aaGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAg                                                                                                 >  1:2826473/1‑150 (MQ=255)
                                                     aaGTTCGTATTTCTGAAAGCCCGCACGTCGGCGCGCCACCGCCCCCGTAGGAGGGGGCTGGCCCGCTAACACCGGCACCGCTTGTGCCATCTCGGGGGTTGCCTGCTGGGGGGGAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAg                                                                                                 <  2:2959408/150‑1 (MQ=255)
                                                     aaGTGCGTCTTTCTGAACGCCCGCACGTCGGCGCGCCCCCGCCCCCGTCGGCGCGGGGTGGCCCGCGACCACCGGCAACGCCTGTGCCATCTCGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAg                                                                                                 <  2:6728176/150‑1 (MQ=255)
                                                      tgTTCGTCATTCTGAAAGCCCGCATGTCGGGGCGCCACCGCCCCCGTAGGAGGGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCAATCTCGCGGGTTGCCTGATTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAgg                                                                                                <  2:5482342/149‑1 (MQ=255)
                                                         tCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGccc                                                                                             >  2:4890628/1‑150 (MQ=255)
                                                           gTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCCGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCccgc                                                                                           >  1:1199784/1‑150 (MQ=255)
                                                             cATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCccgccg                                                                                         >  1:6584391/1‑150 (MQ=255)
                                                               ttCTGAACGCCCGCACGTCGGCGCGCCACCGCCCCCTTTGGTGGGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTCGCGGGTTGCCTGCTTGGTGGGCATGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCgccgcc                                                                                       <  2:1876970/150‑1 (MQ=255)
                                                                tCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCAGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGccc                                                                                      >  2:1154450/1‑150 (MQ=255)
                                                                  tGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGccccc                                                                                    >  1:5010673/1‑150 (MQ=255)
                                                                   gAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCa                                                                                   >  1:93701/1‑150 (MQ=255)
                                                                    aacgCCCGCACGGCGGCCCCCCCCCCCCCCGGTGGGGGGGGGCTGCCCCGCCAACCCCGCCAACGCGTGTGCCATCGCGCGCGGTGCCTTCTTGGTGGCCACCCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCAt                                                                                  <  1:6302467/147‑1 (MQ=255)
                                                                         ccGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGc                                                                                                                                                   >  1:3322652/1‑82 (MQ=255)
                                                                           gCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCGGTGCCAACTAGCGCGGTGCCTGCTTCGTGGGCAAGCCCGCGGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCAGGAACAGc                                                                           >  1:3101692/1‑150 (MQ=38)
                                                                                tAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCAt                                                                      >  2:75699/1‑150 (MQ=255)
                                                                                tAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCAt                                                                      >  2:3488461/1‑150 (MQ=255)
                                                                                   gcgcgcCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATcac                                                                   >  2:4072142/1‑150 (MQ=255)
                                                                                   gcgcgcCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATcac                                                                   >  2:6852788/1‑150 (MQ=255)
                                                                                     gcgcCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCAGCTAGCGCGTTGCCTGCGGCGTGGGCAAGCCCGCTGCCGCTCGGGATCAGTGTGTGGTCAGGCCCGCCGCCCCCCTGAACCGCCGCATcacac                                                                 >  2:5173257/1‑150 (MQ=38)
                                                                                         ccccGCCCCGTTCGGTGGGGGTTTCCCCGCTAACACGGCAACCGCGTGTGCCTTCTCGCGCGTTGCCTGTTTGGTGGGAACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGgc                                                             <  1:5545610/148‑1 (MQ=255)
                                                                                         cACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGgc                                                             >  1:6673715/1‑150 (MQ=255)
                                                                                              cccccGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCca                                                        >  2:2564133/1‑150 (MQ=255)
                                                                                                 ccGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCacac                                                     >  2:2330696/1‑150 (MQ=255)
                                                                                                   gTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACAccc                                                   >  1:6700554/1‑150 (MQ=255)
                                                                                                          gGGGCTGGCCCGCTAACACCGGAAACGCCTGTGCAATCTCGCGCTTTGCCTGCTTCGTGGGAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTg                                            >  2:5837226‑M1/2‑150 (MQ=255)
                                                                                                           gggTTTGCCCGTGAACAGGGGCATCGCCTGTGCATTGTTGCGGTTTGCCTGGTTGGTGGGCTTGCCGGCTGCCTCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCGCATGATCAGCCGCATCACACAGGCCGCCTCACCCAGCGCTGc                                           <  2:6469439/150‑1 (MQ=255)
                                                                                                           GGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGc                                           <  1:6597540‑M1/150‑1 (MQ=255)
                                                                                                             GCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGa                                         <  1:4542933‑M1/150‑1 (MQ=255)
                                                                                                             GCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGa                                         <  1:4291901‑M1/150‑1 (MQ=255)
                                                                                                             GCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGa                                         <  1:5003407‑M1/150‑1 (MQ=255)
                                                                                                                  CCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTc                                    <  2:6636319‑M1/150‑1 (MQ=255)
                                                                                                                   CCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCa                                   <  1:1955883‑M1/150‑1 (MQ=255)
                                                                                                                        cacccccgCAACGGCTGGTGCAATCTCGCGGGTTGCCTGGTTGGTGGCAACCCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGc                               <  1:968673/143‑1 (MQ=255)
                                                                                                                        aaCACGGGAAACGCCTGTGCAATCTCGCGGGTTGCCTGATTCGTGGGAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGcc                              <  2:6537877/150‑1 (MQ=255)
                                                                                                                          cacCGGCACCGCCTGTGCCATCTAGGGGGTTGCCTGCTTGGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCa                            <  2:3697292/150‑1 (MQ=255)
                                                                                                                          cacCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGTGTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGccc                             >  1:1295947/1‑150 (MQ=255)
                                                                                                                           acCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAg                           >  2:845099/1‑150 (MQ=255)
                                                                                                                             cgccAACGCTTTTCCAATCTTGCGCTTTGCTTGCTTGGTGGGAACCCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGAt                         <  1:1581009/147‑1 (MQ=255)
                                                                                                                             cGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGAt                         >  1:4238891/1‑150 (MQ=255)
                                                                                                                             cGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGAt                         >  2:5641273/1‑150 (MQ=255)
                                                                                                                               gCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCa                       >  2:3367351/1‑150 (MQ=255)
                                                                                                                                cAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAg                      >  2:1768327/1‑150 (MQ=255)
                                                                                                                                 caccGCCTGTGCCATCTCGCGGGTGGCCTGGTGGGTGGGAACCCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAg                      <  1:750407/147‑1 (MQ=255)
                                                                                                                                 aaCGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGc                     <  2:4532692/150‑1 (MQ=255)
                                                                                                                                     gCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAAc                  >  2:3268395/1‑150 (MQ=255)
                                                                                                                                      cgtgtgCCATCTCGCGGGTTGCGTGCTTGGTGGGAACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACa                 <  1:912098/148‑1 (MQ=255)
                                                                                                                                      ccTGGGCCTTCTAGGGCGTTGCCTGCTTTTTGGGAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACa                 <  1:2753452/150‑1 (MQ=255)
                                                                                                                                               aTCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCaa            >  2:6682637/1‑147 (MQ=255)
                                                                                                                                               aTCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAAccc         >  1:6509444/1‑150 (MQ=255)
                                                                                                                                                  tAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGt      >  1:1233588/1‑150 (MQ=255)
                                                                                                                                                  tAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGt      >  2:2261316/1‑150 (MQ=255)
                                                                                                                                                    gcgcGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCg    >  2:2973818/1‑150 (MQ=255)
                                                                                                                                                      gcgTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCt  <  1:4444352/150‑1 (MQ=255)
                                                                                                                                                      gcgTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCt  >  1:297864/1‑150 (MQ=255)
                                                                                                                                                      gcgTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCt  <  1:4832405/150‑1 (MQ=255)
                                                                                                                                                      gcgTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCt  <  1:1002121/150‑1 (MQ=255)
                                                                                                                                                      gcgTTGCCTGCTCCGTGGGCACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCt  <  2:1970119/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
ACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCT  >  NC_002947/5278893‑5279190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: