Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,540,239 T→A L26Q (CTG→CAG)  feoA → ferrous iron transporter, protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,540,2390TA94.3% 174.7 / 2.3 53L26Q (CTG→CAG) feoAferrous iron transporter, protein A
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base A (33/17);  total (35/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

ATGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCT  >  NC_000913/3540142‑3540335
                                                                                                 |                                                                                                
aTGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGc                                                                                                <  2:1397937/100‑1 (MQ=255)
   agaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCttt                                                                                             >  2:820359/1‑100 (MQ=255)
   agaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCttt                                                                                             >  2:805615/1‑100 (MQ=255)
    gaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTtc                                                                                            >  1:1522606/1‑100 (MQ=255)
    gaAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGAGAAATCAGCCCGGCATATCGCCAAAAACAGCTTtc                                                                                            >  1:163413/1‑100 (MQ=255)
      aaCAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTTGAAAATCACTGGCTTTTCCCATGAAAACAGCCCGGCATATCGCCAAAAACAGCTTtctc                                                                                          >  1:470647/1‑100 (MQ=255)
       aCAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTctct                                                                                         >  2:1461289/1‑100 (MQ=255)
        cAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCtt                                                                                        >  1:1472211/1‑100 (MQ=255)
           gTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGc                                                                                     >  2:298747/1‑100 (MQ=255)
            tAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCa                                                                                    >  1:1208560/1‑100 (MQ=255)
                  ccTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTAc                                                                              >  1:65269/1‑100 (MQ=255)
                  ccTAGGAAATACACTCCAGATACTGCGTGGAAAATCAATGTCTTTTCCCGTGAATTCAGCTCGAAAAAAAAAAAAAAACAGCTTTCTCTTGGCATGTTAc                                                                              <  2:223348/100‑1 (MQ=255)
                    tATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCt                                                                            >  1:411768/1‑100 (MQ=255)
                      tGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTgg                                                                          >  1:1037849/1‑100 (MQ=255)
                       gCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGc                                                                         <  2:1037842/100‑1 (MQ=255)
                       gCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGc                                                                         >  1:73823/1‑100 (MQ=255)
                       gCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGc                                                                         <  1:593097/100‑1 (MQ=255)
                       gCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCAAAAAACAGCTTTCTCTTGGCATGTTACCTGGc                                                                         <  2:523631/100‑1 (MQ=255)
                              acTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCtt                                                                  <  1:1026062/100‑1 (MQ=255)
                                tCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCtttt                                                                <  1:361286/100‑1 (MQ=255)
                                 ccAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTa                                                               >  2:1250565/1‑100 (MQ=255)
                                 ccAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTa                                                               <  1:575038/100‑1 (MQ=255)
                                 ccAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTa                                                               >  1:1201010/1‑100 (MQ=255)
                                    gATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAAtg                                                            >  2:84730/1‑100 (MQ=255)
                                      tACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAGACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAAtgtg                                                          >  2:1006169/1‑100 (MQ=255)
                                           cGTGGAAAATCACTCGCTGTCCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGATATCTCTTGGCATGTTACCTGGCTCCTCTTTTAATctggt                                                        <  2:861285/97‑1 (MQ=255)
                                                aaaaTCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCg                                                <  1:1122743/100‑1 (MQ=255)
                                                   aTCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTc                                             >  2:1305319/1‑100 (MQ=255)
                                                    tCACTTGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTcc                                            <  1:1478905/100‑1 (MQ=255)
                                                      aCTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCAc                                          >  2:562396/1‑100 (MQ=255)
                                                         ggCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCg                                       >  1:457745/1‑100 (MQ=255)
                                                              ccgatctGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGAc                                  <  1:267610/94‑1 (MQ=255)
                                                                cccGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGAccc                                >  1:1394768/1‑100 (MQ=255)
                                                                 ccGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGAcccc                               >  1:397616/1‑100 (MQ=255)
                                                                  cGTGAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCGCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCa                              >  2:966718/1‑100 (MQ=255)
                                                                     gAAATCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGAcagatcgg                           >  2:267609/1‑93 (MQ=255)
                                                                       aaTCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCat                         <  1:879749/100‑1 (MQ=255)
                                                                        aTCAGCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCata                        <  1:1298853/100‑1 (MQ=255)
                                                                            gCCCGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGa                    >  1:51783/1‑100 (MQ=255)
                                                                               cGGCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAAc                 <  1:385722/100‑1 (MQ=255)
                                                                                 gCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAAccc               >  1:996409/1‑100 (MQ=255)
                                                                                 gCATATCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAAccc               <  1:558469/100‑1 (MQ=255)
                                                                                    tatCGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCcgtc            >  2:779590/1‑100 (MQ=255)
                                                                                       cGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCgtg         <  1:779594/100‑1 (MQ=255)
                                                                                       cGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCgtg         >  2:595864/1‑100 (MQ=255)
                                                                                       cGCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCgtg         >  1:423733/1‑100 (MQ=255)
                                                                                        gCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCgtgt        >  2:81187/1‑100 (MQ=255)
                                                                                        gCCAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCgtgt        >  1:634473/1‑100 (MQ=255)
                                                                                         ccAAAATCAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTGTAATGTGGTGCGCGTCGCTCCACTCGGCGTCCCCATTCATATCGAAACCCGTCGtgtg       <  2:98945/100‑1 (MQ=255)
                                                                                         ccAAAAACAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGtgtg       >  2:1212383/1‑100 (MQ=255)
                                                                                           aaaaaCTGCTTTCCCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGtgtg       >  2:1225166/1‑98 (MQ=255)
                                                                                            aaaaCAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGATGCGCGTCGCTCCACTCGGCGACCCCATTCATATCAAAACCCGTCGTGTGa      >  2:98766/1‑98 (MQ=255)
                                                                                              aaCAGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCt  >  2:1422117/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ATGAGAAACAAGTAGGCACCTATGCAATACACTCCAGATACTGCGTGGAAAATCACTGGCTTTTCCCGTGAAATCAGCCCGGCATATCGCCAAAAACTGCTTTCTCTTGGCATGTTACCTGGCTCCTCTTTTAATGTGGTGCGCGTCGCTCCACTCGGCGACCCCATTCATATCGAAACCCGTCGTGTGAGCCT  >  NC_000913/3540142‑3540335

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