Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,506,0200CA39.6% 72.8 / 21.6 48G21C (GGC→TGC) ydcOBenE family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base C (19/10);  new base A (0/19);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 2.05e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAACAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGA  >  NC_000913/1505872‑1506163
                                                                                                                                                     |                                                                                                                                                
agaACAGGTACGCGATCCCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCa                                                                                                                                                 <  2:471139/150‑1 (MQ=255)
  aaCAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAAtt                                                                                                                                              >  2:947023/1‑151 (MQ=255)
  aaCAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAAtt                                                                                                                                              <  2:39059/151‑1 (MQ=255)
   aCAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑a                                                                                                                                            <  2:1044006/151‑1 (MQ=255)
          cgcgATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGc                                                                                                                                     >  2:92052/1‑151 (MQ=255)
                 ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAAc                                                                                                                              <  2:733776/151‑1 (MQ=255)
                    tAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAAcccc                                                                                                                           >  2:631292/1‑151 (MQ=255)
                            gagTGAGGGTACTGGCCCCCATGGCCAGCCCCAGCGCCGTCATCCAGCCAGAGTTTTGTGCGGGGGTGGCCCGGGCGACAATCGCTTTTTGCCAAATTTTTGCCCCTAAAGGGGGG‑TAGCAAATA‑AATACCGCGAAAACCCCCGCCAGTa                                                                                                                    <  1:613435/150‑1 (MQ=255)
                                cggCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCACCCAGCCAGAGTTTTGGGCGTGGGGGGCTCCGGAGAAAACCGCGGTTTGCCAGTTTATTCCCGCTGAACTGGCG‑TAACAAATA‑AATCCCGCCAGAAACCCCGCCAGTAGCgtg                                                                                                               <  1:641350/149‑1 (MQ=255)
                                    gTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCgtgggtg                                                                                                           >  2:1073676/1‑151 (MQ=255)
                                     tACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCgtgggtgg                                                                                                          >  2:787595/1‑151 (MQ=255)
                                      aCTGCCCCCTATTGCCCCCCCCGGCCCCGCCACCAACCAAGAGTTTGGGGGGGTGGTGGCCCGGCCAAAAACCGCTCCTGCCAAGATTATTGCCGCGAAACGGGGG‑TA‑CCAATTAAATCCCGCAAGAACCCCCGCCAGTAGCgtgggtg                                                                                                           <  1:761439/149‑1 (MQ=255)
                                        tGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggagg                                                                                                       >  1:339541/1‑151 (MQ=255)
                                         gACCCCCATGGCCCCCCCCGGCCCCTTCACCCAGCCAGAGATTTGTGCAGGGGTGGCTCCGGCGACACCCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTggagg                                                                                                       <  1:461723/150‑1 (MQ=255)
                                           cGCCCATGGCCAGCCCCGGCGCCGCCCCCCAGCGAGAAATTTGTGCAGGGGTGGCTCCGGCAACAATCGCTTCTTGCCAGATTATTCCCGCTGATCGGGCG‑TAGCCATTT‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTggagg                                                                                                       <  2:164003/148‑1 (MQ=255)
                                                        ccccAGCGCCTTCACCCAGCCAGAGTTTTGTGCGGTGGGGGCCCGGGCAACAACCGCTGCTGCCCAGATTATTCCCGCTGAACGGGGG‑TAGCCATTT‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGa                                                                                          <  2:975086/148‑1 (MQ=255)
                                                              cgcCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                                                                                 >  2:799404/1‑151 (MQ=255)
                                                              ccccGCCACCCGGCCAGAGTTTTGTGCGGGGGGGGCCCCGCCAACATTCGCTGTTTCCAAGATTTTGCCCGCTAACGTGGGG‑TAGCAAATT‑AAACCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGt                                                                                  <  2:816961/148‑1 (MQ=255)
                                                                  gTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt                                                                             >  1:887339/1‑151 (MQ=255)
                                                                       ccagcaagAGATTTGTGCGGTGGGGGCTCCGGCGAAAATCGCTGCTTGCAAGATTATTGCCGCTGAACTGGCG‑TACCAAATA‑AATCCCCCCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg                                                                           <  2:366732/148‑1 (MQ=255)
                                                                          gccagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt                                                                     >  2:1039516/1‑151 (MQ=255)
                                                                           ccagAGTTTTGTGCAGGGGTGGCCCCGGCGACAATCCTTGCTTCCCATATTATTGCCGCTGAACTGGCG‑TAGCAAATT‑AACCCCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTta                                                                    <  2:69684/151‑1 (MQ=255)
                                                                           caggagtTTTGTGCGGTGGTGGCCCGGGAGAAATCCGCTGTTTGCCAGATTATTGCCGCTGACCTGGCG‑TACCAAATA‑AATCCCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGt                                                                      <  2:409491/146‑1 (MQ=255)
                                                                              gagtTTTGGGCGGTGGTGCCCCCGGCGACACTCGCTGCTTCCAAGTTTATTGCCG‑TGGACGGGGGTTACCAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt                                                                     <  1:841901/147‑1 (MQ=255)
                                                                              gagtTTTGGGCAGTGGTGGCTCCGGCGACAATCGCTGTTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCAAATT‑AATCCCGCAAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                   <  2:475361/149‑1 (MQ=255)
                                                                              gagaTTTGTGCGGGGGTGGCTCGGGCAACAACCGCTGCTGGCCAGATTATTGCCGCTGAACTGGCG‑TAGCAAATA‑AACCCCGCAAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                   <  1:1019169/149‑1 (MQ=255)
                                                                                gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                               >  1:286070/1‑151 (MQ=255)
                                                                                gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                               >  1:837251/1‑151 (MQ=255)
                                                                                  tttGTGCAGTGGTGGCCCCGGCGACAATCGCTGCTTGCCAGATTATTCCCGTTGAACGGGCG‑TAGCAAATA‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt                                                                     <  2:335558/143‑1 (MQ=255)
                                                                                     gtgCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCAGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                            >  2:161613/1‑149 (MQ=255)
                                                                                            ggtggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGt                                                   >  2:1009542/1‑151 (MQ=255)
                                                                                            ggtggCTCCGGCGAAAATCGCTGTTTGCCAGATTATTGCGGCTGAACGGGCG‑TAGCAAATA‑AATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAcc                                                     <  1:359496/149‑1 (MQ=255)
                                                                                                  tCCGGCAACAATCGCTGTTTGCCAGTTTTTTGCCGCTGAACGGGTG‑TCGCCAATT‑AATCCCGCAAAACACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAg                                                 <  2:194514/147‑1 (MQ=255)
                                                                                                    cGGCGAAAATCGCTGTTTGCTATTTTTTTGCCGCATAACGGGCG‑TAACAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGc                                                <  1:640615/146‑1 (MQ=255)
                                                                                                     ggCGACATTCGCTGCTTGCCAGTTTATTGCCGCTAACCTGGCG‑TAGCAATTT‑AATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt                                          <  2:1099319/151‑1 (MQ=255)
                                                                                                     ggCAACAATCGCTGCTTGCCAGATTATTCCCGCGGAACGGGGG‑TGGCAAATA‑AACACCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt                                          <  1:486331/151‑1 (MQ=255)
                                                                                                         aCAACCGCGGCTTGCAAGTTTATTGCGGGTGAGCGGGCG‑TACCCAATA‑AATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc                                        <  1:813974/149‑1 (MQ=255)
                                                                                                                 tgctTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc                              >  1:149056/1‑151 (MQ=255)
                                                                                                                  gctTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc                              >  1:658790/1‑150 (MQ=255)
                                                                                                                       ccAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt                        >  2:651538/1‑151 (MQ=255)
                                                                                                                       caaattttttGCCGTTAAACTGGGT‑TATCAATTT‑AAAACCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                         <  2:498763/142‑1 (MQ=255)
                                                                                                                       caaaatttttCCCGCTGAGCGGGCG‑TACCCAATT‑AACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca                           <  2:770481/144‑1 (MQ=255)
                                                                                                                        cAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTg                       >  1:932917/1‑151 (MQ=255)
                                                                                                                                         aCTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc      >  2:1011177/1‑151 (MQ=255)
                                                                                                                                           tGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg    >  2:888143/1‑151 (MQ=255)
                                                                                                                                            ggCG‑TACCAAATA‑AATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg    <  2:780697/150‑1 (MQ=255)
                                                                                                                                             gCG‑TAGCAAATA‑AACCCCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgga  <  2:439759/151‑1 (MQ=255)
                                                                                                                                             gCG‑GAGCAAATA‑AACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt          <  2:649971/143‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGAACAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCG‑TAGCCAATT‑AATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGA  >  NC_000913/1505872‑1506163

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