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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,735,4450CG13.4% 75.5 / 6.3 30R374P (CGG→CCG) ‡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 = 5.19e-01

GCTGTTCCCAGCCTTTGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGG  >  NC_002947/1735297‑1735593
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                            <  1:764773/113‑4 (MQ=255)
                                              cgcCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACGGGTATcgctgg                                                                                                                                            >  2:764773/1‑110 (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)
                                                                                                                                                    |                                                                                                                                                    
GCTGTTCCCAGCCTTTGCCCTGTTGAAACGCTTCCAGCGCTTGGTGCGCCTGTTCCGGCAGGCTCATCGCAACAGCTCCAGGCCTTGGTCCAGCGCGTCGAAGAAACGCTGCACGTCATCGCTGTCGTTGTACAAACCCAGCGATACCCGTATCGCGCCCTCCAGCCCCAGCCCCTTCAGCAGTGGCATGGCGCAGTGGTGACCGGCACGCACGGCAATGCCCTGTTCAGTCAGCAAATGAGCAATATCGGCGTTGTGCACGCCCTCGATGACGAAACTGGCCAGGGCCGTCTGCGG  >  NC_002947/1735297‑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: