Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 9,533 G→A 9.3% intergenic (‑587/‑9) rpmH ← / → dnaA 50S ribosomal protein L34/chromosomal replication initiator protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029479,5330GA 9.3% 117.5 / 2.8 43intergenic (‑587/‑9)rpmH/dnaA50S ribosomal protein L34/chromosomal replication initiator protein
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base A (2/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

TTCACAGAGGCGAAAAGCCGTTTTGGTCGATAAATGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTGCGCGTCTATGCGCCTAACCGTTTCGT  >  NC_002947/9396‑9681
                                                                                                                                         |                                                                                                                                                    
ttCACAGAGGCGAAAAGCCGTTTTGGTCGATAAATGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCgt                                                                                                                                            <  4:129493/148‑1 (MQ=255)
   acaGAGGCGAAAAGCCGTTTTGGTCGATAAATGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCgtgt                                                                                                                                          >  2:88596/1‑147 (MQ=255)
    caGAGGCGAAAAGCCGTTTTGGTCGATAAATGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAg                                                                                                                                       >  1:173734/1‑149 (MQ=255)
      gagGCGAAAAGCCGTTTTGGTCGATAAATGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTg                                                                                                                                     <  6:207857/149‑1 (MQ=255)
         gCGAAAAGCCGTTTTGGTCGATAAATGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGaa                                                                                                                                  >  2:225223/1‑149 (MQ=255)
               aGCCGTTTTGGTCGATAAATGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTgg                                                                                                                            <  7:28594/149‑1 (MQ=255)
               aGCCGTTTTGGTCGATAAATGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGAGGATATCCGTGTCAGTGGAACTTTgg                                                                                                                            >  5:196619/1‑149 (MQ=255)
                                aaTGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGc                                                                                                            <  1:225223/148‑1 (MQ=255)
                                 aTGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCtt                                                                                                          <  1:259001/149‑1 (MQ=255)
                                 aTGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCtt                                                                                                          <  2:173734/149‑1 (MQ=255)
                                      tGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCgt                                                                                                                                            >  5:181350/1‑110 (MQ=255)
                                      tGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCgt                                                                                                                                            <  6:181350/110‑1 (MQ=255)
                                          ttGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGa                                                                                                  <  1:88596/148‑1 (MQ=255)
                                               gTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATg                                                                                                <  5:275419/145‑1 (MQ=255)
                                               gTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATg                                                                                                >  6:275419/1‑145 (MQ=255)
                                                tGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAActg                                                                                           >  4:105184/1‑149 (MQ=255)
                                                   ttCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGAGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAActgcct                                                                                        <  6:196619/149‑1 (MQ=255)
                                                         aCGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCa                                                                                   <  7:209186/148‑1 (MQ=255)
                                                         aCGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGAGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAg                                                                                  >  7:62002/1‑149 (MQ=255)
                                                                cacaTGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTc                                                                            <  3:105184/148‑1 (MQ=255)
                                                                cacaTGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTc                                                                            >  1:171572/1‑148 (MQ=255)
                                                                        gATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTgg                                                                   <  2:12077/149‑1 (MQ=255)
                                                                          tAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGAt                                                                 <  7:26345/149‑1 (MQ=255)
                                                                           aaCTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTgg                                                                                                               <  4:255076/102‑1 (MQ=255)
                                                                           aaCTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTgg                                                                                                               >  3:255076/1‑102 (MQ=255)
                                                                              tGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGAtccgt                                                             >  4:112072/1‑149 (MQ=255)
                                                                                    cTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAActgcct                                                                                        >  4:239723/1‑116 (MQ=255)
                                                                                    cTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAActgcct                                                                                        <  3:239723/116‑1 (MQ=255)
                                                                                     tCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGc                                                         <  4:139687/146‑1 (MQ=255)
                                                                                     tCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGc                                                         >  3:139687/1‑146 (MQ=255)
                                                                                          cGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCaa                                                                          >  4:238331/1‑124 (MQ=255)
                                                                                          cGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCaa                                                                          <  3:238331/124‑1 (MQ=255)
                                                                                             tACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGaa                                              >  1:262189/1‑149 (MQ=255)
                                                                                                          tgtttTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGcgac                                  <  2:262189/148‑1 (MQ=255)
                                                                                                                  ccTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTgcgc                         <  4:15083/149‑1 (MQ=255)
                                                                                                                  ccTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTgcgc                         >  3:15083/1‑149 (MQ=255)
                                                                                                                   cTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGACGCCGAAGGCGACGAGTTGcgcg                        >  2:206572/1‑149 (MQ=255)
                                                                                                                   cTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTGcgcg                        >  7:236220/1‑149 (MQ=255)
                                                                                                                    tCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGc                                            >  2:77969/1‑128 (MQ=255)
                                                                                                                      aTCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTGCGCGTc                      >  8:233011/1‑148 (MQ=255)
                                                                                                                            cTTTCAAACTCAGAGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTGCGCGTCTATgcgc               <  8:62002/149‑1 (MQ=255)
                                                                                                                             tttCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTGCGCGTCTATGCGcc              >  5:69561/1‑149 (MQ=255)
                                                                                                                                       aGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTGCGCGTCTATGCGCCTAACCGTTTc    >  1:88765/1‑149 (MQ=255)
                                                                                                                                         gggATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGaa                                                                                              <  2:52100/57‑1 (MQ=255)
                                                                                                                                         gggATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGaa                                                                                              >  1:52100/1‑57 (MQ=255)
                                                                                                                                         gggATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGc                                     >  7:109103/1‑114 (MQ=255)
                                                                                                                                         gggATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGc                                     <  8:109103/114‑1 (MQ=255)
                                                                                                                                         gggATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTGCGCGTCTATGCGCCTAACCGTTTCGt  >  6:238160/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
TTCACAGAGGCGAAAAGCCGTTTTGGTCGATAAATGGCTGTTTTGTCGTGGTTCCTAACGTGTCCACATGTGGATAACTGAACGCTCGACCGGTACAATGGCGGTTTGTTTTTGCCTCATCCGGCTTTCAAACTCAGGGGATATCCGTGTCAGTGGAACTTTGGCAGCAGTGCGTGGAGCTTCTGCGCGATGAACTGCCTGCCCAGCAATTCAACACCTGGATCCGTCCGCTACAGGTCGAAGCCGAAGGCGACGAGTTGCGCGTCTATGCGCCTAACCGTTTCGT  >  NC_002947/9396‑9681

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