Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,173,819 G→A 17.4% V43M (GTG→ATG)  bfr → bacterioferritin, iron storage and detoxification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,173,8190GA17.4% 53.2 / 5.1 23V43M (GTG→ATG) bfrbacterioferritin, iron storage and detoxification protein
Reads supporting (aligned to +/- strand):  ref base G (10/9);  new base A (2/2);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GTCGCAATCAATCAGTACTTTCTCCATGCCCGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATATTGAGCGCATTCTTTTTCTGG  >  W3110S.gb/4173750‑4173888
                                                                     |                                                                     
gTCGCAATCAATCAGTACTTTCTCCATGCCCGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAatgatg                                                                       <  1:1744308/70‑1 (MQ=255)
  cGCAATCAATCAGTACTTTCTCCATGCCCGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGAtgtg                                                                     >  1:1075515/1‑70 (MQ=255)
               tACTTTCTCCATGCCCGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATc                                                       <  1:1530764/71‑1 (MQ=255)
                  tttCTCCATGCCCGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATCCa                                                     <  1:439108/70‑1 (MQ=255)
                    tctcCATGCCCGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATATGGAGTATCATg                                                           >  1:954189/1‑62 (MQ=255)
                        cATGCCCGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATCCATtgat                                                <  1:296584/69‑1 (MQ=255)
                         aTGCCCGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATCCATTgatgag                                             <  1:1586986/71‑1 (MQ=255)
                             ccGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGTAGTATCATGAATCCATTgatgagatga                                         <  1:1082098/71‑1 (MQ=255)
                               gAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATCCATTGATGAGATGaaa                                       >  1:935301/1‑71 (MQ=255)
                               gAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATCCATTGATGAGATGaaa                                       >  1:154600/1‑71 (MQ=255)
                                aaTGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATCCATTGATGAGATGaaa                                       >  1:57677/1‑70 (MQ=255)
                                 aTGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTAt                                                               >  1:1201329/1‑45 (MQ=255)
                                  tGTTTAAAAACTGGGGTCTCAAACGTCTCAATGAtgtg                                                                     >  1:1172253/1‑38 (MQ=255)
                                    tttAAAAACTGGTGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATCCATTGATGAGATGAAACACGc                                  >  1:437397/1‑71 (MQ=255)
                                    tttAAAAACTGGGGTCTCAAACGTCTCAATGATATGTAGTATCATGAATc                                                       <  1:1699781/50‑1 (MQ=255)
                                    tttAAAAACTGGGGTCTCAAACGTCTCAATGATATGGAGTATCATGAATc                                                       <  1:829948/50‑1 (MQ=255)
                                      tAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGa                                                                   >  1:1651891/1‑36 (MQ=255)
                                                  tctcAAACGTCTCAATGATGTGGAGTATCATGAATc                                                       <  1:1316007/36‑1 (MQ=255)
                                                  tctcAAACGTCTCAATGATATGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATAtt                     >  1:1469177/1‑70 (MQ=255)
                                                     cAAACGTCTCAATGATGTGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGt                           <  1:1014077/61‑1 (MQ=255)
                                                      aaaCGTCTCAATGATGTGGAGTATCATGAATCCATtgatg                                               >  1:761122/1‑40 (MQ=255)
                                                      aaaCGTCTCAATGATGTGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATATTGAgcg                >  1:613684/1‑71 (MQ=255)
                                                      aaaCGTCTCAATGATGTGGAATATCATGAATCCATTgatgaga                                            >  1:878873/1‑43 (MQ=255)
                                                             tcAATGATGTGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATATTGAGCGCATTCtt         >  1:868264/1‑71 (MQ=255)
                                                               aaTGATGTGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATATTGAGCGCATTCttt        <  1:301806/70‑1 (MQ=255)
                                                               aaTGATGTGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATATTGAGCGCATTCttt        <  1:119848/70‑1 (MQ=255)
                                                                atgatgTGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATATTGAGCGCATTCttttt      <  1:32250/71‑1 (MQ=255)
                                                                 tgatgTGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATAtt                     <  1:81713/55‑1 (MQ=255)
                                                                    tgtgGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATATTGAGCGCATTCTTTTTCTgg  <  1:1529117/71‑1 (MQ=255)
                                                                     gtgGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATATTGAGCGCATTCTTTTTCTgg  <  1:330535/70‑1 (MQ=255)
                                                                     |                                                                     
GTCGCAATCAATCAGTACTTTCTCCATGCCCGAATGTTTAAAAACTGGGGTCTCAAACGTCTCAATGATGTGGAGTATCATGAATCCATTGATGAGATGAAACACGCCGATCGTTATATTGAGCGCATTCTTTTTCTGG  >  W3110S.gb/4173750‑4173888

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