Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,8370AG44.4% 42.5 / 13.3 36I572T (ATC→ACC) fhuEferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (5/15);  new base G (16/0);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 3.22e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCC  >  NC_000913/1159700‑1159976
                                                                                                                                         |                                                                                                                                           
tATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTACCgg                                                                                                                                 >  1:555091/1‑150 (MQ=255)
        gCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTg                                                                                                                          <  1:495976/149‑1 (MQ=255)
         cAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGGCGGGGGTAGGGGTACCCGTggggtggg                                                                                                                        >  1:388664/1‑144 (MQ=255)
          aaTATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGGAGGTGTACCGGGGGGCTGAGCg                                                                                                                      >  1:7162/1‑151 (MQ=255)
          aaTATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCg                                                                                                                      <  1:439883/151‑1 (MQ=255)
             ataGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTTCTGCCGGGGGTAGGTGTACCGGGGGACTGAgaggc                                                                                                                    >  2:521449/1‑146 (MQ=255)
                gcgcGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTTCTGCCGGGGGTAGGGGTAACGGTGGGCTGAGCGACCtt                                                                                                                 >  1:621269/1‑150 (MQ=255)
                 cgcgTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGGGGACTGAGCGACATTAt                                                                                                               >  2:209462/1‑151 (MQ=255)
                    gTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCt                                                                                                            <  1:252008/151‑1 (MQ=255)
                          ccAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCaa                                                                                                      <  2:621280/151‑1 (MQ=255)
                                       gCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATg                                                                                           <  1:627554/149‑1 (MQ=255)
                                                   tAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTTCTGCCGGGGGTAGGGGGACCGGTGGACTTAGCGGCATTTTCCTGCTCCATACGGAAGATGGGTAACGgggggg                                                                              >  1:565362/1‑150 (MQ=255)
                                                    aaTTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTc                                                                            <  1:416431/151‑1 (MQ=255)
                                                           ccGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACgg                                                                      <  1:531913/150‑1 (MQ=255)
                                                              ttGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGGTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTCACGTGGGGGGCAGACGGCTAt                                                                  >  2:428293/1‑151 (MQ=255)
                                                                    tCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTTTCCTGCTTAATACGGAAGATGGCTAACGGGGGGGGCAGACGGCTATTCTTcc                                                            >  2:53259/1‑151 (MQ=255)
                                                                    tCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATc                                                             <  1:393614/150‑1 (MQ=255)
                                                                       aaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTACCGGGGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGGGGTGGGCAg                                                                          >  2:34980/1‑134 (MQ=255)
                                                                       aaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACg                                                                       <  1:34981/137‑1 (MQ=255)
                                                                           ccACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCCCCCTTGCTGCCCGGGGTAGGTGTACCGGTGGAATGAGCGACATTACCCTGCTTCATACGGGAGATGGCTAACGTGGTGGGGCGACGGCTAtttctcccggcc                                                       >  1:424888/1‑139 (MQ=255)
                                                                             aCCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGAtt                                                   <  2:392435/151‑1 (MQ=255)
                                                                                   ttACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCTCCGTTGCTGCCGGGGGGAGGGGTACCCGGGGAATGAGGGACATTTTTCTGCTCCATTCGGAAGATGGGTAACGTGGGGGTCAGACGGGTATTCATCCAGTCCGATTTTAGa                                              >  1:162318/1‑150 (MQ=255)
                                                                                     actgactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGAATGAGCGACATTTTCCTGCTCCATACGGGAGATGGCTTACGGGGGGGGCAGACGGGTATTCATCCAGTCCGATTTTAGAccc                                           >  1:235234/1‑151 (MQ=255)
                                                                                          ctgTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTATCCTTCTCAATAAGGAAGATGGGTAACGTGGGGGTCAGAAGGCTATTCATCCAGTCCGATTTCAGACCCAGCTc                                      >  1:578560/1‑151 (MQ=255)
                                                                                            gTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTTCTTCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGGGGTCAGACGGCTATTTATCCAGTCCGATTTCAGACCCCGCTc                                      >  2:30420/1‑149 (MQ=255)
                                                                                                cATCCACCGCTTTATAGGCGGTTTTGCCGTTTCTGCCGGGGGTAGGGGTACCGGGGGACTGAGCGACCTTTTCCTTCTCCAAACGGGAGAAGGCTAACGGGGGGGGCAGAAGGCTGTTCATCCCGTCCGGTTTTAGACCCAGCTCGGAgt                                 >  1:98879/1‑150 (MQ=255)
                                                                                                 aTCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGGCTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAACGTGGGGGGCAGACGGCTATTCATCCAGTCCGATTTTAGACCCAGCTCGTAggt                                >  1:469032/1‑148 (MQ=255)
                                                                                                    cACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGGCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAgttgtt                             >  2:165375/1‑150 (MQ=255)
                                                                                                           ttATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTGTCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGt                       <  2:232679/149‑1 (MQ=255)
                                                                                                             ataGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGAtt                   <  2:218653/151‑1 (MQ=255)
                                                                                                                ggCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGa                <  2:469036/151‑1 (MQ=255)
                                                                                                                    gTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCa            <  2:437423/151‑1 (MQ=255)
                                                                                                                       tCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGAAGTTGTTACCGGTGATTGGAg               >  1:494021/1‑145 (MQ=255)
                                                                                                                         gCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGatat       >  1:272343/1‑151 (MQ=255)
                                                                                                                              tgctgcCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGcc  >  2:8482/1‑151 (MQ=255)
                                                                                                                              tgctgcCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGc   <  2:578565/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
TATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCC  >  NC_000913/1159700‑1159976

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