Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,873,279 G→A 8.3% A114V (GCG→GTG)  aqpZ ← aquaporin Z

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,873,2790GA 8.3% 131.8 / 2.7 49A114V (GCG→GTG) aqpZaquaporin Z
Reads supporting (aligned to +/- strand):  ref base G (19/25);  new base A (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

GTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGCCGGGAAGCGCCCGCCCACCACCAGCC  >  NC_002947/4873133‑4873411
                                                                                                                                                  |                                                                                                                                    
gTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCa                                                                                                                                    <  4:178074/149‑1 (MQ=255)
gTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCa                                                                                                                                    <  6:253303/149‑1 (MQ=255)
       aCCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCt                                                                                                                             <  7:236782/149‑1 (MQ=255)
                     gggCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTg                                                                                                               <  2:298065/149‑1 (MQ=255)
                      ggCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGc                                                                                                              <  3:111914/149‑1 (MQ=255)
                             gtcgCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCg                                                                                                                     <  4:193873/135‑1 (MQ=255)
                             gtcgCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCg                                                                                                                     >  3:193873/1‑135 (MQ=255)
                             gtcgCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGc                                                                                                       >  2:324923/1‑149 (MQ=255)
                                 cccccATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGAt                                                                                                    >  5:359035/1‑148 (MQ=255)
                                 cccccATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATg                                                                                                   >  8:68221/1‑149 (MQ=255)
                                   cccATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAg                                                                                                 >  5:181740/1‑149 (MQ=255)
                                      atgatgACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTa                                                                                              <  1:324923/149‑1 (MQ=255)
                                          tgtccaccaGGAACATCGCCGTCATCACCACTTCTCTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGAt                                                                                           <  4:177197/145‑1 (MQ=255)
                                              caccaGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAg                                                                                                 >  1:364142/1‑138 (MQ=255)
                                              caccaGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAg                                                                                                 <  2:364142/138‑1 (MQ=255)
                                                    gAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCgg                                                                                                                    <  2:174355/113‑1 (MQ=255)
                                                        aTCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCCCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAg                                                                            >  8:151116/1‑149 (MQ=255)
                                                                      cacTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCa                                                              >  7:140643/1‑149 (MQ=255)
                                                                       acTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATca                                                                                         >  7:81536/1‑121 (MQ=255)
                                                                       acTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATca                                                                                         <  8:81536/121‑1 (MQ=255)
                                                                         ttCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCt                                                           >  4:244048/1‑149 (MQ=255)
                                                                         ttCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCt                                                           >  4:16559/1‑149 (MQ=255)
                                                                           cACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCt                                                           <  3:244048/147‑1 (MQ=255)
                                                                            aCTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCAGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTgg                                                         <  7:236883/148‑1 (MQ=255)
                                                                              tGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGa                                                                 <  1:69834/138‑1 (MQ=255)
                                                                              tGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGa                                                                 >  2:69834/1‑138 (MQ=255)
                                                                              tGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGc                                                       >  3:342201/1‑148 (MQ=255)
                                                                                aCAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGa                                                     <  3:153703/148‑1 (MQ=255)
                                                                                  aaaGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTggg                                                        <  8:284930/143‑1 (MQ=255)
                                                                                  aaaGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTggg                                                        >  7:284930/1‑143 (MQ=255)
                                                                                  aaaGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCACAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATca                                                  >  7:7378/1‑149 (MQ=255)
                                                                                  aaaGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCACAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATca                                                  >  7:2168/1‑149 (MQ=255)
                                                                                     gCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATcaca                                                >  6:75283/1‑148 (MQ=255)
                                                                                        tGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAg                                            >  3:107376/1‑149 (MQ=255)
                                                                                                        ccACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCTATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGc                            <  7:151116/149‑1 (MQ=255)
                                                                                                        ccACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGc                                                       >  3:290672/1‑122 (MQ=255)
                                                                                                        ccACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGc                                                       <  4:290672/122‑1 (MQ=255)
                                                                                                         cACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGcc                           <  3:16559/149‑1 (MQ=255)
                                                                                                         cACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGcc                           <  6:181740/149‑1 (MQ=255)
                                                                                                         cACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGcc                           >  6:139101/1‑149 (MQ=255)
                                                                                                              ggCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGACCACATAAGGCAGCAGCTCTTTTGCCGTGaa                      <  5:280370/149‑1 (MQ=255)
                                                                                                                 gAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGCCGGGAAgcg                   <  4:107376/149‑1 (MQ=255)
                                                                                                                    ttgTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGCCGGGAAGcgccc                <  6:359035/147‑1 (MQ=255)
                                                                                                                      gTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGCCGGGAAGcgcccg               <  7:68221/148‑1 (MQ=255)
                                                                                                                                 cGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGCCGGGAAGCGCCCGCCCACCACCAGc   >  4:241758/1‑149 (MQ=255)
                                                                                                                                  gTTTGAAGCCAACCCCACAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGCCGGGAAGCGCCCGCCCACCACCAGcc  <  8:2168/149‑1 (MQ=255)
                                                                                                                                  gTTTGAAGCCAACCCCACAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGCCGGGAAGCGCCCGCCCACCACCAGcc  <  8:7378/149‑1 (MQ=255)
                                                                                                                                      gAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCt                                                           >  7:219107/1‑88 (MQ=255)
                                                                                                                                      gAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCt                                                           <  8:219107/88‑1 (MQ=255)
                                                                                                                                                  gCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCa                                                                             <  6:40900/58‑1 (MQ=255)
                                                                                                                                                  gCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCa                                                                             >  5:40900/1‑58 (MQ=255)
                                                                                                                                                  |                                                                                                                                    
GTGCAAAACCTGCCGGTGCGCGGGCGTCGGTCGCCCCCATGATGACCACCAGGAACATCGCCGTCATCACCACTTCACTGACAAAGCCTGCACCCAGCGTATAGCCACCGGGCGAGTGGTCGGCGTAGCCGTTTGAAGCCAACCCCGCAGACAGCTCGAAACCGGCCTTGCCGCTGGCGATGAGGTAGATCACACCGGCCGCCAGGATGGCACCGATCACCTGGGCGATCACATAAGGCAGCAGCTCTTTTGCCGGGAAGCGCCCGCCCACCACCAGCC  >  NC_002947/4873133‑4873411

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑