Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 557,335 Δ1 bp 14.3% 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.14.3% 96.4 / 12.3 35intergenic (+75/‑134)rpsN/rpsH30S ribosomal protein S14/30S ribosomal protein S8
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base . (3/2);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

TCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAG  >  NC_002947/557209‑557480
                                                                                                                              |                                                                                                                                                  
tCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGCCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGTCCCTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTg                                                                                                                            >  1:250750/1‑150 (MQ=255)
                      gTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGc                                                                                                      >  1:270100/1‑150 (MQ=255)
                         ccAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCtt                                                                                                                        <  2:70814/129‑1 (MQ=255)
                         ccAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCtt                                                                                                                        >  1:70814/1‑129 (MQ=255)
                          cAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAc                                                                                                  <  2:262863/150‑1 (MQ=255)
                           aGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAcc                                                                                                 >  2:76857/1‑150 (MQ=255)
                             gTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTg                                                                                               >  2:221042/1‑150 (MQ=255)
                                    cAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg                                                                                        <  1:221042/150‑1 (MQ=255)
                                    cAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg                                                                                        <  1:399080/150‑1 (MQ=255)
                                    cAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGAAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg                                                                                        <  1:231760/150‑1 (MQ=255)
                                     aaGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTggg                                                                                       >  1:401786/1‑150 (MQ=255)
                                       ggCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCt                                                                                                                             <  2:663946/110‑1 (MQ=255)
                                       ggCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCt                                                                                                                             >  1:663946/1‑110 (MQ=255)
                                          ctcctGGTAATCGTTGCAGGCTAGGGGAGGGTGGGAGGACAGAGATCTTCCGTCGTGCGGTGACTTTTCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                  <  2:250750/150‑1 (MQ=255)
                                             ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGa                                                                                              >  2:225422/1‑135 (MQ=255)
                                             ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGa                                                                                              <  1:225422/135‑1 (MQ=255)
                                                gTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCaa                                                                            >  1:523598/1‑150 (MQ=255)
                                                             gCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGc                                                                                >  1:109775/1‑133 (MQ=255)
                                                             gCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGc                                                                                <  2:109775/133‑1 (MQ=255)
                                                              cAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGggtggt                                                              <  2:270100/150‑1 (MQ=255)
                                                               aaGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGtt                                                             >  2:280280/1‑150 (MQ=255)
                                                                           gAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCa                                                 >  2:598068/1‑150 (MQ=255)
                                                                                  aGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaa                                          <  1:410780/150‑1 (MQ=255)
                                                                                                    ggTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAg                        <  1:280280/150‑1 (MQ=255)
                                                                                                           ttGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt                                   >  1:798143/1‑132 (MQ=255)
                                                                                                           ttGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt                                   <  2:798143/132‑1 (MQ=255)
                                                                                                              aaacgaaacaaGCCCCATATAGGGGCTTGTTTGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGc                                                     >  1:628000/6‑110 (MQ=37)
                                                                                                              aaacgaaacaaGCCCCATATAGGGGCTTGTTTGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGc                                                     <  2:628000/105‑1 (MQ=37)
                                                                                                               caacaaACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCa             <  1:124426/150‑1 (MQ=255)
                                                                                                                 acaaacaaGCCCCTATATGGGGCTTGTATCGTGTCTGCCTTATGTCTAGAATGACAGTCTCACCTTAGCCTTGGTTTTGGCCCAAACGGATCGGGTGGTTGTCATAGCCGCACGGCTAATAATTCTTGTATCAGGAGCATCTATCCCATg           >  2:13991/1‑150 (MQ=255)
                                                                                                                  caaacaaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa          <  1:76857/150‑1 (MQ=255)
                                                                                                                   aaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa          >  2:422081/1‑149 (MQ=255)
                                                                                                                   aaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa          <  1:422081/149‑1 (MQ=255)
                                                                                                                     acaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTa       >  2:594396/1‑150 (MQ=255)
                                                                                                                          ccccTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAg  >  1:651891/1‑150 (MQ=255)
                                                                                                                              |                                                                                                                                                  
TCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAG  >  NC_002947/557209‑557480

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