Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,176,690 G→A 10.5% W664* (TGG→TAG)  PP_2787 → transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,176,6900GA10.5% 100.0 / 3.0 38W664* (TGG→TAG) PP_2787transporter
Reads supporting (aligned to +/- strand):  ref base G (14/20);  new base A (2/2);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

CACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCT  >  NC_002947/3176548‑3176836
                                                                                                                                              |                                                                                                                                                  
cacaTCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCtggctggc                                                                                                                                              >  6:147353/1‑149 (MQ=255)
 acaTCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGcc                                                                                                                                             <  7:320859/149‑1 (MQ=255)
      atgatgTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATcaa                                                                                                                                        >  6:299808/1‑149 (MQ=255)
      atgatgTTCGAGAACGGTGGCGGCACTGAATTGCTGATCGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACGAGGCGATGCTCTATAAACTAGACGCCTGGCTGGCCATcaa                                                                                                                                        <  1:345761/149‑1 (MQ=255)
        gatgTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATcaac                                                                                                                                       >  2:86754/1‑148 (MQ=255)
                       tGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGc                                                                                                                       <  2:282820/149‑1 (MQ=255)
                        ggcggcACTGACTTGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCa                                                                                                                      <  1:86754/149‑1 (MQ=255)
                              aCTGACTTGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCAc                                                                                                                >  6:191336/1‑149 (MQ=255)
                                    ttGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGc                                                                                                          <  2:339653/149‑1 (MQ=255)
                                         gAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGcc                                                                                                                                             <  3:278497/109‑1 (MQ=255)
                                         gAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGcc                                                                                                                                             >  4:278497/1‑109 (MQ=255)
                                             gTGCTCAATGCCGAGTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCAATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCCGTACGcc                                                                                                 >  8:209368/1‑149 (MQ=255)
                                              tGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCCCCCATGCCGGTACGCCg                                                                                                >  8:5784/1‑149 (MQ=255)
                                                   aaTGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACat                                                                                           <  8:187875/149‑1 (MQ=255)
                                                     tGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACatat                                                                                         <  5:147353/149‑1 (MQ=255)
                                                               tagAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCATGCGATGCTCGATAAACTAGACGCCTAGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGgt                                                                               >  4:353412/3‑149 (MQ=255)
                                                               tCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTAGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGgt                                                                               >  4:347903/1‑149 (MQ=255)
                                                                 gAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGgt                                                                               <  5:191336/147‑1 (MQ=255)
                                                                       ggTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAgg                                                                       <  8:48126/149‑1 (MQ=255)
                                                                       ggTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAgg                                                                       <  6:132313/149‑1 (MQ=255)
                                                                             gCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCg                                                                                   <  3:168538/131‑1 (MQ=255)
                                                                             gCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCg                                                                                   >  4:168538/1‑131 (MQ=255)
                                                                                 ttCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCa                                                             <  3:283414/149‑1 (MQ=255)
                                                                                   cACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAg                                                            <  3:49076/148‑1 (MQ=255)
                                                                                       aTGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTAGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCa                                                       <  3:347903/149‑1 (MQ=255)
                                                                                       aTGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTAGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCa                                                       <  3:353412/149‑1 (MQ=255)
                                                                                               ctctGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCt                                                 >  8:173277/1‑147 (MQ=255)
                                                                                               ctctGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGa                                               <  4:159222/149‑1 (MQ=255)
                                                                                                      gcgggcgATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGt                                                                                                      <  2:243767/87‑1 (MQ=255)
                                                                                                        gggcgATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGc                                      >  1:304393/1‑149 (MQ=255)
                                                                                                          gcgATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTg                                    <  5:299808/149‑1 (MQ=255)
                                                                                                                            tCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGa                                               <  3:130087/120‑1 (MQ=255)
                                                                                                                            tCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGa                                               >  4:130087/1‑120 (MQ=255)
                                                                                                                             cGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTtc                                          >  6:15098/1‑124 (MQ=255)
                                                                                                                             cGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTtc                                          <  5:15098/124‑1 (MQ=255)
                                                                                                                              gATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCtt                >  6:305508/1‑149 (MQ=255)
                                                                                                                              gATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCtt                >  4:368754/1‑149 (MQ=255)
                                                                                                                                          gCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGtt    >  3:202717/1‑149 (MQ=255)
                                                                                                                                           ccTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTc   <  3:368754/149‑1 (MQ=255)
                                                                                                                                            ctggctggCCATCAACAAACCGGCCAACCTGCAGGTCACACATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCt  >  7:229880/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CACATCATGATGTTCGAGAACGGTGGCGGCACTGACTTGCAGAACGTGCTCAATGCCGATTACTCGAAGGCGGTGACGCTGTTCACCATGACCTCCTCTGTGGCGGGCGATTACCAGGCGATGCTCGATAAACTAGACGCCTGGCTGGCCATCAACAAACCGGCCAACCTGCAGGTCACCCATGCCGGTACGCCGTACATATGGACCGGGGTGTTGCAGGAAATCACCCAGGGTCAGGTGCTGAGCTTCTCGCTGGCGTTGCTGGCCGTGACCTTGATGATGATGTTCT  >  NC_002947/3176548‑3176836

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