Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG44.7% 71.2 / 13.2 48N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (4/22);  new base G (21/0);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.47e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTG  >  NC_000913/2677893‑2678188
                                                                                                                                                   |                                                                                                                                                      
gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGTGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTc                                                                                                                                                    <  1:200845/151‑1 (MQ=255)
     ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGCGGCCGAGGgttttctcc                                                                                                                                                >  1:416251/1‑145 (MQ=255)
      tttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCGAGGGGGAGGgttttctcccc                                                                                                                                              >  2:94522/1‑144 (MQ=255)
       ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCccccc                                                                                                                                               >  1:779492/1‑149 (MQ=255)
         gcgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTg                                                                                                                                            <  2:618014/150‑1 (MQ=255)
          cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTc                                                                                                                                          <  1:182735/151‑1 (MQ=255)
                gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTTACCGAGGCGAAGGTTTTCCGCCAcctctagc                                                                                                                                       >  1:534006/1‑141 (MQ=255)
                 cTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGGGGCAGAGGGTTTTCCCCCCTGTCAGcccc                                                                                                                                    >  1:926556/1‑149 (MQ=255)
                 cTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGcc                                                                                                                                      >  1:201243/1‑148 (MQ=255)
                   tAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTc                                                                                                                                  <  1:235005/150‑1 (MQ=255)
                    aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTcc                                                                                                                                 >  1:954472/1‑150 (MQ=255)
                          gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAc                                                                                                                         <  2:754107/151‑1 (MQ=255)
                           aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACt                                                                                                                        <  1:94949/151‑1 (MQ=255)
                                aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCCAGGCAGAGGGTTTTCCCCCCTGTCAGCCCCTCCCCCAGCACTTa                                                                                                                      >  2:729955/1‑149 (MQ=255)
                                       cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCGCTGTCAGCCCTTCCACCAGACTTAACTCCTGa                                                                                                              <  2:735272/149‑1 (MQ=255)
                                        aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCAGGCAGAGGGTTTTCGCCACTCTCCGCCCCTCCACCCGACTTAACACCTg                                                                                                               >  2:527422/1‑147 (MQ=255)
                                         ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATaaa                                                                                                          <  2:78465/151‑1 (MQ=255)
                                         ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATaaa                                                                                                          <  1:1016984/151‑1 (MQ=255)
                                          gACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGAACGAGGCAGAGGGTTTTCGCCACTTTCAGCCCCTCCCCCCGCCCTAAATccctgtaaa                                                                                                          >  2:710599/1‑143 (MQ=255)
                                           aCTCACGTCAACGCCCAGCGACTGTAAGCACCGTTGGGTATCCTCCGCCTTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCGCTGACCGAGGCAGAGGTTTACCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATaaa                                                                                                          <  2:419635/149‑1 (MQ=255)
                                             tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGGAGAGGGTTTTCGCCACTGTCAGCCCCTCCCCCAGACTTAACTCCTGATAAAcc                                                                                                        >  1:100263/1‑149 (MQ=255)
                                              aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCGGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGACTTAAATCCTGATAAAcccc                                                                                                      >  2:27134/1‑148 (MQ=255)
                                                 gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCGGGCAGAGGGTTTTTGCCCCTGTCTGCCCTTTCCCCCGACCTAACTCCCGATAAACCgcgcgc                                                                                                  >  2:617866/1‑151 (MQ=255)
                                                 gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCgcgcgc                                                                                                  <  1:33637/151‑1 (MQ=255)
                                                 gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTCCGCCCCTGTCAGCCCTTCCACCCGACTTAACTCCTGATAAAAcccgc                                                                                                    >  1:519114/1‑149 (MQ=255)
                                                    aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGACTTAACTCCTGATAAACCGCGCGcccc                                                                                               >  2:954936/1‑151 (MQ=255)
                                                       gCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGAGGGTTTTCCCCACTGTCCGCCCCTCCACCAGACTTTAATCCTGATAAACCCCGCGCcccccc                                                                                            >  2:1018840/1‑151 (MQ=255)
                                                            gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCCGGGAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGACTTAACTCCTGATAAAACCCGCGCCCCCCCCGtt                                                                                        >  2:234477/1‑150 (MQ=255)
                                                              gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGAGGGTTTTTGCCACTGTCAGCCCCTCCCCCAGAATTAAATCCTGCTAAAACACGCCCAccccc                                                                                            >  2:695140/1‑144 (MQ=255)
                                                                cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTCCGCCACTGTCAGCCCCTCCACCCGACTTAACTCCTTATAAACCGCGCGCACCCCCCCTTCAGCg                                                                                   >  1:335571/1‑151 (MQ=255)
                                                                  gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCg                                                                                   <  2:87708/149‑1 (MQ=255)
                                                                                 cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACCGCGAGTc                                                                       >  1:1012882/1‑146 (MQ=255)
                                                                                  aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTa                                                                 <  2:646547/151‑1 (MQ=255)
                                                                                    ccTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGAGGGTTTTCGCCCCTGTCCGCCCTTCCACCCGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCCGCGGCACGGCGAGTCCGc                                                                    >  1:988405/1‑146 (MQ=255)
                                                                                      tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCGCCCCCTTCAGCCCCTCCACCAGACCTAACTCCCGATAAAACGCGCGCACCCCCCGGTCAGCGGCACGGCGAGGCAGCGTAg                                                                >  2:448785/1‑148 (MQ=255)
                                                                                       gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTc                                                            <  1:1011891/151‑1 (MQ=255)
                                                                                        cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCg                                                           <  2:46209/151‑1 (MQ=255)
                                                                                          ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGGGGCAGAGGGTTTTCCCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCAGCCGAGCTTCGt                                                          >  1:314100/1‑150 (MQ=255)
                                                                                          ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                         <  2:1039180/151‑1 (MQ=255)
                                                                                            aTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCa                                                        <  1:625403/150‑1 (MQ=255)
                                                                                               tGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGGCGAGGATTTTCCCCACAGGCAGACCTTCCACCCCACTTTACTCCTTGTAAACCCCGCGCACCCCCACGTTAACCGCACGGCGGGGCAGCGTAGCTTCCACACCtt                                                    >  1:881194/1‑151 (MQ=255)
                                                                                                            tcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCCCGCCCACCCCCCGTTCAGCGGCACCGCGAGGCAGCGGAGCTTCGTCACCCTCCAGCCGCGtc                                         >  1:690114/1‑149 (MQ=255)
                                                                                                              aTGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCGCCCCTGTCAGCCCTTTCACCAGACTTTAATCCTGATAAAACGCGCGCACCCCCCGTTTAGCGGCACGGCGAGTCAGCGGAGCTTCGTCACCCTCCCGTCGCGCCTCGc                                     >  1:965092/1‑151 (MQ=255)
                                                                                                                      gcggcgGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAg                             <  2:736126/151‑1 (MQ=255)
                                                                                                                           ggggCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                         <  2:954295/150‑1 (MQ=255)
                                                                                                                                  cTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                  <  2:686556/150‑1 (MQ=255)
                                                                                                                                    gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc               >  1:589851/1‑151 (MQ=255)
                                                                                                                                         aGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc          >  1:756177/1‑151 (MQ=255)
                                                                                                                                               ggggTTTTCCGCCACTGTCAACCCTTCCGCCAGACTTAACTCCTGAGAAACCGCGCGCGCCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCTCCGCTATCCGGGCGTTc      <  2:719311/148‑1 (MQ=255)
                                                                                                                                                 ggTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt   <  1:262467/151‑1 (MQ=255)
                                                                                                                                                   ttttCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTg  >  2:189305/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTG  >  NC_000913/2677893‑2678188

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: