Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 557,335 Δ1 bp 15.0% intergenic (+75/‑134) rpsN → / → rpsH 30S ribosomal protein S14/30S ribosomal protein S8

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947557,3350T.15.0% 110.9 / 15.9 40intergenic (+75/‑134)rpsN/rpsH30S ribosomal protein S14/30S ribosomal protein S8
Reads supporting (aligned to +/- strand):  ref base T (15/19);  new base . (3/3);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GCCTGGGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGAC  >  NC_002947/557187‑557483
                                                                                                                                                    |                                                                                                                                                     
gCCTGGGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAACAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCta                                                                                                                                                      <  2:1088963/150‑1 (MQ=255)
gCCTGGGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGAAGTAACAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCta                                                                                                                                                      >  2:839866/1‑150 (MQ=255)
   tGGGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTatat                                                                                                                                                   <  2:588288/150‑1 (MQ=255)
                gctgcGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtttcg                                                                                                                                     <  2:584033/150‑1 (MQ=255)
                                 tGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTa                                                                                                                    <  2:964041/150‑1 (MQ=255)
                                 ggggcgGTGGCGTACCAGGTATGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAAGCTTCCGACCGCCGGTGATCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTa                                                                                                                    <  2:186942/147‑1 (MQ=255)
                                    gcgGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGaa                                                                                                                 >  2:1059664/1‑150 (MQ=255)
                                                               cctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCATATAGGGGCTTgtttgttg                                                                                                                                   <  2:87114/105‑5 (MQ=38)
                                                               cctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCATATAGGGGCTTgtttgttg                                                                                                                                   >  1:87114/1‑101 (MQ=38)
                                                                ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTtgt                                                                                                                         >  2:1073354/1‑114 (MQ=255)
                                                                ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTtgt                                                                                                                         <  1:1073354/114‑1 (MQ=255)
                                                                   ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGcc                                                                                  >  2:1116504/1‑150 (MQ=255)
                                                                              tGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg                                                                                           <  2:642605/130‑1 (MQ=255)
                                                                              tGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg                                                                                           >  1:642605/1‑130 (MQ=255)
                                                                                  ggCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGggtg                                                                   <  1:223663/150‑1 (MQ=255)
                                                                                     aaGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGtt                                                                >  2:581443/1‑150 (MQ=255)
                                                                                       ggCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAGGTGAGccgg                                                                                            >  1:684965/1‑118 (MQ=255)
                                                                                       ggCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAGGTGAGccgg                                                                                            <  2:684965/120‑3 (MQ=255)
                                                                                         cgcgGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGTGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCa                                                            <  1:1059664/150‑1 (MQ=255)
                                                                                             ggTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGc                                                        <  1:1116504/150‑1 (MQ=255)
                                                                                              gTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGcc                                                       >  2:1182283/1‑150 (MQ=255)
                                                                                                     gggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtt                                                                                                                                        >  2:520849/1‑62 (MQ=255)
                                                                                                     gggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtt                                                                                                                                        <  1:520849/62‑1 (MQ=255)
                                                                                                     gggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGc                                                <  1:1049708/150‑1 (MQ=255)
                                                                                                       gAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTACCCAAACAAATCGGGTGGTTGTCATAGCCGCAAGGCta                                              >  2:561587/1‑150 (MQ=255)
                                                                                                                  cGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTTGTTTCCGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGt                                   <  1:1113641/150‑1 (MQ=255)
                                                                                                                  cGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTg                                    >  1:6066/1‑149 (MQ=255)
                                                                                                                  cGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTg                                    <  2:6066/149‑1 (MQ=255)
                                                                                                                               ccTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCt                      <  1:839866/150‑1 (MQ=255)
                                                                                                                                cTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAACTGAAGCTGGGTTTTTGCCCAAACGGTTAGAAGAGTGTGCAGAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTa                     >  1:570850/1‑150 (MQ=255)
                                                                                                                                      aacaaacaAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAg                                                  <  1:35421/115‑1 (MQ=255)
                                                                                                                                      aacaaacaAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAg                                                  >  2:35421/1‑115 (MQ=255)
                                                                                                                                        gaaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                                     >  2:419813/2‑126 (MQ=255)
                                                                                                                                        gaaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                                     <  1:419813/125‑1 (MQ=255)
                                                                                                                                         aaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                     <  1:773343/77‑1 (MQ=37)
                                                                                                                                         aaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                     >  2:773343/1‑77 (MQ=37)
                                                                                                                                               gCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGtt                                                                         >  2:884911/1‑83 (MQ=255)
                                                                                                                                               gCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGtt                                                                         <  1:884911/83‑1 (MQ=255)
                                                                                                                                               gCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGt                                                                          >  1:848806/1‑82 (MQ=255)
                                                                                                                                               gCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGt                                                                          <  2:848806/82‑1 (MQ=255)
                                                                                                                                                   cTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGAc  >  2:598277/1‑150 (MQ=255)
                                                                                                                                                   cTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGAc  <  1:581443/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCTGGGCCGTAACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGAC  >  NC_002947/557187‑557483

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