Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,067,7450AG50.0% 35.2 / 20.7 48L511P (CTG→CCG) yihQputative sulpholipid alpha‑glucosidase; alpha‑glucosyl fluoride glucosidase
Reads supporting (aligned to +/- strand):  ref base A (18/5);  new base G (0/23);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 2.39e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTC  >  NC_000913/4067597‑4067891
                                                                                                                                                     |                                                                                                                                                    
aCGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAgg                                                                                                                                                    >  1:495351/1‑151 (MQ=255)
         aTGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTaa                                                                                                                                         <  2:76714/151‑1 (MQ=255)
                aTGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCa                                                                                                                                   >  1:589694/1‑150 (MQ=255)
                             cGCCGTAAAACTGCCTGTTGTCGCCAGGAGGGTTCCCTTGGGGGGTGCGCATCATCGGCGTGAGGGCGCTGAAACCGCACCAGCGCAGCAGCAGCTTTTTGCTGCGCTTCATCCCAAACGGGGTGGTGAAACCGCCAATGTCGCTgtgg                                                                                                                       <  2:382584/149‑1 (MQ=255)
                              gCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGc                                                                                                                    >  1:574237/1‑151 (MQ=255)
                                agTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGc                                                                                                                    >  1:579789/2‑149 (MQ=255)
                                 gTCAAACTGCCGGTTGCCGCAAGGACGTTTACCTTCGTGGGGGCCCTTCATCGGGGTGAGGGCGCTAAAACCGCCCCGGCGCAGCACCGGCTTTTTGCTGCGTTTCATCCAAAACGGGGTGTTGTAACCGCCAATGTCGCTGTGGTGCAg                                                                                                                  <  2:217978/150‑1 (MQ=255)
                                     aaCTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGc                                                                                                                >  1:482244/1‑148 (MQ=255)
                                        tGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTc                                                                                                          >  1:410516/1‑151 (MQ=255)
                                              ttGCCGCCGGGACGGTTCCCTTGGGGGGTGCGCTTCTTCGGGGGGAAGGCGCTGAACTCGCACCGGCGCGGAAGCATTCTTTTGGTGCGTTTCATCCAAAACGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTc                                                                                                     <  1:737688/150‑1 (MQ=255)
                                                       ggcggtttACCTCCGTGGGTGCTCTTATTCGGCGTGAAGGCGCGAAACTCGCACCACCGCAGCAGCAGCTTTTTCCTGCGTTTCTCCTCAAAGAGGGTGGTTTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCg                                                                                           <  2:592617/145‑1 (MQ=255)
                                                       ggACGGTTTCCTTGGTGGGTGCGCATCAGCGGCGTGAAGGCGCTAAAACCGCACCAGCGCAGCAGCAGCTTTTTGCTGCGCTTCATCAAAAAGAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCg                                                                                           <  2:301059/151‑1 (MQ=255)
                                                       ggACGGTTCCCTTCGGGGGGGCGAATCATGGGGGTGAGGGCGGTAAAATCGCCCCGGCGCAGCAGAAGCTTTTTGCTGCTCTTCATCCAAAAGGGGGGGGTTTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGc                                                                                                <  2:12214/146‑1 (MQ=255)
                                                           ggTTCCCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTAAAATCGCCCCAGCGCAGCAGCAGCTCTTTGCTGCGTTTCATCTCAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACa                                                                                        <  2:624112/150‑1 (MQ=255)
                                                                ccTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTAAAATCGCCCCAGCGCAGCAGCACCTCTTTGCTGCGCTTCATCTCAAACGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGcgcc                                                                                   <  2:621913/150‑1 (MQ=255)
                                                                   tCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGccgccg                                                                               >  1:169618/1‑151 (MQ=255)
                                                                       gggTGCGCATCACCGGCGGGAGGGCGCTAAAACCGCCCCAGCGCAGCGGCGGCTTTTTGCTGCGCTCCATCTAAAACGGGGTGGTTAACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCgccgcc                                                                                <  1:727872/146‑1 (MQ=255)
                                                                                tcatcGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCCTTGCCAGCGACAGCGCCGCCGGGACAACCGACGc                                                                  >  1:348914/1‑151 (MQ=255)
                                                                                 catcGGGGTGAAGGCGCTGAAACCGCCCCAGCGAAGCAGCAGCTTTTTGCTGCGCTTCATCTAAAACGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAAc                                                                        <  2:334264/144‑1 (MQ=255)
                                                                                    cGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAgg                                                              >  2:663927/1‑151 (MQ=255)
                                                                                     ggCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGc                                                             >  1:766712/1‑151 (MQ=255)
                                                                                         tGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCAtc                                                         >  2:488930/1‑151 (MQ=255)
                                                                                           agggCGCGAAAATCGCACCGGCGCAGCGGCGCTTTTTGGCGGCGCTTCATCACAAAGGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCAt                                                          <  1:655611/146‑1 (MQ=255)
                                                                                            aggccatGAAATAGCACCAGCGCAGTAACAGCCCTTTGCTGCGCTTCATCTAAAACAGAGTAGTCAAACCCCGAATGACGCTGCGGTGCAGGCCATGGCCGGTCATTGCCAGCGACAGCGCCGCCGCGACAACCGACCCCAGGCCAtcgt                                                       >  1:677237/7‑150 (MQ=255)
                                                                                             ggCGCTAAAACCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACGGGGGGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCAtcgtc                                                      <  2:387238/150‑1 (MQ=255)
                                                                                             gccgcTAAAACCGCACCGGGGCGGCAGCTGCTTTTTGCGGGGTTTACTCACAAACGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCAtcgtc                                                      <  1:513676/148‑1 (MQ=255)
                                                                                                    aaaCCGCACCGGCGCAGCGGAATCTTTTGGCTGCGATTCATCACAAACGGGGTTGTGAACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTca                                              <  1:542097/151‑2 (MQ=255)
                                                                                                     aaTCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTcc                                              >  1:297435/1‑150 (MQ=255)
                                                                                                          cACCAGCGCAGCACCAGCTTTTTGCTGCGCTTCATCTCAAACAGGGTGTTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCg                                         <  2:227844/150‑1 (MQ=255)
                                                                                                            ccAGCGCAGCAGCAGCTTTTGCCTGCGCTTCATCTAAAAGGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCg                                         <  2:128664/148‑1 (MQ=255)
                                                                                                               gcgcAGCGCCAGCTCTTTGCTGCTCTTCATCACAAACGGGGGGGTTTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACg                                      <  1:390677/148‑1 (MQ=255)
                                                                                                                     cagcagcTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcg                             >  1:5319/1‑151 (MQ=255)
                                                                                                                        catctctTTGCTGCGTTTCATCAAAAAGGGGGTGGTGTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTg                                 <  1:732404/141‑1 (MQ=255)
                                                                                                                        cagcTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgccc                          >  2:581173/1‑151 (MQ=255)
                                                                                                                        cagcTCTTGGCTGCTCTGAATCAAAACCGGGTGGGTGAAGCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTgg                                <  1:391530/145‑1 (MQ=255)
                                                                                                                         agcTCTTTGCGGCGTTTCTTCTAAAAAGGGGTGTTTTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGt                               <  1:244254/145‑1 (MQ=255)
                                                                                                                         acctttttGCGGCGTTTCATCAAAAAAGGGGTGGTGAAGCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTgg                                <  2:754833/139‑1 (MQ=255)
                                                                                                                            tttttGCTGCGTTTCATCTCAAACGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTc                              <  2:64749/141‑1 (MQ=255)
                                                                                                                            tctTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCgcccgccc                      >  1:187701/1‑151 (MQ=255)
                                                                                                                                  gctgcGCTTCTTCTAAAACGGGGTGGTATAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCAcat                  <  2:595257/150‑1 (MQ=255)
                                                                                                                                     gtgtttCATCTCAATCGGGTGGTTGTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCAcatcat               <  1:764589/146‑1 (MQ=255)
                                                                                                                                         ttCACCTAAAACGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCAcatcat               <  2:58249/146‑1 (MQ=255)
                                                                                                                                            atccaaaaCGGGGTGTTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGa         <  1:453899/144‑1 (MQ=255)
                                                                                                                                                cAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTAttt   >  2:439950/1‑151 (MQ=255)
                                                                                                                                                cAAACAGGGAGCGGCAACCGCCAATGGCGCTGTGGGGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACCTCATGGTGGAGTAttt   >  1:526473/1‑151 (MQ=255)
                                                                                                                                                 aaaCGGGGTGGTTTACCCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAg        <  1:291143/145‑1 (MQ=255)
                                                                                                                                                 aaaCGGGGTGGTGAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTc  <  2:692661/151‑1 (MQ=255)
                                                                                                                                                 aaaCAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGt       <  2:287900/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTC  >  NC_000913/4067597‑4067891

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