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

AGAACAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGA  >  NC_000913/1505872‑1506163
                                                                                                                                                     |                                                                                                                                                
agaACAGGTACGCGATCCCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCa                                                                                                                                                 <  2:471139/150‑1 (MQ=255)
  aaCAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAAtt                                                                                                                                              >  2:947023/1‑151 (MQ=255)
  aaCAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAAtt                                                                                                                                              <  2:39059/151‑1 (MQ=255)
   aCAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTa                                                                                                                                            <  2:1044006/151‑1 (MQ=255)
          cgcgATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGc                                                                                                                                     >  2:92052/1‑151 (MQ=255)
                 ccATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAAc                                                                                                                              <  2:733776/151‑1 (MQ=255)
                    tAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAAcccc                                                                                                                           >  2:631292/1‑151 (MQ=255)
                            gagTGAGGGTACTGGCCCCCATGGCCAGCCCCAGCGCCGTCATCCAGCCAGAGTTTTGTGCGGGGGTGGCCCGGGCGACAATCGCTTTTTGCCAAATTTTTGCCCCTAAAGGGGGGTAGCAAATAAATACCGCGAAAACCCCCGCCAGTa                                                                                                                    <  1:613435/150‑1 (MQ=255)
                                cggCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCACCCAGCCAGAGTTTTGGGCGTGGGGGGCTCCGGAGAAAACCGCGGTTTGCCAGTTTATTCCCGCTGAACTGGCGTAACAAATAAATCCCGCCAGAAACCCCGCCAGTAGCgtg                                                                                                               <  1:641350/149‑1 (MQ=255)
                                    gTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCgtgggtg                                                                                                           >  2:1073676/1‑151 (MQ=255)
                                     tACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCgtgggtgg                                                                                                          >  2:787595/1‑151 (MQ=255)
                                      aCTGCCCCCTATTGCCCCCCCCGGCCCCGCCACCAACCAAGAGTTTGGGGGGGTGGTGGCCCGGCCAAAAACCGCTCCTGCCAAGATTATTGCCGCGAAACGGGGGTACCAATTAAATCCCGCAAGAACCCCCGCCAGTAGCgtgggtg                                                                                                           <  1:761439/149‑1 (MQ=255)
                                        tGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTggagg                                                                                                       >  1:339541/1‑151 (MQ=255)
                                         gACCCCCATGGCCCCCCCCGGCCCCTTCACCCAGCCAGAGATTTGTGCAGGGGTGGCTCCGGCGACACCCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTggagg                                                                                                       <  1:461723/150‑1 (MQ=255)
                                           cGCCCATGGCCAGCCCCGGCGCCGCCCCCCAGCGAGAAATTTGTGCAGGGGTGGCTCCGGCAACAATCGCTTCTTGCCAGATTATTCCCGCTGATCGGGCGTAGCCATTTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTggagg                                                                                                       <  2:164003/148‑1 (MQ=255)
                                                        ccccAGCGCCTTCACCCAGCCAGAGTTTTGTGCGGTGGGGGCCCGGGCAACAACCGCTGCTGCCCAGATTATTCCCGCTGAACGGGGGTAGCCATTTAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGa                                                                                          <  2:975086/148‑1 (MQ=255)
                                                              cgcCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTc                                                                                 >  2:799404/1‑151 (MQ=255)
                                                              ccccGCCACCCGGCCAGAGTTTTGTGCGGGGGGGGCCCCGCCAACATTCGCTGTTTCCAAGATTTTGCCCGCTAACGTGGGGTAGCAAATTAAACCCCCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGt                                                                                  <  2:816961/148‑1 (MQ=255)
                                                                  gTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCt                                                                             >  1:887339/1‑151 (MQ=255)
                                                                       ccagcaagAGATTTGTGCGGTGGGGGCTCCGGCGAAAATCGCTGCTTGCAAGATTATTGCCGCTGAACTGGCGTACCAAATAAATCCCCCCAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCg                                                                           <  2:366732/148‑1 (MQ=255)
                                                                          gccagAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt                                                                     >  2:1039516/1‑151 (MQ=255)
                                                                           ccagAGTTTTGTGCAGGGGTGGCCCCGGCGACAATCCTTGCTTCCCATATTATTGCCGCTGAACTGGCGTAGCAAATTAACCCCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTta                                                                    <  2:69684/151‑1 (MQ=255)
                                                                           caggagtTTTGTGCGGTGGTGGCCCGGGAGAAATCCGCTGTTTGCCAGATTATTGCCGCTGACCTGGCGTACCAAATAAATCCCCCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGt                                                                      <  2:409491/146‑1 (MQ=255)
                                                                              gagtTTTGGGCGGTGGTGCCCCCGGCGACACTCGCTGCTTCCAAGTTTATTGCCGTGGACGGGGGTTACCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt                                                                     <  1:841901/147‑1 (MQ=255)
                                                                              gagtTTTGGGCAGTGGTGGCTCCGGCGACAATCGCTGTTTGCCAGATTATTGCCGCTGAACTGGCGTAGCAAATTAATCCCGCAAAAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                   <  2:475361/149‑1 (MQ=255)
                                                                              gagaTTTGTGCGGGGGTGGCTCGGGCAACAACCGCTGCTGGCCAGATTATTGCCGCTGAACTGGCGTAGCAAATAAACCCCGCAAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTtat                                                                   <  1:1019169/149‑1 (MQ=255)
                                                                                gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                               >  1:286070/1‑151 (MQ=255)
                                                                                gaTTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCg                                                               >  1:837251/1‑151 (MQ=255)
                                                                                  tttGTGCAGTGGTGGCCCCGGCGACAATCGCTGCTTGCCAGATTATTCCCGTTGAACGGGCGTAGCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGtt                                                                     <  2:335558/143‑1 (MQ=255)
                                                                                     gtgCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCAGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTAc                                                            >  2:161613/1‑149 (MQ=255)
                                                                                            ggtggCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGt                                                   >  2:1009542/1‑151 (MQ=255)
                                                                                            ggtggCTCCGGCGAAAATCGCTGTTTGCCAGATTATTGCGGCTGAACGGGCGTAGCAAATAAATCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCAcc                                                     <  1:359496/149‑1 (MQ=255)
                                                                                                  tCCGGCAACAATCGCTGTTTGCCAGTTTTTTGCCGCTGAACGGGTGTCGCCAATTAATCCCGCAAAACACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAg                                                 <  2:194514/147‑1 (MQ=255)
                                                                                                    cGGCGAAAATCGCTGTTTGCTATTTTTTTGCCGCATAACGGGCGTAACAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGc                                                <  1:640615/146‑1 (MQ=255)
                                                                                                     ggCGACATTCGCTGCTTGCCAGTTTATTGCCGCTAACCTGGCGTAGCAATTTAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt                                          <  2:1099319/151‑1 (MQ=255)
                                                                                                     ggCAACAATCGCTGCTTGCCAGATTATTCCCGCGGAACGGGGGTGGCAAATAAACACCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGt                                          <  1:486331/151‑1 (MQ=255)
                                                                                                         aCAACCGCGGCTTGCAAGTTTATTGCGGGTGAGCGGGCGTACCCAATAAATCCCGCCAGAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTcc                                        <  1:813974/149‑1 (MQ=255)
                                                                                                                 tgctTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc                              >  1:149056/1‑151 (MQ=255)
                                                                                                                  gctTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGc                              >  1:658790/1‑150 (MQ=255)
                                                                                                                       ccAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGt                        >  2:651538/1‑151 (MQ=255)
                                                                                                                       caaattttttGCCGTTAAACTGGGTTATCAATTTAAAACCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAg                         <  2:498763/142‑1 (MQ=255)
                                                                                                                       caaaatttttCCCGCTGAGCGGGCGTACCCAATTAACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGcaca                           <  2:770481/144‑1 (MQ=255)
                                                                                                                        cAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTg                       >  1:932917/1‑151 (MQ=255)
                                                                                                                                         aCTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGc      >  2:1011177/1‑151 (MQ=255)
                                                                                                                                           tGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg    >  2:888143/1‑151 (MQ=255)
                                                                                                                                            ggCGTACCAAATAAATCCCGCAAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCg    <  2:780697/150‑1 (MQ=255)
                                                                                                                                             gCGTAGCAAATAAACCCCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCgga  <  2:439759/151‑1 (MQ=255)
                                                                                                                                             gCGGAGCAAATAAACCCCGCCAAAACCCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGtt          <  2:649971/143‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGAACAGGTACGCGATACCATAATGTCAGAGTCAGCGTACTGACGCCCATTGCCAGCCCCAGCGCCGTCATCCAGCCAGAGATTTGTGCAGTGGTGGCTCCGGCGACAATCGCTGCTTGCCAGATTATTGCCGCTGAACTGGCGTAGCCAATTAATACCGCCAGAAACCCCGCCAGTAGCGTGGGTGGAGGAATAGAAAACAGACGCATAGTCACCTCGTGCGTTATAGCGTACAGGCACCGTAGCATTTGTCCGTTATAACGCACAAGTGATAAACTTCCGTTTTGCCGGA  >  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: