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

GCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTC  >  NC_000913/1159687‑1159985
                                                                                                                                                      |                                                                                                                                                    
gCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGa                                                                                                                                                      <  2:573088/151‑1 (MQ=255)
   tttCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCTTTGACTTCAAATTCCACCCCTTTACTGACTTTCCCATCCCCCGCTTTATAGGCGGTTTCGCCGTTGCTGCCTGGGATa                                                                                                                                                    <  2:680215/150‑1 (MQ=255)
             tATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTGCCGGGGGTAGGTGTAcccg                                                                                                                                          >  1:400835/1‑148 (MQ=255)
              aTCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGCTTCTGCCGGGGGGAGGGGTTCCgggg                                                                                                                                        >  2:4865/1‑149 (MQ=255)
                    tGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGGAGGGGTACCGgggggatg                                                                                                                                   >  2:599273/1‑146 (MQ=255)
                     gCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACt                                                                                                                                    <  1:414225/148‑1 (MQ=255)
                       aaTATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCg                                                                                                                               <  1:110474/151‑1 (MQ=255)
                        atatAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGa                                                                                                                              <  2:56701/151‑1 (MQ=255)
                                     cgcCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCCGTGGGCTGAGCGACATTTTCCt                                                                                                                     >  2:485495/1‑147 (MQ=255)
                                         aaaTGTCAGCTGCCTGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCGCCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCaaca                                                                                                             <  2:555070/151‑3 (MQ=255)
                                                 gCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGATTTGCCGTTGCTGCCGGGGGTAGGTGTACCGGTGGACTGAGCGAAATTATTCTGGTCAATACGGAAGa                                                                                                      >  1:403642/1‑150 (MQ=255)
                                                    gCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGc                                                                                                  <  1:211330/151‑1 (MQ=255)
                                                      cAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCCATACGGAAGATgg                                                                                                   >  1:665888/1‑148 (MQ=255)
                                                        gTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGCTGGTGCCGGGGGTTAGTGGACAGGTGGACTTCGGGACATTTACCTGGTcgct                                                                                                              >  2:392063/1‑132 (MQ=255)
                                                             cGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCCGTCCTGCCGCGGGAAGGAGTAACGGCGGGCCGAGCGACATTATTCTGATCCACACGGAAggagagtaaaggtg                                                                                          >  2:139292/1‑137 (MQ=255)
                                                              ggTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACgtggtg                                                                                        <  2:439039/151‑1 (MQ=255)
                                                                              aGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCCTTTCTGCCGGGGGTAGGGGTACCCGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGGCAGACGGCTAtt                                                                          >  2:50775/1‑149 (MQ=255)
                                                                                ttCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCCCCCTTTCTGCCCGGGGTAGGGGTACCCGTGGGCCGGGCGGCATTTTCCTTCTCAATACGGGAGATGGGTTACGGGGTGGGCAGACgggcatt                                                                          >  2:657085/1‑142 (MQ=255)
                                                                                  cAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCCGTGCTGCCGGGGGTAGGGGGACCGGTGGACTGAGCGACATTTTCCTGCTCAATACGGGAGATGGCTAACGTGGTGGTCAGACGGCTATtcctcc                                                                     >  1:111205/1‑150 (MQ=255)
                                                                                       tCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTcc                                                                <  1:232608/150‑1 (MQ=255)
                                                                                         cACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCCGTGGACTGGGCGACATTATTCCGCTCAATACGGGAGATGGCTGACGGGGGGGGCAGACGGCCATTCATcccg                                                                   >  1:291738/1‑143 (MQ=255)
                                                                                            cccTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGGGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGGGGGCAGACGGCTATTCATCCAGTCCGAttt                                                           >  1:189191/1‑150 (MQ=255)
                                                                                             ccTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTATCCTGCTCCATACGGGAGATGGCTAACGTGGGGGTCAGACGGCTATTCATCCAGTCCGAttt                                                           >  1:488676/1‑149 (MQ=255)
                                                                                                ttACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGGGGGCAGACGGCTATTCATCCCGTCCGATTTTAGa                                                       >  2:268447/1‑150 (MQ=255)
                                                                                                   ctgactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGGGGACTGAGCGACATTATCCTGCTCCATACGGGAGATGGCTAAAGTGGGGGTCAGACGGGTATTCATCCAGTCCGATTTCAGAccc                                                    >  2:495414/1‑150 (MQ=255)
                                                                                                     gactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTTCCGGGGATAGGCGTACCGGGGGACTGGGCGACATTATCCTGCTCCATACGGGAGATGGCTAACGAGGTGGTCAGACGGCTATTCATCCAGTCCGCTTTCAGACCCAg                                                  >  1:36866/1‑150 (MQ=255)
                                                                                                        tgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTATCCTGCTCCATACGGGAGATGGGTAACGTGGGGGGCAGACGGCTATTCATCCAGTCCGATTTCAGAccc                                                    >  2:143227/1‑145 (MQ=255)
                                                                                                           cccATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTa                                            <  2:626126/150‑1 (MQ=255)
                                                                                                                    cGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCg                                  <  1:442313/151‑1 (MQ=255)
                                                                                                                                ggTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGcc                      >  2:512021/1‑151 (MQ=255)
                                                                                                                                    tCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAg                    <  1:541017/149‑1 (MQ=255)
                                                                                                                                      gCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGatat                <  1:322619/151‑1 (MQ=255)
                                                                                                                                        cGTTGCTGCCGGGGATAGGTGTGCCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATAttt              <  1:158414/151‑1 (MQ=255)
                                                                                                                                                   gggATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGt   <  2:439766/151‑1 (MQ=255)
                                                                                                                                                    ggATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc  >  1:142211/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTC  >  NC_000913/1159687‑1159985

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