Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,079,4680AG50.9% 58.3 / 81.0 53D55G (GAT→GGT) putPproline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (3/23);  new base G (27/0);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 4.17e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-10

CGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACA  >  NC_000913/1079321‑1079614
                                                                                                                                                   |                                                                                                                                                         
cGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGatat                                                                                                                                                        <  1:632880/151‑1 (MQ=255)
     ttGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGCg                                                                                                                                                  >  1:498253/1‑151 (MQ=255)
       ggTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCggc                                                                                                                                                <  2:211366/151‑1 (MQ=255)
       ggTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCgg                                                                                                                                                 <  2:582081/150‑1 (MQ=255)
        gTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCggct                                                                                                                                               <  1:779913/151‑1 (MQ=255)
            cATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCgcctgg                                                                                                                                            <  1:128118/149‑1 (MQ=255)
             aTTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTGGGGTTTGGGCgggtgg                                                                                                                                            >  2:791799/1‑148 (MQ=255)
               tttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGAGATGAGCGgctggctggt                                                                                                                                       >  1:75138/1‑149 (MQ=255)
                ttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGAGTCGGGTATCAGCggcgggtcgtt                                                                                                                                       >  2:155369/1‑145 (MQ=255)
                   tgtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGGGGCGGGCTGTTTAt                                                                                                                                  >  1:657476/1‑150 (MQ=255)
                   tgtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGGCggcgggcg                                                                                                                                          >  1:194106/1‑144 (MQ=255)
                     tCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGGATGAGCGGCTGGGTGTTAAttg                                                                                                                                >  1:153134/1‑148 (MQ=255)
                       tataTCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGTATGAGCGGGCGGGTGTTTATtggg                                                                                                                              >  2:282567/1‑150 (MQ=255)
                        atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATgggg                                                                                                                              <  2:604590/149‑1 (MQ=255)
                        atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGtt                                                                                                                            <  1:554935/151‑1 (MQ=255)
                              ttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCg                                                                                                                        <  2:351024/149‑1 (MQ=255)
                               tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGGGCGGGTGGCTGTTAATTGGGGTGCCggg                                                                                                                      >  1:632926/1‑150 (MQ=255)
                                   aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGCGGGTGGGTGGTTATTGGGTTGGCggggg                                                                                                                    >  2:796408/1‑146 (MQ=255)
                                     gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGCGGCCGGGTGTTTATTGGGGTGGCGGGGGCTGtt                                                                                                            >  1:172236/1‑151 (MQ=255)
                                        attgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGttttt                                                                                                         <  1:849489/151‑1 (MQ=255)
                                           gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTtttctt                                                                                                      <  2:364497/151‑1 (MQ=255)
                                           gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTtttct                                                                                                       <  2:905385/150‑1 (MQ=255)
                                           gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGCCGGGGGTCGGCGGCCGGCTGTTAATTGGGGTTTCCGGGCCCGtttttttt                                                                                                      >  1:452715/1‑151 (MQ=255)
                                                  tttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGATATGAGAGGGTGGGTGTTTATGGGGGTGGCGGGGGCTGTtttttttttccg                                                                                                 >  2:206096/1‑145 (MQ=255)
                                                   ttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGa                                                                                               <  1:754039/150‑1 (MQ=255)
                                                    tATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGAt                                                                                              <  2:704177/150‑1 (MQ=255)
                                                            ggCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAg                                                                                     <  2:495833/151‑1 (MQ=255)
                                                            ggCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAg                                                                                     <  2:178079/151‑1 (MQ=255)
                                                              cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGAGCGGCCGGCTGTTTATGGGGGTGCCGGGGGGTGTTTTTTTTTTCGGGGTTTctccaagct                                                                                   >  1:851639/1‑143 (MQ=255)
                                                              cGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATTAGGGGCGGGCTGTTTATGGGGGTGCCGGGGGCCGTTTTTTTTTTTTGGGTTTTCGaaa                                                                                      >  1:279557/1‑148 (MQ=255)
                                                               gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTAGGGAATGGGCGGGCGGCTGTTTATTGGGGTTTAGGGGGGCGTTTTTTTTTTTGGGGGTTTCCgga                                                                                       >  1:478598/1‑145 (MQ=255)
                                                                 tCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATTGGCGGGCGGCTGTTTATGGGGGTGTCGGGGGCTGTTTTTTTTTTCGGGGGTTTCGAAAACttg                                                                                 >  1:873701/1‑148 (MQ=255)
                                                                  cAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGa                                                                                <  1:791721/150‑1 (MQ=255)
                                                                   aaCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGAGGGGGCGGCTGGTTATTGGGTGGCCGGGGGGGGTTTTTTTTTCCGGGGTTTTCGGAAGCTgggtg                                                                              >  2:887571/1‑148 (MQ=255)
                                                                    aCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGGGGGTGGGTGTTTATTGGGGTGGCGGGCGCTGTTTTTTTTTTCGGGGTTTTCGAAAGCTggggcg                                                                             >  2:144025/1‑147 (MQ=255)
                                                                        aaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCAt                                                                         <  2:326233/151‑1 (MQ=255)
                                                                        aaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCCGGTTTGAGCGGGTGGGCGGTTATTGGGGTGGCGGGGGCTGGTTTTTTTTTCCGGGATTTCGGaaggtgggacgcc                                                                            >  2:544999/1‑138 (MQ=255)
                                                                               tGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGGGGGGGTGGGTGGTTATTGGGGTGCCGGGGGCTGTTTTTTTTTTCGGGGTTTTCGGAAGGTGGGTCGCCCTTGGcc                                                                    >  1:629223/1‑149 (MQ=255)
                                                                                  cgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGAGGGGGTGGGTGTTTATTGGGGTGTCGGGGGGTGTTTTTTTTTTCGGGGTTTCCGAAAACTGGATCGCCATTTGCCTGAc                                                                >  2:588387/1‑150 (MQ=255)
                                                                                  cgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGGCGGGTGGTTATGGGGGTGGCGGGCGGCGGTTTTTTTTTCGGGGTTTTCGGAAACTGGGTTGCCCTTGGCCTggcc                                                               >  1:735093/1‑148 (MQ=255)
                                                                                   gacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGCCGGCCGTTTATGGGGGTGGCGGGGGCTGTTTTTTTTTTCCGGGTTTCCGAAAGCTGGGGCGccct                                                                         >  1:177279/1‑138 (MQ=255)
                                                                                           tCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGgcgc                                                       <  2:194409/150‑1 (MQ=255)
                                                                                           tCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATTGGCGGGTGGCTGTTAATGGGGGTGGCGGGGGCCGTTTTTTTTTTCCGGGGTTTCGAAAACTGGGTCGCCCTTGGCCCGGCATTATGcgcg                                                      >  2:112204/1‑151 (MQ=255)
                                                                                            cTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGCGGGTGGGTGTTTATGGGGGTGGCGGGGGCTGTTTTTTTTTTCCGGGGTTTCGAAAGCTGGGTCGCCCCTTGCCCGGCCTTTGGcgcg                                                      >  2:401947/1‑150 (MQ=255)
                                                                                                gCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGAt                                                 <  2:260699/151‑1 (MQ=255)
                                                                                                   gTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGGGGGTGGCTGTTTATTGGGGTGGCGGGCGCTGTTTTTTTTTTCGGGGTTTTCGAAAGCTGGGTCGCCCTTTGCCCGGCATTAGGCGCGGGGGtt                                                >  2:143105/1‑149 (MQ=255)
                                                                                                    tCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGCGGCTGGCTGTTAATGGGGGTGCCGGGGGCTGTTTTTTTTTTCGGGGTTTTCCAAAGCTGGGTCGCCATTGGCCCGACATTAGGGGGGGGGGTTaa                                              >  2:620717/1‑150 (MQ=255)
                                                                                                         gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGaa                                        <  2:916716/151‑1 (MQ=255)
                                                                                                                   aTTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGcc                              <  2:241439/151‑1 (MQ=255)
                                                                                                                     tCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCg                             <  2:71198/150‑1 (MQ=255)
                                                                                                                         gACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCg                         <  1:239326/150‑1 (MQ=255)
                                                                                                                           cGGCATTATCGGCGGGTGCGTCGGGTATGAGGGGCTGGCTGTTTATTGGGGTGCCGGGCGCTGTTTTTTTTTTCGGGGTTTCCGGAAGCTGGGTCGCCCTTGGCCTGACATTAGGCGCGGGGGTTTAATGGGAGCTGGGGGGCggggg                         >  1:242599/1‑146 (MQ=255)
                                                                                                                                       cggGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAc          <  2:252787/151‑1 (MQ=255)
                                                                                                                                                cGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATaca  >  2:258933/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                         
CGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACA  >  NC_000913/1079321‑1079614

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: