Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,082,896 T→A R218R (CGT→CGA efeO → inactive ferrous ion transporter EfeUOB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,082,8960TA97.4% 118.6 / ‑5.0 38R218R (CGT→CGAefeOinactive ferrous ion transporter EfeUOB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (17/20);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GACGCGCCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGT  >  NC_000913/1082767‑1083024
                                                                                                                                 |                                                                                                                                
gACGCGCCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGaaaaa                                                                                                                         >  1:720453/1‑139 (MQ=255)
 aCGCGCCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAg                                                                                                                        >  1:675304/1‑139 (MQ=255)
   gcgcCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCa                                                                                                                      <  2:1171150/139‑1 (MQ=255)
    cgcCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGCATGGCAGCATTGACGCCCGTGCAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCAACTGGAAAAAGCAt                                                                                                                     >  2:838145/1‑139 (MQ=255)
        aGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGtt                                                                                                                 <  2:675304/139‑1 (MQ=255)
              aTGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGa                                                                                                           <  2:1124765/139‑1 (MQ=255)
               tGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGac                                                                                                          <  2:720453/139‑1 (MQ=255)
                          aaCCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCaa                                                                                               <  1:1062279/139‑1 (MQ=255)
                                 tGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGggatg                                                                                        <  2:1205527/139‑1 (MQ=255)
                                          gTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCGCCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGt                                                                                 >  2:726972/1‑137 (MQ=255)
                                          gTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGt                                                                                 <  1:726972/137‑1 (MQ=255)
                                                cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGag                                                                         <  1:194826/139‑1 (MQ=255)
                                                 gATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGagc                                                                        <  1:635919/139‑1 (MQ=255)
                                                   tCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGagcag                                                                      <  1:1138460/139‑1 (MQ=255)
                                                     tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAacac                                                                                                     >  2:970237/1‑106 (MQ=255)
                                                     tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAacac                                                                                                     <  1:970237/106‑1 (MQ=255)
                                                       gATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATc                                                                                    >  1:1244729/1‑121 (MQ=255)
                                                       gATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATc                                                                                    <  2:1244729/121‑1 (MQ=255)
                                                             agcaTTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCAACCGACTGGAAAAAGAATGGTTTGCCAACAACACCACCAAAGGGAGGGATCGGAAGGCTGAGCGGCTTAATACCg                                                            >  1:1164059/1‑139 (MQ=255)
                                                              gcaTTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGacaa                                                                                                        <  1:353495/94‑1 (MQ=255)
                                                              gcaTTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGacaa                                                                                                        >  2:353495/1‑94 (MQ=255)
                                                                                         aGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGaaa                                                                                                                           >  2:798969/1‑48 (MQ=255)
                                                                                         aGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGaaa                                                                                                                           <  1:798969/48‑1 (MQ=255)
                                                                                          gCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCa                                                                                   >  2:565246/1‑87 (MQ=255)
                                                                                          gCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCa                                                                                   <  1:565246/87‑1 (MQ=255)
                                                                                           cAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGc                                                                                                                       >  1:933876/1‑50 (MQ=255)
                                                                                           cAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGc                                                                                                                       <  2:933876/50‑1 (MQ=255)
                                                                                              aaaGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAAc                           >  1:583070/1‑139 (MQ=255)
                                                                                                    gccgACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCtt                     >  2:799676/1‑139 (MQ=255)
                                                                                                        aCCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTccc                 <  1:492615/139‑1 (MQ=255)
                                                                                                                  cACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCa                                                                                   >  1:912453/1‑63 (MQ=255)
                                                                                                                  cACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACTCCACCAAAGGGATGGATCa                                                                                   <  2:912453/63‑1 (MQ=255)
                                                                                                                   aCTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAg      >  2:1283360/1‑139 (MQ=255)
                                                                                                                    cTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCt                                                                    <  2:913325/76‑1 (MQ=255)
                                                                                                                    cTGGTTTCCACCGACAGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCt                                                                    >  1:913325/1‑76 (MQ=255)
                                                                                                                       gTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGgga                                                                                          <  2:573641/51‑1 (MQ=255)
                                                                                                                       gTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGgga                                                                                          >  1:573641/1‑51 (MQ=255)
                                                                                                                       gTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGt  >  2:1240854/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                
GACGCGCCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGT  >  NC_000913/1082767‑1083024

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