Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,684,468 G→A 8.0% A349A (GCG→GCA)  PP_1478 → xenobiotic reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,684,4680GA 8.0% 143.6 / 2.6 50A349A (GCG→GCA) PP_1478xenobiotic reductase
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base A (2/2);  total (27/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.82e-01

CATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTGAGGATGTCACTTACCCATTGGTTAATGCTCTTTTCCGCCAGTTCGGCTTTA  >  NC_002947/1684321‑1684611
                                                                                                                                                   |                                                                                                                                               
cATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATTATCTGGGTGGCAACAAGGCCACTTGGCCGgcg                                                                                                                                                 >  2:319482/1‑148 (MQ=255)
     tCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTAc                                                                                                                                           <  5:372402/149‑1 (MQ=255)
                 cccAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGa                                                                                                                                 <  6:255606/147‑1 (MQ=255)
                       ccGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCg                                                                                                                                  <  4:273104/140‑1 (MQ=255)
                       ccGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCg                                                                                                                                  >  3:273104/1‑140 (MQ=255)
                             acaTTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAAt                                                                                                                   >  7:48404/1‑149 (MQ=255)
                               aTTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATcc                                                                                                                 >  4:42226/1‑149 (MQ=255)
                                  cTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGa                                                                                                               <  4:193741/148‑1 (MQ=255)
                                     cGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATTATCTGGGTGGCACCAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATa                                                                                                             <  1:319482/147‑1 (MQ=255)
                                           gggCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTc                                                                                                     <  3:42226/149‑1 (MQ=255)
                                                       gTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGa                                                                                         >  3:143591/1‑149 (MQ=255)
                                                         gATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGc                                                                                                                                   >  6:13528/1‑105 (MQ=255)
                                                         gATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGc                                                                                                                                   <  5:13528/105‑1 (MQ=255)
                                                          aTCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGa                                                                                                        >  1:110271/1‑131 (MQ=255)
                                                          aTCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGa                                                                                                        <  2:110271/131‑1 (MQ=255)
                                                          aTCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGAtt                                                                                       >  2:99501/1‑148 (MQ=255)
                                                               ccTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAg                                                                                    <  8:48404/146‑1 (MQ=255)
                                                               ccTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTa                                                                                 >  3:188104/1‑149 (MQ=255)
                                                                 tGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACt                                                                               >  5:8708/1‑149 (MQ=255)
                                                                 tGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGACACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACt                                                                               >  7:60948/1‑149 (MQ=255)
                                                                  ggCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGCCCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTg                                                                              >  2:230433/1‑149 (MQ=255)
                                                                               ctgcGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGa                                                                   <  1:99501/147‑1 (MQ=255)
                                                                               ctgcGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAAt                                                                 >  8:155812/1‑149 (MQ=255)
                                                                                 gcgcGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCACAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGc                                                               >  5:242391/1‑149 (MQ=255)
                                                                                 gcgcGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCACAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGc                                                               >  5:248128/1‑149 (MQ=255)
                                                                                  cgcgATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCg                                                              >  1:9463/1‑149 (MQ=255)
                                                                                       tCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTg                                                         <  5:391795/149‑1 (MQ=255)
                                                                                         cTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCACAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGc                                                       <  6:248128/149‑1 (MQ=255)
                                                                                         cTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCACAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGc                                                       <  6:242391/149‑1 (MQ=255)
                                                                                          ttACTGTCCGCTGCATGCTTATTATGATCTGTGTTGCAACAAGGCCACTTGGCCTGCGCAGTACCAAAGGGCTACCATCCGCGCCAATCCGATACATGAGTCGGACTTGCgg                                                                                           <  6:108060/112‑1 (MQ=255)
                                                                                          ttACTGGCCGCTGCATGCCGATGATGATCTGGTTGGCAAAAATGCCACTTGGCCGGCGCAGTACCAAAGGGAGACCAGCCGCGCCAATCCGATACATGAGTCGAACTTGctg                                                                                           >  5:108060/1‑110 (MQ=255)
                                                                                          ttACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCt                                                      >  5:255788/1‑149 (MQ=255)
                                                                                            aCTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTc                                                                         >  5:251732/1‑128 (MQ=255)
                                                                                            aCTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTc                                                                         <  6:251732/128‑1 (MQ=255)
                                                                                               ggCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGTGAATAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTGAGGa                                                 <  7:310494/149‑1 (MQ=255)
                                                                                                        aTGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGGTGAATGCGGCGGGGCTGAGGATGTCACTTa                                        >  6:420820/1‑149 (MQ=255)
                                                                                                          gCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGc                                                               >  2:340382/1‑124 (MQ=255)
                                                                                                          gCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGc                                                               <  1:340382/124‑1 (MQ=255)
                                                                                                             gatgatgatCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAATGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTg                                                     <  1:359925/131‑1 (MQ=255)
                                                                                                             gatgatgatCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTg                                                     >  2:359925/1‑131 (MQ=255)
                                                                                                                   gatCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACAt                                                                                                          <  7:108557/72‑1 (MQ=255)
                                                                                                                   gatCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACAt                                                                                                          >  8:108557/1‑72 (MQ=255)
                                                                                                                         ggTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGcagta                                                           >  8:167577/1‑109 (MQ=255)
                                                                                                                         ggTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGcagta                                                           <  7:167577/113‑5 (MQ=255)
                                                                                                                         ggTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGTGGCCTGGCTGAGGATGTCACTTACCTATTGGTTAATg                          <  8:60948/146‑1 (MQ=255)
                                                                                                                            ggCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTGAGGATGTCACTTACCCATTGGTTAATGCTCttt                    >  7:68767/1‑149 (MQ=255)
                                                                                                                              caacaaGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGCATGCGGCCTGGCGGAGGATGTCCCTTACCCATGGGTTAATGCTCTTGTc                  >  8:282189/1‑149 (MQ=255)
                                                                                                                              caacaaGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTGAGGATGTCACTTACCCATTGGTTAATGCTCTTTTc                  <  2:9463/149‑1 (MQ=255)
                                                                                                                                          ttGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTGAGGATGTCACTTACCCATTGGTTAATGCTCTTTTCCGCCAGTTCGGc      >  1:287282/1‑149 (MQ=255)
                                                                                                                                                ggCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTGAGGATGTCACTTACCCATTGGTTAATGCTCTTTTCCGCCAGTTCGGCTTTa  >  6:172045/1‑147 (MQ=255)
                                                                                                                                                ggCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTGAGGATGTCACTTACCCATTGGTTAATGCTCTTTTCCGCCAGTTCGGCTTTa  <  5:172045/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CATGATCACTGAACCTGCCCAGGCCGAGCACATTCTGCGTACCGGGCAGGCCGACGTGATCTTCCTGGCTCGGGAGCTGCTGCGCGATCCTTACTGGCCGCTGCATGCCGATGATGATCTGGGTGGCAACAAGGCCACTTGGCCGGCGCAGTACCAAAGGGCGACCAGCCGCGCCAATCCGATACATGAGTCGGACTTGCGGGATTAAGCTACTGCATTCAGCTGAATGCGGCCTGGCTGAGGATGTCACTTACCCATTGGTTAATGCTCTTTTCCGCCAGTTCGGCTTTA  >  NC_002947/1684321‑1684611

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