Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,161,403 G→A A50V (GCC→GTC)  fhuE ← ferric coprogen/ferric rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,161,4030GA95.2% 49.9 / ‑6.0 21A50V (GCC→GTC) fhuEferric coprogen/ferric rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (7/13);  total (7/14)
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.95e-01

ATATAGTTATCGATCTGGAATCCGCGGGAATAATAAAGAGCACGATCGGAATCCGCCTGACTTTTGCTGATCCCCAGCGTGTTTTCCATCACTTCGCCCAGCGTTTGTAACTGCTGATCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGATATTGATTATCCCTGTTAAATTGTGTTGAAAGCATTTGAAATCTCTTGGTTTATGTTGTCGTCAGGCTCATGCCGGAGGAGGGTGATGAGCGATTATATTTAGTGTGCGAAT  >  NC_000913/1161160‑1161628
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
atatAGTTATCGGTCTGGAATCCGCGGGAATAATAAAGAGCACGATCGGAATCCGCCTGACTTTTGCTGATCCCCAGCGTGTTTTCCATCACTTCGCCCAGCGTTTGTAACTGCTGCTCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGCGGAAGACGTTACACTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGaaac                                                                                                                                                                                                                            >  2:227110/1‑249 (MQ=255)
           gATCTGGAATCCGCGGGAATAATAAAGAGCACGATCGGAATCCGCCTGACTTTTGCTGATCCCCAGCGTGTTTTCCATCACTTCGCCCAGCGTTTGTAACTGCTGATCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATc                                                                                                                                                                                                                  >  1:1161/1‑250 (MQ=255)
           gATCTGGAATACGCGGGCATAATAAAGAGCACGATCGGAATCCGCCTGACTTTTGCTGATCCCCAGCGTGTTTTCCATAACTTCGCCCAGCGTTTGTAACTGCTGATCTTCCATACGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGCGTCATCTGCATTTTGGTACCCGCAGAGGCAGACGTTACGCTGTAATCATTTTCGCAATCATCTGGCGCTGTGACTGAACACGCAACAATc                                                                                                                                                                                                                  >  2:350359/1‑250 (MQ=255)
                ggAATCCGCGGGAATAATAAAGAGCACGATCGGAATCCGCCTGACTTTTGCTGATCCCCAGCGTGTTTTCCATCACTTCGCCCAGCGTTTGTAACTGCTGATCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCAcc                                                                                                                                                                                                               <  1:227110/248‑1 (MQ=255)
                                                    ccgtcTGATTTTTGCTGATCCCCATCGGTTATTCCTTCACTTCACCCAGCGTTTGTAACGGCTGTTCTTCCATCCGCTGCTGGATAACAATAGTGACCGTCTGAGGAATATCTCTGTGAGTCAGCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGTCATCAGCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTGTAGTGGCTGGTGCAGCAAAAGCGGTAGAAgg                                                                                                                                                                         <  2:251154/246‑1 (MQ=255)
                                                                         ccAGCGTGTTTTCCATCACTTCGCCCAGCGTTTGTAACTGCTGATCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCaa                                                                                                                                                   <  1:268574/251‑1 (MQ=255)
                                                                         ccAGCGTGTTTTCCATCACTTCGCCCAGCGTTTGTAACTGCTGATCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCaa                                                                                                                                                   <  1:350359/251‑1 (MQ=255)
                                                                           aGCGTGTTTTCCATCCATTCGCCCAGCGTTTGTAACTGCTGATCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAcc                                                                                                                                                 <  1:208580/251‑1 (MQ=255)
                                                                                  tgttCCATTACTTCTCCCAGCGTGTTTAACGGCTGATCTTCCATCCGCGGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGGGTCAGTTGCATATGGGTACCCGCAGAGGTAGACGTGACGCTGTAATCAGTTTCGCCATCATCTGGAGCTGTGACGGAACCCGCAACAATAACCGTTGCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCTGTGTTATGCAACCTGCaa                                                                                                                                            <  2:154968/248‑1 (MQ=255)
                                                                                        aTCACTTCGCCCAGCGTTTGTAACTGCTGATCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGGAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTg                                                                                                                                      >  1:289451/1‑250 (MQ=255)
                                                                                          cACTTCGCCCATCTGTTGTAACTGCGGATCGTCCCGCCGCTGCTGGCGACCATTAGTGACCTACTGAGGAAGATCACGTTTAGTCATCTGCAGTTTGGTTCCCGCAGAGGTATACGTTTCGCTGTAATCATTTTCGCCATCATCTGTAGCTGTGACTGATCCCTCAACAATCACCGTTTCTTCAGTGGCTTGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACg                                                                                                                                   <  2:296005/251‑1 (MQ=255)
                                                                                                            tAACTGCTGATCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCGGCAAAAGCGGCAGAAGGAAATAGGGCCAGTGCTATGCAACCGGCGAGTAGTGACGGTTTGGTGATGGCTTgtt                                                                                                                 >  1:845/1‑249 (MQ=255)
                                                                                                                          ccATCCGCTGCTGGCTAACAATAGTGACCGTCTGATGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCGTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCa                                                                                                                                                              <  2:339949/192‑1 (MQ=255)
                                                                                                                          ccATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCa                                                                                                                                                              >  1:339949/1‑192 (MQ=255)
                                                                                                                                 ctgctgGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGGGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAAACGGCAAGTAGTGACGGTTTGGTGATGGCTTGATATTGATTATCCCTGTTAAAtt                                                                                            >  1:357888/1‑251 (MQ=255)
                                                                                                                                            cAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGATATTGATTATCCCTGTTAAATTGTGTTGAAAGc                                                                                 <  1:91670/251‑1 (MQ=255)
                                                                                                                                                                atatCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGATATTGATTATCCCTGTTAAATTGTGTTGAAAGCATTTGAAATCTCTTGGTTTa                                                             <  1:180484/251‑1 (MQ=255)
                                                                                                                                                                  atCACGTGGAGTCAGCTGTATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAAACATTTTCGCCATCATCTGGAGCTGTGACTGAACCCTCAACAATCACCGTTGCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGTCGGTTTGGTGATGGCTTGATATTGATTATCCCTGTTAAATTGTGTTGAAAGCATTTGAAATCTCTTGGTTTAtg                                                           <  2:289451/251‑1 (MQ=255)
                                                                                                                                                                                                             ggtACGCTGTTATCATTTTCGCCATCATCTTGAGCTGTGACTGAACCCTCAACAATCTCCGGTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTATTGACGGTTTGGTTATGGCTTGATATTGATGATCCCTGGTAAATTGTGTTGAAAGCATTTGAAATCTCTTGGGTTATGTTGTCGGCAGGCTCATTCCGGAGGAGGGTGATGATCGATtata                <  2:360083/249‑1 (MQ=255)
                                                                                                                                                                                                                        aTCATTTGCGACATCATCTGGAGCTTGGACTGACCCCTCAAGAATCACCGGTTCTTCAGTGGCGGGTGCAGCAAAAGCGACATAAGGTAATAGTGCCAGTGCTATGCAAGCGTCAAGTAGTGACGGTTTGGTTATGGCTTGATATTGATTTTCCCTGTTAAATTGTGTTGAAAGGATTTGAAAGCTCTTGTGTTATGGTGTCGTCAGGCTCATGCCGGAGGAGGGTGATGAGCGATTATATGTAGTGTGCg     <  2:359562/251‑1 (MQ=255)
                                                                                                                                                                                                                             tttCGCCATCATCTGGAGCTGTGACTTAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCATTGCTATGCAACCGGCAAGTAGTTACGTTTTGGTGATGGCTTGATATTGATTATCCCTGTTAAATTGTGTTGAAAGCATTTGAAATCTCTTGGTTTATGTTGTCGTCAGGCTCATGCCGGAGGATGGTGATGAGCGATTATATTTAGTGTGCGaat  <  1:12226/249‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                  
ATATAGTTATCGATCTGGAATCCGCGGGAATAATAAAGAGCACGATCGGAATCCGCCTGACTTTTGCTGATCCCCAGCGTGTTTTCCATCACTTCGCCCAGCGTTTGTAACTGCTGATCTTCCATCCGCTGCTGGCTAACAATAGTGACCGACTGAGGAATATCACGTTGAGTCATCTGCATTTTGGTACCCGCAGAGGTAGACGTTACGCTGTAATCATTTTCGCCATCATCTGGAGCTGTGGCTGAACCCTCAACAATCACCGTTTCTTCAGTGGCTGGTGCAGCAAAAGCGGCAGAAGGTAATAGTGCCAGTGCTATGCAACCGGCAAGTAGTGACGGTTTGGTGATGGCTTGATATTGATTATCCCTGTTAAATTGTGTTGAAAGCATTTGAAATCTCTTGGTTTATGTTGTCGTCAGGCTCATGCCGGAGGAGGGTGATGAGCGATTATATTTAGTGTGCGAAT  >  NC_000913/1161160‑1161628

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