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

ATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAA  >  NC_000913/1079323‑1079618
                                                                                                                                                 |                                                                                                                                                                  
aTGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGa                                                                                                                                                              <  1:543174/151‑1 (MQ=255)
  gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGCTATGAgg                                                                                                                                                            >  1:976478/1‑150 (MQ=255)
  gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGc                                                                                                                                                            <  1:548244/151‑1 (MQ=255)
  gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGc                                                                                                                                                            <  1:2792/151‑1 (MQ=255)
     ggTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCgg                                                                                                                                                         <  2:472076/150‑1 (MQ=255)
          cATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTGGGGTCTGAggcgcagg                                                                                                                                                      >  2:233160/1‑142 (MQ=255)
           aTTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTGGGGTATGAgggg                                                                                                                                                         >  2:168251/1‑144 (MQ=255)
               tgtgTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGGgggtgggcggtt                                                                                                                                               >  1:63580/1‑144 (MQ=255)
                   tCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGTCGTGTCGGGTGTGGGGGgctgggtggtta                                                                                                                                             >  1:407328/1‑143 (MQ=255)
                      atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATggg                                                                                                                                      <  1:189711/148‑1 (MQ=255)
                             tGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCggg                                                                                                                             <  2:702812/150‑1 (MQ=255)
                              ggCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGgcg                                                                                                                          <  1:275159/151‑1 (MQ=255)
                              ggCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGgcg                                                                                                                          <  2:639766/151‑1 (MQ=255)
                                 aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGGTGGCTGTTTATGGGGGTGGCggg                                                                                                                             >  1:83644/1‑146 (MQ=255)
                                 aTGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGCGGGGGGTTTATGGGGGTGGCggg                                                                                                                             >  1:597910/1‑146 (MQ=255)
                                    atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGTATGAGGGGCTGGGTGTTTCTGGGGTTGTCggggg                                                                                                                           >  2:946835/1‑143 (MQ=255)
                                    atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTTTGTGCGGGCGGGCGTTTATTGGGGTTGcgg                                                                                                                               >  2:1073555/1‑140 (MQ=255)
                                      attgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGttttt                                                                                                              <  1:252573/151‑1 (MQ=255)
                                        tgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTtttc                                                                                                             <  1:491831/150‑1 (MQ=255)
                                                tttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTCTGAGCGGCGGGGGGCTGATGGGGGTGGCGGGCGGCGGTTTTTTTTTCCggg                                                                                                     >  2:336864/1‑151 (MQ=255)
                                                 ttATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGa                                                                                                    <  2:71300/150‑1 (MQ=255)
                                                  tATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGGGGGTGGGTGTTTATGGGGGTGCCGGGGGCCGGTTTTTTTTTCCgggtt                                                                                                   >  1:490088/1‑149 (MQ=255)
                                                       ccTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGa                                                                                             <  2:53237/151‑1 (MQ=255)
                                                         tGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGGTGGGTGTTAATGGGGGTTGCGGGGGGGGGTTTTTTTTTCGGGGTTTTCGGa                                                                                             >  1:862618/1‑150 (MQ=255)
                                                         tGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATTAGCGGGGGGGGGTTGATTGGGGTTTCGGGGGCCGGTTTTTTTTTCCGGGTTTCCCGa                                                                                             >  2:508743/1‑150 (MQ=255)
                                                           gCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGc                                                                                         <  1:250255/151‑1 (MQ=255)
                                                             gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATTGGGGGGTGGGTGTTTATGGGGGTGGCGGGGGGGGGTTTTTTTTTCGGGGTTTTCGTAAACt                                                                                        >  2:441775/1‑150 (MQ=255)
                                                             gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATTGGGGGGCGGGTGTTTATTGGGGTGGCGGGGGGGGTTTTTTTTTTTCGGGTTTTCCaaaa                                                                                           >  2:662894/1‑148 (MQ=255)
                                                                cAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGa                                                                                     <  1:758345/150‑1 (MQ=255)
                                                                 aaCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATGAGGGGGTGGCTGTTTATGGGGGTGTCGGGGGCGGGTTTTTTTTTCCGGGTTTTCGGaaaatgggtt                                                                                    >  1:1067943/1‑144 (MQ=255)
                                                                   cGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGCGCGTCGGCTATGAGCGGGGGGCTGGTTAAGGGGGTTCCGGGGGCTGTTTTTTTTTTCGGGGATTTCGaaaagag                                                                                       >  2:1027580/1‑141 (MQ=255)
                                                                     aaaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATTAGCGGCTGGGTGTTTATTGGGGTGTCCGGGGCCGGTTTTTTTTTTCGGGGTTTCGAAAGCTGGGTCGcc                                                                               >  2:989832/1‑150 (MQ=255)
                                                                           tttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGCGGGGGGGTGTTTAAGGGGGTGGCGGGCGCTGTTTTTTTTTTCGGGGATTTCGGAAACTGGGTTGCCATTgggc                                                                        >  1:870565/1‑149 (MQ=255)
                                                                              gacgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGa                                                                     <  1:323977/151‑1 (MQ=255)
                                                                                  acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATTAGGGGCCGGGTGGTTATTGGGGTGGCGGGGGGTGTTTTTTTTTTCGGGGTTTTCGAAAGCTTGGGCGCCCCTGGCCTGGCAt                                                                  >  1:440673/1‑150 (MQ=255)
                                                                                   cTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGATATGAGGGGCGGGGTGTTTATGGGGGTGCGCGGCGCGGTTTTTTTTTTCGGGGGTTCCCCAAGCTGGATCGCCCTTTGCCtctca                                                                   >  1:256354/1‑144 (MQ=255)
                                                                                     atatTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGAGGGGGTGGGTGTTTATGGGGGTGGCGGGGGGCGTTTTTTTTTTCGGGGGTTTCGAAAACTGGGTCGCCCTTTGCCTGGCATTatg                                                              >  1:239305/1‑149 (MQ=255)
                                                                                                       gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGaa                                            <  2:199318/151‑1 (MQ=255)
                                                                                                       gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATTGGCGGGCGGGTGTTTATTGGGGTTGCGGGGGCCGTTTTTTTTTTCGGGGGTTTCCGAAGGTGGGTTGCCACTTGGCTGGCATTTGgcggggggggttaactg                                              >  2:777027/1‑135 (MQ=255)
                                                                                                        tCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATGAGGGGGTGGCTGGTTATGGGGGTGGCGGGGGCCGGTTTTTTTTTCGGGGTTTTCGAAAGCTGGGTCGCCCTTGGCCTGGCATTTGGCGCGGGGGTTAAATGGaa                                            >  1:754833/1‑150 (MQ=255)
                                                                                                           tGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCt                                         <  2:84356/150‑1 (MQ=255)
                                                                                                               ccATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtggtgg                                    <  1:978890/151‑1 (MQ=255)
                                                                                                               ccATTCGTGACGGCATTATCGGCGGGGGCGACGGGTATGAGGGGCGGGGGGTTTATGGGGGGGCCGGGGGGTGTGTTTTTTTTCTGGGATTTTGGAAGATCGGGTGCCAGTTTCCCGCCAATAGGCGCGGGGGTTTACTGGAAGGtgg                                       >  1:353881/1‑148 (MQ=255)
                                                                                                                  ttCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGGGGCTGGGTGTTTATGGGGGTGCCGGGCGCTGTTTTTTTTTTCGGGGTTTTCGAAAGCTGGGGCGGCCGTTGGCTGTCATTATGGGCGGGGGTTTACTGGGAGCtggggg                                    >  1:539196/1‑148 (MQ=255)
                                                                                                                    cGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCggg                               <  1:689155/151‑1 (MQ=255)
                                                                                                                          ggCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTg                         <  2:63573/151‑1 (MQ=255)
                                                                                                                                aTCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCgtgt                    <  2:13314/150‑1 (MQ=255)
                                                                                                                                      gggTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATAcc             <  2:178762/151‑1 (MQ=255)
                                                                                                                                          gCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAAt         <  2:791828/151‑1 (MQ=255)
                                                                                                                                            gTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATac       <  2:1003606/151‑1 (MQ=255)
                                                                                                                                              cGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATgca      <  2:664652/150‑1 (MQ=255)
                                                                                                                                              cGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaa     <  1:650698/151‑1 (MQ=255)
                                                                                                                                                 atatGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATAcaacaa  <  2:239243/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                                  
ATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAACAA  >  NC_000913/1079323‑1079618

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: