Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,080 G→C 6.9% intergenic (‑223/+272) lysO ← / ← aqpZ L‑lysine exporter/water channel AqpZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,0800GC 6.9% 130.4 / 3.5 53intergenic (‑223/+272)lysO/aqpZL‑lysine exporter/water channel AqpZ
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base C (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

TTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCC  >  NC_000913/914938‑915219
                                                                                                                                              |                                                                                                                                           
ttAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATg                                                                                                                                             >  2:401901/1‑143 (MQ=255)
          aTGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGTGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCt                                                                                                                                   <  2:853786/143‑1 (MQ=255)
            gTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGa                                                                                                                                 <  1:72650/143‑1 (MQ=255)
            gTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCGTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGGGGGGTGa                                                                                                                                 >  2:326187/1‑143 (MQ=255)
             tgtagCAAAAGGACGAGTATAGTATTCTGGTTCATTTAAAGATATTGTTGCGCTTTAATTACAAACTGCACCGATGTTTGTGGCGTCAAAATCGCCGAGGCGTTCACTGAAGGCCAGGGCAGCCCACATGGATGTGGGCTGAg                                                                                                                                <  2:333184/139‑1 (MQ=21)
                   aTAAGGAGCAGTAGAGTATTCGTTTTGATTCAAAGAAATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGACGTCGATATCGCCGAGGCGTTCCCTGAAGGCCGGGGCGGCCCAGATGGAGGTGGGGTGAGGGcgcg                                                                                                                          >  2:485337/1‑143 (MQ=25)
                   aTAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGcgcg                                                                                                                          >  2:181401/1‑143 (MQ=255)
                        gAGCAGTAGAGTATTCGTTTTCATTTAAAGCTATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTa                                                                                                                     <  2:368021/143‑1 (MQ=255)
                            tgTAGAGGATTCGTTTTCATTTAGAGATATTCTTGCGCTTTAATTACAAGCTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAgg                                                                                                                 <  2:780547/142‑1 (MQ=255)
                             gTAGAGTATTCGTTTTCATGTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAggg                                                                                                                <  2:488330/143‑1 (MQ=255)
                               agagTATTCGTTTTCATTGAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGAt                                                                                                              <  2:688335/143‑1 (MQ=255)
                                 attATTCGTTTTCATTTAACGTTACTCCTGCGCTTTAACTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGt                                                                                                            <  2:972269/141‑1 (MQ=18)
                                   tATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCGCATGGACGGGGGCTGAGGGCGCGGGTTACAGGGGTGttg                                                                                                          >  1:73198/1‑142 (MQ=25)
                                    aTTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAc                                                                                                         >  2:284868/1‑143 (MQ=255)
                                     ttCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCCGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAcc                                                                                                        >  1:453459/1‑143 (MQ=255)
                                     ttCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTAcc                                                                                                        <  2:377032/143‑1 (MQ=255)
                                        gTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTcg                                                                                                     >  1:780526/1‑143 (MQ=255)
                                         cttgcATTGAAGGATATTCTTGCGCTCTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAAGCGCCGAGGCGTTCCCTCAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTcgc                                                                                                    <  1:403675/139‑1 (MQ=18)
                                           ttCATTTAAAGATATGCTGGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCgcgc                                                                                                  <  2:483253/143‑1 (MQ=255)
                                              aTTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCg                                                                                               <  2:949280/143‑1 (MQ=255)
                                               tttAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCt                                                                                                                                   <  2:347117/106‑1 (MQ=255)
                                               tttAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCt                                                                                                                                   >  1:347117/1‑106 (MQ=255)
                                                 tAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGAcc                                                                                            <  2:631222/143‑1 (MQ=255)
                                                      atatTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTa                                                                                       >  1:132568/1‑143 (MQ=255)
                                                                   tgaaTCACAACCGGCCCCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATa                                                                          <  1:153830/141‑1 (MQ=255)
                                                                        tACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCg                                                                     >  2:103373/1‑143 (MQ=255)
                                                                         aCAAACTGCACCGCCGTTGGTGGCGTCAAAACCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGa                                                                    <  1:116324/143‑1 (MQ=255)
                                                                                    cGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCgcgc                                                         >  2:146799/1‑143 (MQ=39)
                                                                                      aTGTTTGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACCGGGCAGTGGGCTGAGGGCGCGTCTCACAGGGCTGTTACCTCCCGCCCGACCCGGTGGCCGTAAGGGGTGAGCCGAGGGCACCGCGCaa                                                       >  2:787642/1‑142 (MQ=11)
                                                                                         ttGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcgg                                                    <  1:401901/143‑1 (MQ=35)
                                                                                         ttGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcgg                                                    <  1:428782/143‑1 (MQ=35)
                                                                                                 gTCAAAGTCGCGGTGGCGGTCCCTGAAGGCCGGGGCAGCCGACCTGGATGTGGGCTGAGGGCGCGTTTTACAGGGCTGTTACCTCGCGCCCGACCCGGTAGCCGTATGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTg                                            <  2:19219/143‑1 (MQ=11)
                                                                                                 gTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGTGGGCGCGTTATACAGGTATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGCTCAGTCGAGGGCACCGCGCAGCGGCGATTTTg                                            >  2:961705/1‑143 (MQ=11)
                                                                                                  tCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGGTTTTGt                                           >  2:125473/1‑143 (MQ=25)
                                                                                                  tCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGt                                           >  1:338002/1‑143 (MQ=35)
                                                                                                    aaaaGCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTc                                         <  1:422770/143‑1 (MQ=21)
                                                                                                     aaaTCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGATGAGGACTTGAGTCACGGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCg                                        <  1:719535/143‑1 (MQ=11)
                                                                                                          gCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCaga                                   <  2:881048/143‑1 (MQ=32)
                                                                                                          gCCGAGGCGTGCCCTGAGGGCGGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCaga                                   <  2:148816/143‑1 (MQ=12)
                                                                                                              aGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTCTTCGCCAGAGccc                               >  1:712382/1‑143 (MQ=17)
                                                                                                                     cccTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTg                        >  1:928255/1‑143 (MQ=32)
                                                                                                                         ggaggggggggCAGCCGACATGGATGGGGGCTGAGGGGGTGTTTTACAGGGATGTTACCTCGGGCCCGACCCGGTAGCCGTAAGGGATAAGTGGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAgg                    <  2:472213/136‑1 (MQ=11)
                                                                                                                         gAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAgg                    >  2:866029/1‑143 (MQ=32)
                                                                                                                         gAAGGCCGGAGCAGCCCACATGGATGTAGGCTGAGGGCGCGTTTTACAGGGATGGTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATCAGTCGAGGGCGCCGCGCAGCGGCGACGTTGTTCGCGAGAGGCCGGGGGAGCAgg                    >  2:360014/1‑143 (MQ=11)
                                                                                                                            gcggggggCAGCCGAGAGGGGTGTGGGCTGAGGGGGCGTTTGAGAGGGATGTTACCGCGCGCCCGACCCGGTAGCGGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAggggg                 <  1:866091/139‑1 (MQ=11)
                                                                                                                              ccGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTCTGTTCGCCCGAGCCCGGGGGTGCAGGGggcg               >  1:822473/1‑143 (MQ=17)
                                                                                                                                 gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcggcg            >  1:284340/1‑143 (MQ=32)
                                                                                                                                 gggCAGCCCACATGGATGTGGGCTGAGAGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAGGTCGAGGGCACGGCGCAGCGGCGATTTTGTTCGCCAGAGCCGGGGGGTGCAGGAggcggag            >  2:84844/1‑141 (MQ=11)
                                                                                                                                 gggCAGCCCACATCGATGTGGCCTGAGAGCGCGTATTACAGGGACGTTACCTCGCGCGCGCCGCGGTAGCCGTGAGGGGCAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGAGCGAGGagcggag            >  2:215654/1‑141 (MQ=11)
                                                                                                                                   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcggcggc          >  1:844057/1‑143 (MQ=32)
                                                                                                                                   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTGACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcggcggc          <  1:181401/143‑1 (MQ=17)
                                                                                                                                       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCAGGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtg      >  1:263628/1‑142 (MQ=17)
                                                                                                                                       cccACATGGATGTGGGCTGAGCGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCGCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt      >  1:197807/1‑143 (MQ=17)
                                                                                                                                           caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGgcc  >  1:12186/1‑143 (MQ=32)
                                                                                                                                              |                                                                                                                                           
TTAAAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCC  >  NC_000913/914938‑915219

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