Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,079,4600TG43.6% 89.1 / 14.1 55G52G (GGT→GGGputPproline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base T (4/27);  new base G (24/0);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 9.40e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACC  >  NC_000913/1079314‑1079607
                                                                                                                                                  |                                                                                                                                                           
aGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCg                                                                                                                                                          <  2:987915/150‑1 (MQ=255)
  cacaCCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCg                                                                                                                                                       <  2:289131/151‑1 (MQ=255)
  cacaCCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCg                                                                                                                                                       <  1:95774/151‑1 (MQ=255)
    cacCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTggggtggg                                                                                                                                                     >  1:827233/1‑144 (MQ=255)
        gATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGgcggcg                                                                                                                                                       >  2:1121623/1‑145 (MQ=255)
           gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGc                                                                                                                                              <  1:137510/151‑1 (MQ=255)
           gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCgggggggat                                                                                                                                                  >  2:633482/1‑139 (MQ=255)
              ggTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCgg                                                                                                                                           <  2:370086/150‑1 (MQ=255)
               gTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAgggg                                                                                                                                           >  1:23099/1‑149 (MQ=255)
                    aTTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGGGGGGATGAGCgg                                                                                                                                           >  1:960742/1‑144 (MQ=255)
                       ttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGGATGAGAggccggggg                                                                                                                                    >  2:813406/1‑145 (MQ=255)
                           gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATgg                                                                                                                            <  2:956832/151‑1 (MQ=255)
                              tataTCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGGGGGTGGGTGTTTATtggg                                                                                                                          >  1:1063342/1‑150 (MQ=255)
                               atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATggg                                                                                                                           <  1:263228/148‑1 (MQ=255)
                                          aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTg                                                                                                          <  2:1018873/151‑1 (MQ=255)
                                          aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATGAGCGGGGGGGTGTTTATTGGGGTGTCGGgcg                                                                                                               >  2:876198/1‑148 (MQ=255)
                                             atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGttt                                                                                                     <  1:631871/151‑1 (MQ=255)
                                              tatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGAGGGGGTGGGGGTTTATTGGGGTGGCGGGGGGTGgttt                                                                                                    >  1:156358/1‑151 (MQ=255)
                                                   attgGGTTTATCGCCTGGCGATCAATGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGTGGCGTCGGGATTGAACGGGTGGCTTTTGATATGGTCTCCGCGCGGTGTgtttgttt                                                                                               >  2:85879/1‑151 (MQ=255)
                                                   attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTtttcttt                                                                                               <  1:169543/151‑1 (MQ=255)
                                                   attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTAAGGGGGGCGGGGGGGTGATTGGGGTGGCGGGGGGTGGtttttttt                                                                                               >  2:1046436/1‑151 (MQ=255)
                                                       ggTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGGGGCGGGGTGGTGATTGGGGTGGCGGGGGGTGGtttttttt                                                                                               >  2:201107/1‑147 (MQ=255)
                                                          ttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGCGGGTGGGGGTTTATGGGGGTGGCGGGGGCTGTttttttttttcgg                                                                                          >  2:1110949/1‑144 (MQ=255)
                                                            aTCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGAtt                                                                                      <  2:232532/150‑1 (MQ=255)
                                                             tCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGTGGGGATGGGGGGGCGGCCGGTTATTGGGGTGGCGGGGGGTGGTTTTTTTTCCGGGGttt                                                                                     >  1:423946/1‑150 (MQ=255)
                                                               gCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCg                                                                                  <  1:348787/151‑1 (MQ=255)
                                                                 cTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGGGGGTATGAGCGGGTGGCTGTTTATGGGGGTGCCGGGGGCGGGTTTTTTTTTCGGGGGTTTc                                                                                    >  1:166731/1‑148 (MQ=255)
                                                                    gCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGAATGTGAGGGAGGGTGTTTATGGGGGTGGCGGGGGGCGGTTTTTTTTTCGGGGGTTTcttaa                                                                                >  1:308142/1‑145 (MQ=255)
                                                                     cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCt                                                                            <  1:698056/151‑1 (MQ=255)
                                                                     cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTGTGTACGGCCGGCGGTTCATTGGGGTGCCCGGCGCTGTGTTTTTTTTTCGGGGttttgtaaca                                                                              >  2:563224/1‑143 (MQ=255)
                                                                      gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTg                                                                           <  2:1075703/151‑1 (MQ=255)
                                                                      gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGGGGGGGGGTGTTTATTGGGGTGGCGGGGGCGGTTTTTTTTTTCCGGGTTTTCCaaagatg                                                                            >  1:317656/1‑147 (MQ=255)
                                                                               aaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTGGGGGATGAGGGGGTCGGTGGTAATAGGGGTGGCGGGGGGTGGTTTTTTTTTCCGGGTTTTCGAAAACTGGGTCGcc                                                                    >  1:1048807/1‑149 (MQ=255)
                                                                                 aaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTg                                                                <  2:209912/151‑1 (MQ=255)
                                                                                      tGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCt                                                            <  1:1154209/150‑1 (MQ=255)
                                                                                       gacgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTg                                                           <  2:211105/150‑1 (MQ=255)
                                                                                       gacgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGGGGGGGGGTTATGGGGGTGCGGGGGGGGGGTTTTTTTTTCGGGGGTTTTCAAAGGTGGATCGCCCtttg                                                               >  2:73045/1‑144 (MQ=255)
                                                                                         cgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGGGGGTATGTGAGGGGGGCTGTTTTTGGGGGTGCCGGGGGCGGGTTTTTTTTTCGGGGGTTTCGAAAGGCGGGTAGCCAtttggcatac                                                         >  1:786387/1‑142 (MQ=255)
                                                                                              atatTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTa                                                     <  1:868095/149‑1 (MQ=255)
                                                                                              atatTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAg                                                    <  2:584969/150‑1 (MQ=255)
                                                                                              atatTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTTTGAGGGGCCGGGCGGTTATTGGGTTGGCGGGGGGGGTTTTTTTTTTCCGGGTTTTCGAAAGCGGGATCGCCCTTGGGCTTACCttttg                                                   >  2:210631/1‑148 (MQ=255)
                                                                                                   cTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGAATGTGCGGCCGGGTGGTTATGGGGGTGGCGGGGGGGGTTTTTTTTTCTGCGGTTTTCCGCAGCTGGGCCGCCCTTGGcccgccctt                                                      >  2:403134/1‑136 (MQ=255)
                                                                                                       gCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGAt                                          <  2:436741/151‑1 (MQ=255)
                                                                                                               aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGa                                  <  1:26135/151‑1 (MQ=255)
                                                                                                               aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGCGGGGTCCGATATTAGCGGGCGGGCGTTTATTGGGGTGGCGGGGGCGGGGTTTTTTTTCGGGGTTTTCCGAAGGTGGGGCGCCCTTTGCCTGTCATTTTGGGCGTgg                                            >  2:579798/1‑141 (MQ=255)
                                                                                                                gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGCGGCTGGGCGTTAATGGGGGTGGCGGGGGGTGGTTTTTTTTTCGGGGGTTTCCAAAGACGGGTCGGCATTTGGCTGTCATTTGGGGCCGGGGTTTACTgg                                   >  2:1170883/1‑149 (MQ=255)
                                                                                                                 tCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGGGGGTGGGTGTTTATTGGGGTGCCCGGGGCTGTTTTTTTTTTCGGGGGTTTCGGAAGCTGGGTCGCCCTTTGGCTGTCATTATGGGCGTGGGGTAAATgg                                   >  1:1150098/1‑148 (MQ=255)
                                                                                                                      ggCCATTCGTGGCGGCAGTATCGGCGGGTGCGTCCGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtg                             <  2:559931/149‑1 (MQ=255)
                                                                                                                             cGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCggg                    <  1:171094/151‑1 (MQ=255)
                                                                                                                              gTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGc                   <  1:316673/151‑1 (MQ=255)
                                                                                                                                 aCGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGt                <  2:284452/151‑1 (MQ=255)
                                                                                                                                      aTTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgt           <  1:96680/151‑1 (MQ=255)
                                                                                                                                        tATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtg          <  1:654460/150‑1 (MQ=255)
                                                                                                                                             gcggGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCAt     >  2:524548/1‑150 (MQ=255)
                                                                                                                                               gggTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAcc  >  2:751292/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                           
AGCACACCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACC  >  NC_000913/1079314‑1079607

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