Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 915,081 G→C 10.5% intergenic (‑224/+271) 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,0810GC10.5% 112.9 / 4.8 45intergenic (‑224/+271)lysO/aqpZL‑lysine exporter/water channel AqpZ
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base C (2/3);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCG  >  NC_000913/914941‑915220
                                                                                                                                            |                                                                                                                                           
aaaGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGAGGCGTCAAAATCGCCGACGCGTTCCCTGAAGGCCGGGGAAGCCAACATGGAt                                                                                                                                           >  2:417819/1‑143 (MQ=255)
    gAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTgg                                                                                                                                       <  1:225001/143‑1 (MQ=255)
       aTGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCt                                                                                                                                    <  2:604677/143‑1 (MQ=255)
        tGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATCCATGTGGGCTg                                                                                                                                   <  1:587690/143‑1 (MQ=255)
              gCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGgcg                                                                                                                             >  1:701365/1‑143 (MQ=255)
              gCAGAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATGTTCTTGCGCTTTAATTACAATCTGCACCGATGTTGGTGGCGTCAAAATCGCAGAGGCGTTCCCTGAAGGCTGGGGCAGCCCATATGGATGTGGGCTGAGGgcg                                                                                                                             <  2:551353/143‑1 (MQ=37)
               cATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAg                                                                                                                                 >  2:732540/1‑138 (MQ=255)
               cATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAg                                                                                                                                 <  1:732540/138‑1 (MQ=255)
                        cAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAg                                                                                                                   >  2:448211/1‑143 (MQ=255)
                         aGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAgg                                                                                                                  >  1:850422/1‑143 (MQ=255)
                           tAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGGGCTGCTTTTATGCTGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGa                                                                                                                >  1:1043976/1‑143 (MQ=25)
                           tAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGa                                                                                                                >  2:273486/1‑143 (MQ=255)
                           tAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGa                                                                                                                >  2:203921/1‑143 (MQ=255)
                            agagTATTCGTTTTCATTTAAAGAAATTCTTGCGCTTGAATTACAAACTGCACCGATGTTGGTGAGGTAATATGCACCAAGGCGTTCCATGAACGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGAt                                                                                                               <  1:244828/143‑1 (MQ=21)
                              agTATTCGTTTGCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGt                                                                                                             <  1:967013/143‑1 (MQ=255)
                                        ttCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGCATGTCGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCgcgc                                                                                                   <  1:400102/143‑1 (MQ=255)
                                          cATTTAAAGCTATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGccc                                                                                                 <  1:331926/143‑1 (MQ=255)
                                               aaaGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGAccc                                                                                            >  2:593616/1‑143 (MQ=255)
                                                 aGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCgg                                                                                          >  1:700183/1‑143 (MQ=255)
                                                   gtgttCTTGCGCTGGAATTACAAACTGCACCGATGTTGGTGGCGTCAAGATCGCCGAGGCGTTCCCTGAGGGCCGGGGCAGCCGACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTa                                                                                        <  2:451530/140‑1 (MQ=17)
                                                   atatTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTAGAGGGATGTTACGTCGCGCCCGACCCGGTa                                                                                        >  1:262033/1‑143 (MQ=255)
                                                    tatTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTACCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAg                                                                                       <  1:1045857/143‑1 (MQ=255)
                                                        cgtGCGCTTTAATTACACACCGCACCGCTGGGGGGGGCGTCAAAATCGCCGAGGCGTGCCCGGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGt                                                                                   <  2:215300/141‑1 (MQ=14)
                                                          tgccgtttCATTCCAGAGTGCACCGGTGGTGGTGGCGTCAAACTCGCCGAGGGGTCCCCTGAAGGCCGGGGCAGCCCACAGGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTaa                                                                                 <  2:831358/138‑1 (MQ=11)
                                                          tGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTaa                                                                                 <  2:845538/143‑1 (MQ=255)
                                                          tGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTaa                                                                                 <  2:547917/143‑1 (MQ=255)
                                                                  aaTTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAg                                                                         >  2:315137/1‑143 (MQ=255)
                                                                  aaTTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGGAGCCGTAAGGGATAAg                                                                         >  2:201998/1‑143 (MQ=255)
                                                                        aaccTGCACCGATGTTGGGGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACTCCGTAGCCGTAAGGGATAAGTCGAgg                                                                   <  2:922632/140‑1 (MQ=255)
                                                                          aCTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGc                                                                 >  1:180841/1‑143 (MQ=255)
                                                                           cTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCACTGAAGGCGGGGGGAGAGCACATGGATGTGGGAATGGAGTGGGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAggggg                                                                >  1:489752/1‑141 (MQ=14)
                                                                                 cGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCgcgc                                                          <  1:203921/143‑1 (MQ=40)
                                                                                   aTGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAg                                                        >  2:488847/1‑143 (MQ=38)
                                                                                        ggtggCGTCAAAATCGCCGAGGCGGTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAgcggcg                                                   <  1:519816/143‑1 (MQ=25)
                                                                                          tggCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGAGGTGGGCTGAGGGCGCGTTTGACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATGCGTCGAGGGCACCGCGCAgcggcggt                                                 >  1:787700/1‑141 (MQ=255)
                                                                                          tggCGTCAAAATCGCCGAGGCGTTCCCGGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAt                                                 >  1:48447/1‑143 (MQ=25)
                                                                                                cacaaTCGGCGAGGCGTTCCCTGCAGGCCGGGGCAGCCATCATGGATGTGGGCTGAGGGCGCGTGTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGtt                                           <  2:479589/140‑1 (MQ=11)
                                                                                                   aaTCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGc                                        >  1:272617/1‑143 (MQ=34)
                                                                                                       gCCGAGGCGTCCCCGGAAGGCCGGGGCAGCCGACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCaga                                    <  1:210409/143‑1 (MQ=12)
                                                                                                            ggCGTTCGCTGAAGGCCGGGGCAGCCCACATGGAGGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCg                               <  1:634313/143‑1 (MQ=17)
                                                                                                               gTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCCCCAGAGccccgcc                            >  1:857242/1‑139 (MQ=17)
                                                                                                                             ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTCCCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGGTAGCCGCACGCCGGACGTGCAGGGggcggc              >  2:395680/1‑143 (MQ=11)
                                                                                                                               ggcggcCCACATGGATGGGGCCTGAGGGCGGGTAGTACAGGGACGTTATATCGCCAAACGTCTGGTGTCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGggcggcgg            <  1:486426/143‑1 (MQ=2)
                                                                                                                                  aGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCCGCGGCGCTTTTGTGCGCCAGAGCCCGGGGGGGCGGGGGgcggcggcgc         >  2:840027/1‑142 (MQ=11)
                                                                                                                                  aGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGACGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGACAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGCGGGCGGCGGCGa         >  2:101643/1‑143 (MQ=11)
                                                                                                                                      cgcagggATGGGGGCTGAGGGCGCGTTTTACAGGGATGGTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg     <  2:448548/138‑1 (MQ=12)
                                                                                                                                         atggatggGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTTCAGGGGGCGGCGGCGATTGgccg  <  2:642050/143‑1 (MQ=17)
                                                                                                                                            |                                                                                                                                           
AAAGGAAATGTTATGCATAAGGAGCAGTAGAGTATTCGTTTTCATTTAAAGATATTCTTGCGCTTTAATTACAAACTGCACCGATGTTGGTGGCGTCAAAATCGCCGAGGCGTTCCCTGAAGGCCGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGGCCG  >  NC_000913/914941‑915220

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