Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,687 T→C 13.6% intergenic (‑301/+26) sufA ← / ← rydB Fe‑S cluster assembly protein/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,6870TC13.6% 110.3 / 10.9 45intergenic (‑301/+26)sufA/rydBFe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base T (25/13);  new base C (3/3);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

CTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAA  >  NC_000913/1764541‑1764835
                                                                                                                                                  |                                                                                                                                                     
cTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCt                                                                                                                                                     <  4:327480/149‑1 (MQ=255)
  aaGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCtgt                                                                                                                                                   <  3:558097/149‑1 (MQ=255)
             aTCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGc                                                                                                                                       >  3:261469/1‑149 (MQ=255)
                 cATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTaa                                                                                                                                   <  2:390746/149‑1 (MQ=255)
                 cATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCG‑CACAGACTTaa                                                                                                                                   <  1:173422/149‑1 (MQ=255)
                     cttAGATAATAACCATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGt                                                                                                                               >  1:689588/2‑149 (MQ=255)
                          aTAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGt                                                                                                                               >  4:521317/1‑144 (MQ=255)
                          aTAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGGTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGt                                                                                                                               <  3:521317/144‑1 (MQ=255)
                           tAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTAc                                                                                                                         <  4:261469/149‑1 (MQ=255)
                                 ttattaTCTAAGGAGCATTAAGGCTGTTATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAGATAAGATAACCCTTAAGTCTGTG‑CGACACAGGCTTaa                                                                                                                                   <  1:23571/133‑1 (MQ=255)
                                 ttattaTCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTaa                                                                                                                                   >  2:23571/1‑133 (MQ=255)
                                       tCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGgcgc                                                                                                             >  4:682990/1‑149 (MQ=255)
                                       tCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGgcgc                                                                                                             >  4:514799/1‑149 (MQ=255)
                                       tCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTaa                                                                                                                                   <  4:461700/127‑1 (MQ=39)
                                       tCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTG‑TCGCACAGACTTaa                                                                                                                                   >  3:461700/1‑127 (MQ=39)
                                            ggAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTAtc                                                                                                        >  3:225025/1‑149 (MQ=255)
                                                 aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTATGGGTTTCTACCCCATCCGGCGCTTATCTCCgg                                                                                                   >  4:564692/1‑149 (MQ=255)
                                                  ttAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGc                                                                                                  >  4:603593/1‑149 (MQ=255)
                                                                aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCtt                                                                                    >  4:205931/1‑149 (MQ=255)
                                                                            ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAgg                                                                        >  3:185785/1‑149 (MQ=255)
                                                                                  tGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGGCTTTG‑CGACACAGGCTTAAGGGTTTCTACCCCCACCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGAt                                                                  <  3:682990/149‑1 (MQ=255)
                                                                                    aaaaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATaa                                                                <  2:689588/149‑1 (MQ=255)
                                                                                           cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAACAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGGTCTTGAAGGGGCGATAAGaataat                                                         <  3:651224/149‑1 (MQ=255)
                                                                                            tATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtc                                                        <  3:414966/149‑1 (MQ=255)
                                                                                              tttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAAtctc                                                      <  3:514799/149‑1 (MQ=255)
                                                                                                   gagaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATaaa                                                 >  3:154278/1‑149 (MQ=255)
                                                                                                           tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAccc                                         >  1:306549/1‑149 (MQ=255)
                                                                                                             aaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgc                                       >  1:287381/1‑149 (MQ=255)
                                                                                                              aaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc                                      >  4:149974/1‑149 (MQ=255)
                                                                                                                     cAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTT‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                               <  3:564692/149‑1 (MQ=255)
                                                                                                                          ttCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaa                          >  3:610047/1‑149 (MQ=255)
                                                                                                                          ttCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaa                          >  2:367534/1‑149 (MQ=255)
                                                                                                                             aGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACt                       >  3:172841/1‑149 (MQ=255)
                                                                                                                              gAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTg                      >  2:25766/1‑149 (MQ=255)
                                                                                                                              gAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTCTAACACAAACTg                      >  3:427084/1‑149 (MQ=255)
                                                                                                                               aaaTAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                     <  4:154278/149‑1 (MQ=255)
                                                                                                                                     gAAAACCCTTAAGTCTGTG‑TGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg               >  2:446908/1‑149 (MQ=255)
                                                                                                                                     gAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGGTTAg               >  1:669261/1‑149 (MQ=255)
                                                                                                                                     gAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg               >  2:27676/1‑149 (MQ=255)
                                                                                                                                        aaCCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtat            >  4:12117/1‑149 (MQ=255)
                                                                                                                                           ccTTAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGtattat         >  1:358279/1‑149 (MQ=255)
                                                                                                                                            cTTAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTa                <  4:305076/141‑1 (MQ=255)
                                                                                                                                            cTTAAGCCTGTG‑TCGCACAGACTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTa                >  3:305076/1‑141 (MQ=255)
                                                                                                                                             ttAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt       >  3:101795/1‑149 (MQ=255)
                                                                                                                                               aaGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttttt     <  4:610047/149‑1 (MQ=255)
                                                                                                                                                  tCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGGTTAGTATTATTTTTGaa  <  2:669261/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTG‑CGACACAGGCTTAAGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAA  >  NC_000913/1764541‑1764835

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑