Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,329,979 C→T S216S (AGC→AGT sohB → S49 peptidase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,329,9790CT98.0% 160.3 / NA 50S216S (AGC→AGTsohBS49 peptidase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/23);  minor base G (1/0);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAG  >  NC_000913/1329863‑1330109
                                                                                                                    |                                                                                                                                  
aGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGca                                                                                                             <  1:255831/140‑1 (MQ=255)
      ggTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGc                                                                                                       <  1:1230758/140‑1 (MQ=255)
        ttACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCg                                                                                                     <  1:885280/140‑1 (MQ=255)
        ttACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCg                                                                                                     <  2:717339/140‑1 (MQ=255)
        ttACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCg                                                                                                     <  1:734606/140‑1 (MQ=255)
         tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCgg                                                                                                    <  2:1392329/140‑1 (MQ=255)
          aCATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCggg                                                                                                   <  1:508341/140‑1 (MQ=255)
          aCATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCggg                                                                                                   <  2:147969/140‑1 (MQ=255)
                 ggCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGtata                                                                                            >  1:1283851/1‑140 (MQ=255)
                  gCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATaa                                                                                           >  1:1391013/1‑140 (MQ=255)
                     tgtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAgcg                                                                                        <  1:709008/140‑1 (MQ=255)
                     tgtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAgcg                                                                                        <  1:435603/140‑1 (MQ=255)
                     tgtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAgcg                                                                                        <  2:410544/140‑1 (MQ=255)
                        gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCAc                                                                                     <  1:1197765/140‑1 (MQ=255)
                        gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCAc                                                                                     <  1:500624/140‑1 (MQ=255)
                            cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTg                                                                                 >  2:1369138/1‑140 (MQ=255)
                            cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTg                                                                                 >  1:523951/1‑140 (MQ=255)
                            cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACCGCCGGGCAGTATAAGCGCACGcag                                                                                 >  1:288017/1‑138 (MQ=255)
                            cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTACACACTGCCGGGCAGTATAAGCGCACGCTg                                                                                 >  1:1472862/1‑140 (MQ=255)
                            cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGGAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTg                                                                                 >  1:314925/1‑140 (MQ=255)
                                 aaaaTTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGtt                                                                            >  2:566432/1‑140 (MQ=255)
                                           ccGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGaa                                                                  <  2:704457/140‑1 (MQ=255)
                                                   tttGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCga                                                          >  1:260023/1‑140 (MQ=255)
                                                   tttGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCga                                                          >  1:124379/1‑140 (MQ=255)
                                                           tgtgGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGc                                                  >  2:286599/1‑140 (MQ=255)
                                                           tgtgGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGc                                                  >  1:1147885/1‑140 (MQ=255)
                                                             tgGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCgg                                                >  2:98052/1‑140 (MQ=255)
                                                             tgGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCgg                                                >  2:583672/1‑140 (MQ=255)
                                                             tgGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCgg                                                >  2:899266/1‑140 (MQ=255)
                                                             tgGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCgg                                                >  1:1178218/1‑140 (MQ=255)
                                                             tgGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCgg                                                >  1:1356562/1‑140 (MQ=255)
                                                             tgGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCgg                                                >  1:261894/1‑140 (MQ=255)
                                                               ggTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGga                                              >  1:1098517/1‑140 (MQ=255)
                                                               ggTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGga                                              >  1:1360245/1‑140 (MQ=255)
                                                               ggTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGga                                              >  2:332506/1‑140 (MQ=255)
                                                                     aTTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAAtt                                        <  1:925031/140‑1 (MQ=255)
                                                                     aTTGGCGTGGTGGCGCAAATGACCAACTTTAACCGCTTCCTGAAAAGTAAAGGTATTGATAACGAACTGCACACCGCCGGGCAGGATAAGCGCACGCTGACGATGCTGGGCGAAAATACCGAAGAAAGGCGGGAGAAAtt                                        >  2:1442437/1‑140 (MQ=255)
                                                                         gCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCcgc                                    >  2:916371/1‑140 (MQ=255)
                                                                              gtggCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAaga                               >  2:1270680/1‑140 (MQ=255)
                                                                              gtggCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAaga                               >  1:801000/1‑140 (MQ=255)
                                                                               tggCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAagag                              <  1:534478/140‑1 (MQ=255)
                                                                                       aTGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGa                      <  1:419163/140‑1 (MQ=255)
                                                                                       aTGCCCAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGa                      <  2:953104/140‑1 (MQ=255)
                                                                                           ccAACTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACg                  <  1:286599/140‑1 (MQ=255)
                                                                                              aCTTTAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCAt               <  1:899266/140‑1 (MQ=255)
                                                                                                 ttAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAg            <  2:268794/140‑1 (MQ=255)
                                                                                                  tAACCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAgt           <  2:275156/140‑1 (MQ=255)
                                                                                                    aCCGCTTCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAgttg         >  2:201398/1‑140 (MQ=255)
                                                                                                          tCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTaaa   >  2:337524/1‑140 (MQ=255)
                                                                                                          tCCTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTaaa   >  2:1001408/1‑140 (MQ=255)
                                                                                                           ccTGAAAAGTAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAg  <  2:405247/140‑1 (MQ=255)
                                                                                                                    |                                                                                                                                  
AGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCTGAAAAGCAAAGATATTGATATCGAACTGCACACTGCCGGGCAGTATAAGCGCACGCTGACGTTGCTGGGTGAAAATACCGAAGAAGGGCGGGAGAAATTCCGCGAAGAGTTGAACGAAACGCATCAGTTGTTTAAAG  >  NC_000913/1329863‑1330109

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