Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,174,189 C→T 43.7% intergenic (+20/+9) bfr → / ← gspO bacterioferritin, iron storage and detoxification protein/bifunctional prepilin leader peptidase and methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,174,1890CT43.7% 18.9 / 60.6 48intergenic (+20/+9)bfr/gspObacterioferritin, iron storage and detoxification protein/bifunctional prepilin leader peptidase and methylase
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base T (11/10);  total (24/24)
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

TGGGCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAC  >  W3110S.gb/4174123‑4174255
                                                                  |                                                                  
tGGGCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTc                                                                >  1:275100/1‑71 (MQ=255)
tGGGCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTc                                                                >  1:864948/1‑71 (MQ=255)
tGGGCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTc                                                                >  1:1461558/1‑71 (MQ=255)
        cAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGCGGGAGAACCCTCCCGGTtaa                                                        >  1:902015/1‑71 (MQ=255)
        cAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGtt                                                          <  1:1127853/69‑1 (MQ=255)
        cAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGtt                                                          <  1:382562/69‑1 (MQ=255)
        cAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGtt                                                          <  1:876117/69‑1 (MQ=255)
        cAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGt                                                           >  1:1074539/1‑68 (MQ=255)
        cAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTtaa                                                        >  1:58363/1‑71 (MQ=255)
        cAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTta                                                         >  1:401030/1‑70 (MQ=255)
         aaaaTTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGtt                                                          <  1:1039175/68‑1 (MQ=255)
               aTCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGATGGAGAATCCTCCCGGTTAATAagag                                                  <  1:1386226/70‑1 (MQ=255)
               aTCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAagag                                                  <  1:833440/70‑1 (MQ=255)
               aTCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAagag                                                  <  1:148610/70‑1 (MQ=255)
               aTCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAagag                                                  <  1:148976/70‑1 (MQ=255)
               aTCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAagag                                                  <  1:412833/70‑1 (MQ=255)
               aTCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAagag                                                  <  1:1669097/70‑1 (MQ=255)
                     aaGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATc                                               <  1:1951726/67‑1 (MQ=255)
                        cacaGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCa                                        <  1:1317020/71‑1 (MQ=255)
                        cacaGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCa                                        <  1:1313199/71‑1 (MQ=255)
                        cacaGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCa                                        <  1:1154484/71‑1 (MQ=255)
                        cacaGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCa                                        <  1:1086699/71‑1 (MQ=255)
                          caGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:984784/71‑1 (MQ=255)
                          caGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:844726/71‑1 (MQ=255)
                          caGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:399383/71‑1 (MQ=255)
                          caGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:1066414/71‑1 (MQ=255)
                                        agGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTtaata                                                      >  1:1472996/1‑41 (MQ=255)
                                        agGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        >  1:371374/1‑71 (MQ=255)
                                        agGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        >  1:111007/1‑71 (MQ=255)
                                        agGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        >  1:1374813/1‑71 (MQ=255)
                                        agGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        <  1:137217/71‑1 (MQ=255)
                                        agGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        <  1:1530623/71‑1 (MQ=255)
                                        agGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        <  1:1492546/71‑1 (MQ=255)
                                        agGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        <  1:140989/71‑1 (MQ=255)
                                        agGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        <  1:755416/71‑1 (MQ=255)
                                           ttGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCa                      >  1:1911158/1‑70 (MQ=255)
                                           ttGAGCCGTGTCGGGAGGGAGAACCCGCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAta                       >  1:764308/1‑68 (MQ=255)
                                            tGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCAtg                    >  1:683221/1‑71 (MQ=255)
                                            tGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCAt                     >  1:1268283/1‑70 (MQ=255)
                                             gAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATAcccc                                >  1:1540612/1‑58 (MQ=255)
                                                              agaATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAc  >  1:750586/1‑71 (MQ=255)
                                                              agaATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAc  >  1:799893/1‑71 (MQ=255)
                                                              agaATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAc  >  1:50290/1‑71 (MQ=255)
                                                              agaATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAc  >  1:278774/1‑71 (MQ=255)
                                                              agaACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATTATGCTGCCGAGCGCAAGCCAc  >  1:529581/1‑71 (MQ=255)
                                                              agaACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAc  >  1:351055/1‑71 (MQ=255)
                                                              agaACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAc  >  1:454046/1‑71 (MQ=255)
                                                              agaACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAc  >  1:1853748/1‑71 (MQ=255)
                                                                  |                                                                  
TGGGCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAC  >  W3110S.gb/4174123‑4174255

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