Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,8370AG39.1% 64.6 / 12.0 47I572T (ATC→ACC) fhuEferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (7/21);  new base G (18/0);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 2.24e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCAC  >  NC_000913/1159695‑1159987
                                                                                                                                              |                                                                                                                                                      
ttCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTGCCCGGGGTAggag                                                                                                                                                  >  2:747803/1‑147 (MQ=255)
ttCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGATTCGCCGGTGCTGCCGGGGGTAgggg                                                                                                                                                  >  2:294966/1‑147 (MQ=255)
   gTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTAcc                                                                                                                                              >  1:974674/1‑150 (MQ=255)
           ctGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTGGAGTTCAAATTCCACCCCTTTACGGTGGGTCGCATCCACCGCTTTATAGGCGCTTTCGCCGTTGCTGCCCGGGATAGGTGTACCTGTTGAc                                                                                                                                       <  2:703076/149‑1 (MQ=255)
              cAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGGTACCCGTGGACTGa                                                                                                                                    >  2:157140/1‑149 (MQ=255)
               aaTATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCg                                                                                                                                 <  1:999764/151‑1 (MQ=255)
                atatAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGc                                                                                                                                  <  1:88757/149‑1 (MQ=255)
                   taGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTACCGGGGGACTGAGCGGCa                                                                                                                              >  1:779449/1‑150 (MQ=255)
                            gcgcCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTTCTGCCGGGGGTAGGTGTACCGGGGGGCTGAGCGACATTATcc                                                                                                                        >  1:146925/1‑147 (MQ=255)
                                 aaaTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGGCTGAGCGACATTATCCTGCTcca                                                                                                                 >  1:25267/1‑147 (MQ=255)
                                   aTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCCGTGCTGCCGGGGGTAGGGGTACCCGGGGGCTGAGCGACATTATCCTGCTCAATACg                                                                                                             >  2:1080142/1‑151 (MQ=255)
                                       cAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACgg                                                                                                            >  2:294399/1‑148 (MQ=255)
                                       cAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGGGGACTGAGCGACATTATCCTGCTCAATACgg                                                                                                            >  1:219305/1‑148 (MQ=255)
                                        aGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCCCCCTTTGTGACGGGGTATGGGGGACCGGGGGGCTGAGAGACATTATCCTGCTCCAAACg                                                                                                             >  2:1110565/1‑146 (MQ=255)
                                        aGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACGGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCTGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGa                                                                                                        <  2:698850/151‑1 (MQ=255)
                                              cAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGCTTCGCCGGTGATACCGGGGGAAGGTGTACCGGTGGGCTTGGCGAGACTCTCCCCCTCAAAACCGGAGAt                                                                                                       >  1:207540/1‑146 (MQ=255)
                                               aGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTaa                                                                                                 <  2:278698/151‑1 (MQ=255)
                                                gTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCTCCGTTTCTGCCGGGGGTAGGGGTACCCGTGGACTGGGCGACATTATCCTGCTCCATACGGGAGGTGGCTAAc                                                                                                >  2:358893/1‑151 (MQ=255)
                                                 ttgttgGTAATTGCGACGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGGAGGGGTACCGGTGGGCTGGGCGACATTATACTGCTCAATTCGGGAGAAGGGTAACg                                                                                               >  2:745823/1‑151 (MQ=255)
                                                     cGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACgtgg                                                                                            <  2:207869/150‑1 (MQ=255)
                                                              cgcCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTACCGGGGGACTGAGGGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGGGGGCAg                                                                                     >  1:480250/1‑148 (MQ=255)
                                                               gcCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGAc                                                                                   <  2:636263/149‑1 (MQ=255)
                                                                   ttGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGGTACCGGGGGGCTGGGCGACATTTTCCTTCTCCATACGGGAGATGGGTAAAGGGGTGGGCAGACgg                                                                                 >  1:257221/1‑147 (MQ=255)
                                                                    tGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCCTTGCTGCCGGGGATAGGTGTTCCGGTGGACTGCGGGACATTTTCCTGCTCAATACGGGACATGGCTAACGGGGGGGTTAGACgggcatt                                                                            >  1:292409/1‑146 (MQ=255)
                                                                          cAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCCTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCCATACGGAAGATGGGTAACGTGGTGGTCAGACGGCTATTAATcc                                                                       >  2:824949/1‑150 (MQ=255)
                                                                          cAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATcc                                                                       <  2:791769/150‑1 (MQ=255)
                                                                                ccACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCg                                                                 <  2:343376/150‑1 (MQ=255)
                                                                                     ccTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGTGTACCGGTGGGCTGAGCGACATTTTCCTGCTCAATACGGAAGATGGCTAACGTGGGGGGCAGACGGCTATTCATACAGTCCGAttt                                                             >  1:436557/1‑149 (MQ=255)
                                                                                      cTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGGGGGCAGACGGCTATTCATCCAGTCCGAttt                                                             >  1:675430/1‑148 (MQ=255)
                                                                                            tgactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCCTTGTTGCCCGGGGTAGGTGTACCGGCGGAGTGAGCGGCATTATCCTGCTCCAAAAGGGAGAAGGGCAACGAGGTGGGCAGAAGGCTATTCATCCCGTCCGATTTTAGAccc                                                      >  1:384906/1‑149 (MQ=255)
                                                                                             gactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCa                                                     <  1:72579/149‑1 (MQ=255)
                                                                                                tgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGGGGGCTGAGCGACATTATCCTGCTCAAAACGGAAGATGGCTAACGGGGGGGGCAGACGGCTATTCATCCAGTCCGATTTTAGACCCAGc                                                   >  1:479286/1‑148 (MQ=255)
                                                                                                 gTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTc                                                 <  2:417874/149‑1 (MQ=255)
                                                                                                      aTCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAgttg                                          <  1:813784/151‑1 (MQ=255)
                                                                                                        cctccGCTTTATAGGCGGTTTCGCCGTTGCTGCCCGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTGTTCATCCAGTCCGATCTCAGACCCAGCTCGTAgttg                                          <  2:890037/146‑1 (MQ=255)
                                                                                                         cACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCCGTGGACTGAGCGACATTATCCTGCTCCATACGGAAGATGGCTAAAGGGGGGGGCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCCCGTAgttgtt                                        >  2:17462/1‑150 (MQ=255)
                                                                                                           ccGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTAcc                                     <  1:985977/151‑1 (MQ=255)
                                                                                                                    aGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTg                             <  2:252283/150‑1 (MQ=255)
                                                                                                                      gCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAg                          <  1:246804/151‑1 (MQ=255)
                                                                                                                            tCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGc                         >  2:306771/1‑146 (MQ=255)
                                                                                                                             cGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGat                    <  1:287395/150‑1 (MQ=255)
                                                                                                                             cGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGat                    <  1:544922/150‑1 (MQ=255)
                                                                                                                                  ttGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGc              >  2:370967/1‑151 (MQ=255)
                                                                                                                                        ccGGTGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGaa         <  2:95980/150‑1 (MQ=255)
                                                                                                                                         cGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACt       <  2:6567/151‑1 (MQ=255)
                                                                                                                                           gggATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTg      >  2:370154/1‑150 (MQ=255)
                                                                                                                                             gATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCa   <  1:744546/151‑1 (MQ=255)
                                                                                                                                              aTAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCAc  <  1:370305/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
TTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCAC  >  NC_000913/1159695‑1159987

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