Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,605,848 G→T 8.3% L135M (CTG→ATG)  PP_1407 ← sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,605,8480GT 8.3% 151.1 / 4.2 48L135M (CTG→ATG) PP_1407sulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base T (2/2);  total (25/23)
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.73e-01

TGTGACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCGGC  >  NC_002947/1605705‑1605995
                                                                                                                                               |                                                                                                                                                    
tgtgACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCg                                                                                                                                                 >  7:105832/1‑149 (MQ=255)
         gTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCgaga                                                                                                                                        <  3:294802/149‑1 (MQ=255)
           gAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGagaa                                                                                                                                       <  8:105832/148‑1 (MQ=255)
            aaTTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGagaag                                                                                                                                      <  5:97049/148‑1 (MQ=255)
              tttCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGaagaag                                                                                                                                   <  6:43339/149‑1 (MQ=255)
                   cccAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGc                                                                                                                              >  2:296878/1‑149 (MQ=255)
                    ccAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTggtg                                                                                                                             <  8:155559/149‑3 (MQ=255)
                      aGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCCTGATCGGct                                                                                                                                              <  1:40643/130‑2 (MQ=255)
                      aGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCCTGATCGGct                                                                                                                                              >  2:40643/1‑129 (MQ=255)
                        gcgcTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATg                                                                                                                          >  2:26418/1‑148 (MQ=255)
                                    gTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCga                                                                                                                                          >  8:178809/1‑120 (MQ=255)
                                    gTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCga                                                                                                                                          <  7:178809/120‑1 (MQ=255)
                                              cTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCcgc                                                                                                   <  3:76030/149‑1 (MQ=255)
                                                 aGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCgg                                                                                                 >  1:125296/1‑148 (MQ=255)
                                                 aGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCgg                                                                                                 <  3:204276/148‑1 (MQ=255)
                                                   gATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCaa                                                                                              >  2:96724/1‑149 (MQ=255)
                                                    aTGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAg                                                                                                         >  7:165483/1‑137 (MQ=255)
                                                    aTGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAg                                                                                                         <  8:165483/137‑1 (MQ=255)
                                                                cccTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGcc                                                                                 >  8:10544/1‑149 (MQ=255)
                                                                 ccTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGc                                                                                         <  3:76754/140‑1 (MQ=255)
                                                                 ccTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGc                                                                                         >  4:76754/1‑140 (MQ=255)
                                                                 ccTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCg                                                                                >  1:157428/1‑149 (MQ=255)
                                                                  cTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGa                                                                               >  6:49103/1‑149 (MQ=255)
                                                                    cGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGAc                                                                              >  7:194909/1‑148 (MQ=255)
                                                                     gATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGAGCAGGATCGGCCGCGAGAAGAAGCTCGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGAc                                                                              <  4:35875/147‑1 (MQ=255)
                                                                        cTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCa                                                                        >  5:25894/1‑149 (MQ=255)
                                                                        cTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCa                                                                        >  5:30625/1‑149 (MQ=255)
                                                                           gAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACg                                                                      <  1:296878/148‑1 (MQ=255)
                                                                                ccaccaCCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACg                                                                <  7:10544/149‑1 (MQ=255)
                                                                                 caccacCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGa                                                               >  6:309887/1‑149 (MQ=255)
                                                                                  accaccTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGAt                                                              >  3:75297/1‑149 (MQ=255)
                                                                                      ccTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCAcc                                                          <  6:25894/149‑1 (MQ=255)
                                                                                      ccTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCAcc                                                          <  6:30625/149‑1 (MQ=255)
                                                                                                  gCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGTAATCAGCATCACGCCGACC‑AGCACGGTg                                                                   >  7:128074/1‑128 (MQ=255)
                                                                                                  gCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGTAATCAGCATCACGCCGACC‑AGCACGGTg                                                                   <  8:128074/128‑1 (MQ=255)
                                                                                                    ccGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGAt                                            >  8:151000/1‑149 (MQ=255)
                                                                                                    ccGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGAt                                            <  5:309887/149‑1 (MQ=255)
                                                                                                         gATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGAGCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGCCAGTTCGACGATgcgcg                                       >  8:125533/1‑148 (MQ=255)
                                                                                                               gctcaggcCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGt                                 >  1:259098/1‑149 (MQ=255)
                                                                                                                       cGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATg                         <  1:26418/149‑1 (MQ=255)
                                                                                                                                 gAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGcc                      >  1:166413/1‑142 (MQ=255)
                                                                                                                                 gAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGcc                      <  2:166413/142‑1 (MQ=255)
                                                                                                                                    gcggCCGATCCTGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGa                                                        >  5:256137/4‑105 (MQ=255)
                                                                                                                                    gcggCCGATCCTGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGa                                                        <  6:256137/102‑1 (MQ=255)
                                                                                                                                     tAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCAc           <  7:151000/149‑1 (MQ=255)
                                                                                                                                      aGCCGCTCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGAGACGATCACCGACAGATCGACGATCCGCTCCGGGTCCCCCAAGGCCACCGGTGCCACg          >  3:167574/1‑149 (MQ=255)
                                                                                                                                       gCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCAc           <  1:96724/148‑1 (MQ=255)
                                                                                                                                       gCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACgg         >  3:16220/1‑149 (MQ=255)
                                                                                                                                               ggATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATAACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCGGc  >  2:285885/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGTGACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACC‑AGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCGGC  >  NC_002947/1605705‑1605995

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: ‑