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

ATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079323‑1079616
                                                                                                                                                 |                                                                                                                                                        
aTGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGa                                                                                                                                                     <  1:995685/151‑1 (MQ=255)
 tGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGTGGTGGata                                                                                                                                                        >  2:270557/1‑147 (MQ=255)
  gTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGc                                                                                                                                                   <  1:388689/151‑1 (MQ=255)
    tGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTGGGGGATGAGCgg                                                                                                                                                 >  2:278106/1‑151 (MQ=255)
                     tataTCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGGTTAGGGGGTGGCGGGTAATgggg                                                                                                                              >  1:699763/1‑150 (MQ=255)
                     tataTCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTAATAGCGGCTGGCTGTTTATtggg                                                                                                                              >  2:769501/1‑150 (MQ=255)
                      atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGtt                                                                                                                            <  2:738830/151‑1 (MQ=255)
                      atatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGAGGGGGTGGGTGTTAATtggg                                                                                                                              >  2:364575/1‑149 (MQ=255)
                       tatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGGGCGGGGGGGTGGTTATgggg                                                                                                                              >  1:867436/1‑148 (MQ=255)
                       tatCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTTTGAGGgggtggggggttattggggtg                                                                                                                              >  2:559030/1‑136 (MQ=255)
                          cTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGGATGAGGGGGTGGGTGTTTATgggg                                                                                                                              >  1:956595/1‑145 (MQ=255)
                           tttGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGGGGGTGGGTGTTTATGGGGGTGGCgg                                                                                                                       >  2:1058735/1‑151 (MQ=255)
                              ggCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGgc                                                                                                                     <  2:679360/150‑1 (MQ=255)
                                   gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGt                                                                                                               <  1:164066/150‑1 (MQ=255)
                                      attgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGttttt                                                                                                           <  1:302519/151‑1 (MQ=255)
                                         gattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTtttctt                                                                                                        <  2:390650/151‑1 (MQ=255)
                                          attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGGGGGGGTGGGGGTTTTTGGGGGTGCCGGGGGGTGTtttttttt                                                                                                       >  2:1010333/1‑151 (MQ=255)
                                          attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGTGGGTATTTGCGGCTGGGTGTTTATTGGGTTGGCCGGGGCCGTtttttttt                                                                                                       >  2:249028/1‑151 (MQ=255)
                                           ttgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGTGGGGATGAGCGGGTGGGTGTTTATTGGGGTGGCGGGGGGCGGTtttttttt                                                                                                       >  2:445068/1‑151 (MQ=255)
                                              ggTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCgg                                                                                                   <  1:1062350/151‑1 (MQ=255)
                                               gTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGGGGGTGGCTGTTTATGGGGGTGCCGGGGGCTGGTTTTTTTTTCggg                                                                                                  >  2:201586/1‑151 (MQ=255)
                                                   aTCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGAtt                                                                                               <  1:773188/150‑1 (MQ=255)
                                                      gCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCg                                                                                           <  1:674070/151‑1 (MQ=255)
                                                             gATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATTAGCGGGTGGGCGGTAATGGGGGTGGCGGGGGGGGGTTTTTTTTTCGGGGGTTTCCCaa                                                                                         >  1:568781/1‑147 (MQ=255)
                                                              aTCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATTGGCGGGCGGGGGGTTTTGGGGGTGGCGGGGGCTGTTTTTTTTTTCCGGGGTTTccca                                                                                          >  2:622074/1‑143 (MQ=255)
                                                                cAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGAt                                                                                 <  1:117491/151‑1 (MQ=255)
                                                                 aaCGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATc                                                                                <  1:929100/151‑1 (MQ=255)
                                                                      aaaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCAt                                                                           <  1:64180/151‑1 (MQ=255)
                                                                       aaaCTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCAtt                                                                          <  1:412874/151‑1 (MQ=255)
                                                                           tttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGTATGAGCGGCTGGGTGTTAATTGGGGTGGCGGGGGGTGTTTTTTTTTTCGGGGTTTTCGCAAACTGGGTCGCCAtt                                                                          >  2:767109/1‑147 (MQ=255)
                                                                            ttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGGATTAGGGGGTGGGTGTTTATTGGGTTGGCGGGGGCGGTTTTTTTTTTCGGGGTTTTCCAAAGGTGGGTCGCCAtttggc                                                                      >  2:571657/1‑146 (MQ=255)
                                                                            ttGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGAGGGGCTGGCTGGTTATGGGGGTGGCGGGCGCCGTTTTTTTTTTCGGGGTTTTCGCAAGCTGGGGCGccccttg                                                                        >  1:727567/1‑143 (MQ=255)
                                                                              gacgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTTTTAGGGGGCGGGTGTTTATGGGGGTGCCGGGGGGCGTTTTTTTTTTCGGGGTTTTCCAAAACTGGagtggcctttgcctct                                                                   >  1:809324/1‑136 (MQ=255)
                                                                                cgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACa                                                                 <  2:96915/151‑1 (MQ=255)
                                                                                 gacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACAt                                                                <  1:253790/151‑1 (MQ=255)
                                                                                  acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGAGGGGCTGGCGGGTTAAGGGGGTGCCGGGGGCTGGTTTTTTTTTCGCGGTTTTTGAAAGATGGGTCGCCCtttggctctcatt                                                               >  2:578189/1‑140 (MQ=255)
                                                                                     atatTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAg                                                             <  2:288765/150‑1 (MQ=255)
                                                                                       atTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGGGGGTATGGGAGGGCGGCTGTTTATTGGGGTGTCGGGGGGTGGTTTTTTTTTTGGGGGTTTCGCAAGCTGGGTCGCCATTCGCCCGACACTAGgcg                                                          >  2:384741/1‑151 (MQ=255)
                                                                                         tCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGgcgc                                                         <  2:215801/150‑1 (MQ=255)
                                                                                                ggTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGATTAGGGGGGGGGTGGTTATTGGGGTGGCCGGGGGTGGTTTTTTTTTTTGGGTTTTCGCAAAGTGGGTCGCCCTTTGCCTTGCATTTGgcggggggggt                                                  >  1:889507/1‑142 (MQ=255)
                                                                                                 gTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTTTGAGCGGGTGGGTGTTAATGGGGTTGCCGGGCGCCGTTTTTTTTTTCCGGGTTTCCGCAAGATGGGGCGCCCCTGGCCCGGCATTAggggggcgg                                                     >  2:147453/1‑139 (MQ=255)
                                                                                                   cGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACt                                              <  1:95146/151‑1 (MQ=255)
                                                                                                      aGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTgg                                            <  1:610787/150‑1 (MQ=255)
                                                                                                       gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGaa                                          <  1:1055107/151‑1 (MQ=255)
                                                                                                       gTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGaa                                          <  2:613106/151‑1 (MQ=255)
                                                                                                           tGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCt                                       <  1:363987/150‑1 (MQ=255)
                                                                                                               ccATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCtggtgg                                  <  1:416339/151‑1 (MQ=255)
                                                                                                                cATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTTTGAGCGGCTGGGTGGTAATGGGGGTGGCCGGGGGTGTTTTTTTTTTTCGGGGTTTCGGAAACTGGGGAGCCATTGGCCCGACATTAGGGGCGGGGGTTTACTTGAAGGtggtggg                                 >  2:445566/1‑150 (MQ=255)
                                                                                                                   tCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCg                               <  1:828025/150‑1 (MQ=255)
                                                                                                                                  cggcggGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCa               <  1:1062447/151‑1 (MQ=255)
                                                                                                                                    gcggGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATa             <  1:834961/151‑1 (MQ=255)
                                                                                                                                          gCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAAt       >  2:528253/1‑151 (MQ=255)
                                                                                                                                           cGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAAt       >  1:868228/1‑150 (MQ=255)
                                                                                                                                               ggATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaac  >  2:161267/1‑151 (MQ=255)
                                                                                                                                                gATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaac  >  1:394760/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                        
ATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGAGCGGCTGGCTGTTAATGGGGTTGCCGGGCGCTGTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079323‑1079616

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