Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,735,444 C→G 13.4% R374R (CGG→CGC) ‡ csdA ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,735,4440CG13.4% 77.7 / 6.3 30R374R (CGG→CGC) ‡csdAcysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base C (15/11);  new base G (2/2);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

ACGCTGTTCCCAGCCTTTGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGG  >  NC_002947/1735295‑1735593
                                                                                                                                                     |                                                                                                                                                     
aCGCTGTTCCCAGCCTTTGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATAcc                                                                                                                                                       <  1:203603/150‑1 (MQ=255)
aCGCTGTTCCCAGCCTTTGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGCCATCGCTGTCGTTGTACAAACCCAGCGATAcc                                                                                                                                                       <  2:776671/150‑1 (MQ=255)
  gCTGTTCCCAGCCTTTGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCg                                                                                                                                                     >  2:1138381/1‑150 (MQ=255)
  gCTGTTCCCAGACTTTGCATTATGGGAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCg                                                                                                                                                     >  2:1012773/1‑150 (MQ=255)
         ccAGCCTTTGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACGGGTATcgct                                                                                                                                              >  1:636152/1‑149 (MQ=255)
          cAGCCTTTGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACGGGTATcgctg                                                                                                                                             <  2:636152/150‑3 (MQ=255)
                ttGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTcca                                                                                                                                       <  1:122456/150‑1 (MQ=255)
                  gCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTccagc                                                                                                                                     <  2:203815/150‑1 (MQ=255)
                          gAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCcagccccacccc                                                                                                                             >  2:369564/1‑150 (MQ=255)
                                  tCCAGCGCTTGGTGCGCCCGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTcagca                                                                                                                     <  2:382719/150‑1 (MQ=255)
                                    cAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGt                                                                                                                   <  2:1190562/150‑1 (MQ=255)
                                     aGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGACACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTg                                                                                                                  >  2:163441/1‑150 (MQ=255)
                                          ttGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATg                                                                                                             >  1:945252/1‑150 (MQ=255)
                                           tGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATgg                                                                                                            <  1:543402/150‑1 (MQ=255)
                                                cgcCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACGGGTATcgctgg                                                                                                                                            >  2:764773/1‑110 (MQ=255)
                                                cgcCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACGGGTATcgctgg                                                                                                                                            <  1:764773/113‑4 (MQ=255)
                                                 gcCTGTTCCGGAAGGCTCATAGAAACAGCACCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAgt                                                                                                      >  2:12923/1‑150 (MQ=255)
                                                       tCCGGCAGGCTCATCGCAACAGCTCCAGACCTTGGTCCAGCGCGTCGAAGAAACGCTGCCCGTCATCGCTGTCGTGGTACCAACCCGGCGATACCCGTATCGCGCCCTCCAGCCCCCGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGAc                                                                                                <  1:133553/150‑1 (MQ=255)
                                                                 tCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGgcccgcac                                                                                      >  2:516163/1‑150 (MQ=255)
                                                                 tCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTAGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACAGgcacgcac                                                                                      >  1:338875/1‑150 (MQ=255)
                                                                              tCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGt                                                                         >  1:609019/1‑150 (MQ=255)
                                                                               ccAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCTTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGtt                                                                        <  1:932381/150‑1 (MQ=255)
                                                                                 aGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCCGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtca                                                                      <  1:1138381/150‑1 (MQ=255)
                                                                                   gcattGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtcagt                                                                    >  2:350817/4‑150 (MQ=255)
                                                                                    ccTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCCGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtcagtc                                                                   <  2:481494/150‑1 (MQ=255)
                                                                                     cTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTtcagtca                                                                  >  2:418859/1‑150 (MQ=255)
                                                                                                     gaagaaACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATc                                                  >  1:78158/1‑150 (MQ=255)
                                                                                                                           cTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATg                            <  1:418859/150‑1 (MQ=255)
                                                                                                                              tCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACg                         >  2:1028897/1‑150 (MQ=255)
                                                                                                                                          cccAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAgg             <  2:945252/150‑1 (MQ=255)
                                                                                                                                            cAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGc           >  2:268598/1‑150 (MQ=255)
                                                                                                                                              gCGATACCCGGAACGCGCCCCCAAGCCCCAGCACCTTCAGCAGTGGCATGGCGCCGAGGGGACCGGCACGCACGGCACTGCCCTATGCAGTCAGCAAATGAACAAAATCGGCGTAGTGCACGCACTCGATGACGAAACTGGGCAGTGCCg         >  2:382584/1‑150 (MQ=255)
                                                                                                                                                     ccGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCgg  <  1:163441/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGCTGTTCCCAGCCTTTGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGG  >  NC_002947/1735295‑1735593

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