Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,079,4680AG42.6% 78.8 / 14.0 55D55G (GAT→GGT) putPproline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (2/29);  new base G (23/0);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 2.76e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAAT  >  NC_000913/1079323‑1079611
                                                                                                                                                 |                                                                                                                                                       
aTGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGa                                                                                                                                                    <  1:47240/151‑1 (MQ=255)
  gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAg                                                                                                                                                   <  2:125677/150‑1 (MQ=255)
     ggTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGGATGAGCgg                                                                                                                                                >  1:931388/1‑150 (MQ=255)
      gTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCggcc                                                                                                                                              <  1:149466/151‑2 (MQ=255)
          cATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCggctggct                                                                                                                                         <  2:920679/151‑1 (MQ=255)
          cATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCggctggct                                                                                                                                         <  1:836157/151‑1 (MQ=255)
              ttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGTGCgggggggggtt                                                                                                                                       >  2:245358/1‑141 (MQ=255)
               tgtgTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTa                                                                                                                                   <  2:261458/150‑1 (MQ=255)
                      atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGtt                                                                                                                           <  1:297176/151‑1 (MQ=255)
                            ttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCgg                                                                                                                      <  1:549713/150‑1 (MQ=255)
                                cATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCt                                                                                                               <  2:932677/151‑1 (MQ=255)
                                cATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGATGAGGGGGTGGGTGGTTATTGGGGTTGCggggg                                                                                                                    >  2:527535/1‑148 (MQ=255)
                                  tGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGCGGGTGGCTGTTTATGGGGGTGGCggggg                                                                                                                   >  2:721835/1‑145 (MQ=255)
                                     tatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGGCTGAGGGGGGCGGTGTTTATTGGGGTGTCGGgcgggggttt                                                                                                           >  2:723408/1‑144 (MQ=255)
                                     tatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGCGGCGGGCGGTTTATGGGGGTGCCGGGGGCTGtttt                                                                                                         >  2:144497/1‑151 (MQ=255)
                                       ttgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGGTGGGTGTTTATGGGGGTGGCGGGGGCTGgttttt                                                                                                       >  2:29909/1‑150 (MQ=255)
                                         gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTtttct                                                                                                      <  2:148624/150‑1 (MQ=255)
                                          attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTtttctt                                                                                                     <  2:132453/150‑1 (MQ=255)
                                            tgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTtttctttc                                                                                                  <  2:942544/150‑1 (MQ=255)
                                                   aTCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGAGCGGCTGGGTGTTTATGGGGGTGGCGGGGGCTGTttttttttttgggggtt                                                                                           >  2:668867/1‑142 (MQ=255)
                                                       ccTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGAGCGGGGGGCGGTTTATGGGGGTGTCGGGGGCTGTTTTTTTTTTCGGGGttt                                                                                          >  2:681886/1‑147 (MQ=255)
                                                            cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCt                                                                                <  1:303623/151‑1 (MQ=255)
                                                             gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTg                                                                               <  2:183180/151‑1 (MQ=255)
                                                                 aaCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGTCGGGGTTGGGGGGGTGGGTGGTTATTGGGGTTGCGGGGGCGGGTTTTTTTTTCCGGGTTTTCGGaaggttg                                                                               >  2:266925/1‑143 (MQ=255)
                                                                 aaCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATc                                                                           <  1:163110/151‑1 (MQ=255)
                                                                 aaCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGGGGGGGTGGGTGTTTATGGGGGTGGCGGGGGGTGTTTTTTTTTTCGGGGGTTTCGCAAGCTg                                                                               >  2:720497/1‑147 (MQ=255)
                                                                    gAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGTTTTGCGGGGGGCGGTTAATTGGGGGGTCGGGGGCCGGTTTTTTTTTTCGGGTTTTCGGAAAgtggaagg                                                                           >  1:565019/1‑142 (MQ=255)
                                                                       aaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGGGGGTATGAGCGGGCGGCGTTTTATTGGGGTGGCGGGGGGCGGTTTTTTTTTCCGGGGTTTTGGaaacttggtcggcctt                                                                      >  2:686336/1‑138 (MQ=255)
                                                                        aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTg                                                                    <  2:121941/151‑1 (MQ=255)
                                                                        aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTg                                                                    <  1:613282/151‑1 (MQ=255)
                                                                             tGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGGGCGGGTGGGGGTTGATGGGGGGGGCGGGCGGCGGGTTTTTTTTCTGGGTTTTCTGaagatgggttg                                                                           >  1:20759/1‑131 (MQ=255)
                                                                                cgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGCTGGGTGTTTATTGGGGTGTCGGGGGCTGTTTTTTTTTTCGGGGGTTTCGCAAACTGGGGCGCCATTGGCCTGAc                                                             >  2:1045342/1‑150 (MQ=255)
                                                                                cgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGTTGTGCGGCTGGCTGTTTATTGGGGTTGCGGGGGGCGTTTTTTTTTTTTGGGGTTTCCGAAGATGGGTCGccccttggctgtcc                                                            >  1:821728/1‑139 (MQ=255)
                                                                                    tataTTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGTATTAGGGGGCGGGTGGTAATGGGGGGTCCGGGGGGGGTTTTTTTTTTTGGGGTTTTCCCAAGGTGGATCCCCCCTGGCCCGTCAtt                                                          >  2:1022599/1‑149 (MQ=255)
                                                                                       atTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGTGGTCCGGGATTGGCGGGCGGCGGTTTAATGGGGTTGCTGGCGGCGTTTTTTTTTTTTGTGGTTTCCACAGCTGGGTAGCCCTTTGCCTGAccatt                                                         >  1:763338/1‑143 (MQ=255)
                                                                                        ttCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATTGGGGGGTGGGTGGTTATTGGGGTGCCGGGGGGTGGTTTTTTTTTCCGGGTTTTCGGAAGGTGGGTTGGCCTTTGGctggcatttaggg                                                      >  2:693651/1‑138 (MQ=255)
                                                                                         tCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATTAGGGGGTGGGTGTTTATTGGGGTGTCGGGGGCCGTTTTTTTTTTCGGGGGTTTCGGAAACTGGAGCGCCCGTGGGCTGACATTAAGggcg                                                   >  2:605979/1‑151 (MQ=255)
                                                                                             ggCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTg                                                 <  2:16508/149‑1 (MQ=255)
                                                                                              gCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGTATGAGCGGGCGGGTGTTTATGGGGGTGGCGGGGGCGGGTTTTTTTTTCGGGGGTTTCGAAAGCTGGGTCGCCCTTTGCCTGACATTAGGCGCGTggg                                               >  2:211439/1‑149 (MQ=255)
                                                                                              gCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGAt                                              <  2:582022/151‑1 (MQ=255)
                                                                                                ggTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGCGGGGGGCTGGTTATTGGGGGTGCGGGCGCGGTTTTTTTTTCCCGGGGTTTCGGAAACTGGGTGGCCCTTGGCCTGACATTAGGGGCGGGGGtt                                             >  1:263982/1‑150 (MQ=255)
                                                                                                 gTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTaa                                           <  2:244775/151‑1 (MQ=255)
                                                                                                      aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGACCGCTATGAGTGGGGGGGGGATAATGGGGGTCCCCGGCGCTGTTTTTTTTTTTGGGGTTTGCCGAAAGTGGAACGCCCTTTGCCCGCCCTTAAGCGCGGGGGtttaatcg                                       >  2:876772/1‑144 (MQ=255)
                                                                                                             ggCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtggt                               <  1:149725/151‑1 (MQ=255)
                                                                                                              gCCATTCGTGACGGCATTATCGGCGGGTGGGTCCGGTCTGAGAGGGTGGCCGTTCATAGGGGTGCGGGGGGGTGTTTTTTTTTTCGGGGGTTTCGCAACCTGGGGAGCCCTTGGGCTGGAATTAGGCGCGGGGGTTTAATgggaagtgg                                >  2:999271/1‑142 (MQ=255)
                                                                                                                    cGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCgg                         <  1:280807/150‑1 (MQ=255)
                                                                                                                         cGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGtt                   <  2:669725/151‑1 (MQ=255)
                                                                                                                               tATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtg              <  1:969714/150‑1 (MQ=255)
                                                                                                                                aTCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGtgtg            <  2:478076/151‑1 (MQ=255)
                                                                                                                                aTCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGtgtg            <  1:861363/151‑1 (MQ=255)
                                                                                                                                  cggcggGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGc           <  1:797179/150‑1 (MQ=255)
                                                                                                                                    gcggGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCAt         >  1:40002/1‑150 (MQ=255)
                                                                                                                                         tGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGa    >  2:612105/1‑150 (MQ=255)
                                                                                                                                          gCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGaa   <  2:55478/150‑1 (MQ=255)
                                                                                                                                          gCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAAt  <  2:317766/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                       
ATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAAT  >  NC_000913/1079323‑1079611

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