Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,539,816 T→A 10.5% intergenic (+110/‑347) yhgF → / → feoA putative transcriptional accessory protein/ferrous iron transporter, protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,539,8160TA10.5% 154.9 / 6.7 47intergenic (+110/‑347)yhgF/feoAputative transcriptional accessory protein/ferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base A (2/3);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAC  >  NC_000913/3539694‑3539958
                                                                                                                          |                                                                                                                                               
cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGATAT‑TCTCCGGc                                                                                                                                       >  3:432748/1‑132 (MQ=255)
cAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGATAT‑TCTCCGGc                                                                                                                                       <  4:432748/132‑1 (MQ=255)
 aaaaaaCGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCa                                                                                                                     >  4:1432917/1‑149 (MQ=255)
         ttAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCata                                                                                                             <  3:1432917/149‑1 (MQ=255)
                            tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATAT‑TCTCCGGCActt                                                                                                                                   >  2:262236/1‑108 (MQ=38)
                            tACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAATGCCGGAGATAT‑TCTCCGGCActt                                                                                                                                   <  1:262236/108‑1 (MQ=38)
                                                         gCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaat                                                             >  3:592441/1‑149 (MQ=255)
                                                             aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatt                                                                                    >  1:1025970/1‑122 (MQ=255)
                                                             aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAtcatt                                                                                    <  2:1025970/122‑1 (MQ=255)
                                                             aTAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                         <  4:592441/149‑1 (MQ=255)
                                                                aGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAatata                                                                                                        <  3:866898/99‑1 (MQ=255)
                                                                aGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAatata                                                                                                        >  4:866898/1‑99 (MQ=255)
                                                                 ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAata                                                                                             >  3:1097975/1‑109 (MQ=255)
                                                                 ggCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAata                                                                                             <  4:1097975/109‑1 (MQ=255)
                                                                       cACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAg                                                   <  2:1343996/145‑1 (MQ=255)
                                                                       cACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAg                                                   >  1:1343996/1‑145 (MQ=255)
                                                                          gccgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaatt                                                            >  1:579693/1‑133 (MQ=255)
                                                                          gccgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCttaatt                                                            <  2:579693/133‑1 (MQ=255)
                                                                            cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAc                                                                   >  1:373401/1‑124 (MQ=255)
                                                                            cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAAc                                                                   <  2:373401/124‑1 (MQ=255)
                                                                            cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                         >  3:72224/1‑134 (MQ=255)
                                                                            cgcATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                         <  4:72224/134‑1 (MQ=255)
                                                                                  tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCAttttt                                                                                                                                 >  1:479627/1‑56 (MQ=255)
                                                                                  tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCAttttt                                                                                                                                 <  2:479627/56‑1 (MQ=255)
                                                                                  tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAt                                                                                        <  1:561769/97‑1 (MQ=255)
                                                                                  tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTAt                                                                                        >  2:561769/1‑97 (MQ=255)
                                                                                  tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaaa                                                                    <  4:1310970/117‑1 (MQ=255)
                                                                                  tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGaaa                                                                    >  3:1310970/1‑117 (MQ=255)
                                                                                  tGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAAtt                                    >  3:1395075/1‑149 (MQ=255)
                                                                                     cAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCAc                                 <  4:1395075/149‑1 (MQ=255)
                                                                                      aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCataa                                                                                                            >  4:1223470/1‑73 (MQ=255)
                                                                                      aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCataa                                                                                                            <  3:1223470/73‑1 (MQ=255)
                                                                                      aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGc                                                  >  2:1084388/1‑131 (MQ=255)
                                                                                      aacgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGc                                                  <  1:1084388/131‑1 (MQ=255)
                                                                                       acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAat                                                                                                                                                 >  3:1007352/1‑36 (MQ=255)
                                                                                       acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAat                                                                                                                                                 <  4:1007352/36‑1 (MQ=255)
                                                                                       acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                         <  2:1237108/123‑1 (MQ=255)
                                                                                       acgaacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                         >  1:1237108/1‑123 (MQ=255)
                                                                                          aacgCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTc                            <  4:178594/149‑1 (MQ=255)
                                                                                                  gCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa                                                    >  2:493098/1‑117 (MQ=255)
                                                                                                  gCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTa                                                    <  1:493098/117‑1 (MQ=255)
                                                                                                      gACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATAtttt                >  4:560889/1‑149 (MQ=255)
                                                                                                              aatGCCGGAGATAT‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                         <  2:224090/98‑1 (MQ=39)
                                                                                                              aatGCCGGAG‑ATATTCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTaaa                                                         >  1:224090/3‑100 (MQ=39)
                                                                                                                   cGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTa                                                                              >  3:774971/1‑74 (MQ=255)
                                                                                                                   cGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTa                                                                              <  4:774971/74‑1 (MQ=255)
                                                                                                                    ggAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAc  >  1:851499/1‑149 (MQ=255)
                                                                                                                        aaTA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCgg                                          >  3:937732/1‑105 (MQ=255)
                                                                                                                        aaTA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCgg                                          <  4:937732/105‑1 (MQ=255)
                                                                                                                         ata‑tCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTatat                   >  2:216284/1‑127 (MQ=255)
                                                                                                                         ata‑tCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTatat                   <  1:216284/127‑1 (MQ=255)
                                                                                                                          |                                                                                                                                               
CAAAAAACGTTAAACGCCCGTACTGGCCTACGGTTCGAATTTGCACGAAATCGTAGGGCAGATAAGGCGTTCACGCCGCATCTGGCAACGAACGCCTTGCCTGACATAAAAGTGCCGGAGAATA‑TCTCCGGCATTTTTATTCCACAGCCAAACTCATAATATATTCCGGCAATATTTATCATTTCATTAACAACTGAAACCTTAATTAAACATTAGCCAGTCCGGGTAATTCACTATTCGAATTATATTTTCGCTGCGATATAAC  >  NC_000913/3539694‑3539958

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