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

GCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTG  >  NC_000913/2677893‑2678188
                                                                                                                                                   |                                                                                                                                                      
gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGTGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTc                                                                                                                                                    <  1:200845/151‑1 (MQ=255)
     ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGCGGCCGA‑GGgttttctcc                                                                                                                                                >  1:416251/1‑145 (MQ=255)
      tttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCGAGGGGGA‑GGgttttctcccc                                                                                                                                              >  2:94522/1‑144 (MQ=255)
       ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCccccc                                                                                                                                               >  1:779492/1‑149 (MQ=255)
         gcgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTg                                                                                                                                            <  2:618014/150‑1 (MQ=255)
          cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTc                                                                                                                                          <  1:182735/151‑1 (MQ=255)
                gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTTACCGAGGCGAA‑GGTTTTCCGCCAcctctagc                                                                                                                                       >  1:534006/1‑141 (MQ=255)
                 cTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGGGGCAGA‑GGGTTTTCCCCCCTGTCAGcccc                                                                                                                                    >  1:926556/1‑149 (MQ=255)
                 cTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGcc                                                                                                                                      >  1:201243/1‑148 (MQ=255)
                   tAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTc                                                                                                                                  <  1:235005/150‑1 (MQ=255)
                    aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTTCGCCCCTGTCAGCCCTTcc                                                                                                                                 >  1:954472/1‑150 (MQ=255)
                          gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑Ac                                                                                                                         <  2:754107/151‑1 (MQ=255)
                           aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACt                                                                                                                        <  1:94949/151‑1 (MQ=255)
                                aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCCAGGCAGA‑GGGTTTTCCCCCCTGTCAGCCCCTCCCCCAGCACTTa                                                                                                                      >  2:729955/1‑149 (MQ=255)
                                       cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCGCTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGa                                                                                                              <  2:735272/149‑1 (MQ=255)
                                        aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCAGGCAGA‑GGGTTTTCGCCACTCTCCGCCCCTCCACCCG‑ACTTAACACCTg                                                                                                               >  2:527422/1‑147 (MQ=255)
                                         ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATaaa                                                                                                          <  2:78465/151‑1 (MQ=255)
                                         ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATaaa                                                                                                          <  1:1016984/151‑1 (MQ=255)
                                          gACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGAACGAGGCAGA‑GGGTTTTCGCCACTTTCAGCCCCTCCCCCCG‑CCCTAAATccctgtaaa                                                                                                          >  2:710599/1‑143 (MQ=255)
                                           aCTCACGTCAACGCCCAGCGACTGTAAGCACCGTTGGGTATCCTCCGCCTTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCGCTGACCGAGGCAGA‑GGTTTACCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATaaa                                                                                                          <  2:419635/149‑1 (MQ=255)
                                             tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGGAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCCCCAG‑ACTTAACTCCTGATAAAcc                                                                                                        >  1:100263/1‑149 (MQ=255)
                                              aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCGGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCCG‑ACTTAAATCCTGATAAAcccc                                                                                                      >  2:27134/1‑148 (MQ=255)
                                                 gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCGGGCAGA‑GGGTTTTTGCCCCTGTCTGCCCTTTCCCCCG‑ACCTAACTCCCGATAAACCgcgcgc                                                                                                  >  2:617866/1‑151 (MQ=255)
                                                 gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCgcgcgc                                                                                                  <  1:33637/151‑1 (MQ=255)
                                                 gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTCCGCCCCTGTCAGCCCTTCCACCCG‑ACTTAACTCCTGATAAAAcccgc                                                                                                    >  1:519114/1‑149 (MQ=255)
                                                    aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCCG‑ACTTAACTCCTGATAAACCGCGCGcccc                                                                                               >  2:954936/1‑151 (MQ=255)
                                                       gCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGA‑GGGTTTTCCCCACTGTCCGCCCCTCCACCAG‑ACTTTAATCCTGATAAACCCCGCGCcccccc                                                                                            >  2:1018840/1‑151 (MQ=255)
                                                            gcgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCCGGGAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCCG‑ACTTAACTCCTGATAAAACCCGCGCCCCCCCCGtt                                                                                        >  2:234477/1‑150 (MQ=255)
                                                              gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGA‑GGGTTTTTGCCACTGTCAGCCCCTCCCCCAG‑AATTAAATCCTGCTAAAACACGCCCAccccc                                                                                            >  2:695140/1‑144 (MQ=255)
                                                                cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTCCGCCACTGTCAGCCCCTCCACCCG‑ACTTAACTCCTTATAAACCGCGCGCACCCCCCCTTCAGCg                                                                                   >  1:335571/1‑151 (MQ=255)
                                                                  gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCg                                                                                   <  2:87708/149‑1 (MQ=255)
                                                                                 cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCCTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACCGCGAGTc                                                                       >  1:1012882/1‑146 (MQ=255)
                                                                                  aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTa                                                                 <  2:646547/151‑1 (MQ=255)
                                                                                    ccTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGCAGA‑GGGTTTTCGCCCCTGTCCGCCCTTCCACCCG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCCGCGGCACGGCGAGTCCGc                                                                    >  1:988405/1‑146 (MQ=255)
                                                                                      tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGA‑GGGTTTTCGCCCCCTTCAGCCCCTCCACCAG‑ACCTAACTCCCGATAAAACGCGCGCACCCCCCGGTCAGCGGCACGGCGAGGCAGCGTAg                                                                >  2:448785/1‑148 (MQ=255)
                                                                                       gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTc                                                            <  1:1011891/151‑1 (MQ=255)
                                                                                        cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCg                                                           <  2:46209/151‑1 (MQ=255)
                                                                                          ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGGGGCAGA‑GGGTTTTCCCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCAGCCGAGCTTCGt                                                          >  1:314100/1‑150 (MQ=255)
                                                                                          ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                         <  2:1039180/151‑1 (MQ=255)
                                                                                            aTTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCa                                                        <  1:625403/150‑1 (MQ=255)
                                                                                               tGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGGCGA‑GGATTTTCCCCACAGGCAGACCTTCCACCCC‑ACTTTACTCCTTGTAAACCCCGCGCACCCCCACGTTAACCGCACGGCGGGGCAGCGTAGCTTCCACACCtt                                                    >  1:881194/1‑151 (MQ=255)
                                                                                                            tcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGGTTTTCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCCCGCCCACCCCCCGTTCAGCGGCACCGCGAGGCAGCGGAGCTTCGTCACCCTCCAGCCGCGtc                                         >  1:690114/1‑149 (MQ=255)
                                                                                                              aTGCCGTTGCGGCGGGGCCACTGACCGAGGGAGA‑GGGTTTTCGCCCCTGTCAGCCCTTTCACCAG‑ACTTTAATCCTGATAAAACGCGCGCACCCCCCGTTTAGCGGCACGGCGAGTCAGCGGAGCTTCGTCACCCTCCCGTCGCGCCTCGc                                     >  1:965092/1‑151 (MQ=255)
                                                                                                                      gcggcgGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAg                             <  2:736126/151‑1 (MQ=255)
                                                                                                                           ggggCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                         <  2:954295/150‑1 (MQ=255)
                                                                                                                                  cTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                  <  2:686556/150‑1 (MQ=255)
                                                                                                                                    gACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc               >  1:589851/1‑151 (MQ=255)
                                                                                                                                         aGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc          >  1:756177/1‑151 (MQ=255)
                                                                                                                                               ggggTTTTCCGCCACTGTCAACCCTTCCGCCAG‑ACTTAACTCCTGAGAAACCGCGCGCGCCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCTCCGCTATCCGGGCGTTc      <  2:719311/148‑1 (MQ=255)
                                                                                                                                                 ggTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt   <  1:262467/151‑1 (MQ=255)
                                                                                                                                                   ttttCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTg  >  2:189305/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGA‑GGTTTTCCGCCACTGTCAGCCCTTCCACCAG‑ACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTG  >  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: ‑