Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,322,868 G→T 9.5% L330M (CTG→ATG)  PP_1155 ← YaiC/YhcK/ArdA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,322,8680GT 9.5% 124.2 / 3.3 42L330M (CTG→ATG) PP_1155YaiC/YhcK/ArdA family protein
Reads supporting (aligned to +/- strand):  ref base G (21/17);  new base T (2/2);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGACGCTCCATGGCT  >  NC_002947/1322731‑1323013
                                                                                                                                         |                                                                                                                                                  
gcgcGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTg                                                                                                                                         >  8:285356/1‑149 (MQ=255)
  gcgAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCtt                                                                                                                                              <  2:128501/142‑1 (MQ=255)
  gcgAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTgagat                                                                                                                                         >  1:128501/1‑143 (MQ=255)
  gcgAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGcc                                                                                                                                       >  7:28679/1‑149 (MQ=255)
             cTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCa                                                                                                                           <  8:28679/149‑1 (MQ=255)
              ttCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGCTGTACATCGccaag                                                                                                                                             <  4:151323/131‑4 (MQ=255)
              ttCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGCTGTACATCGccaag                                                                                                                                             >  3:151323/1‑128 (MQ=255)
              ttCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGCTGTACATCGccaa                                                                                                                                              <  1:22975/130‑3 (MQ=255)
              ttCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGCTGTACATCGccaa                                                                                                                                              >  2:22975/1‑128 (MQ=255)
               tCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGt                                                                                                                         >  6:1710/1‑149 (MQ=255)
                cAGGCAACCTCCCGCCAGGCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAg                                                                                                                          <  7:168984/147‑1 (MQ=255)
                   gCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTc                                                                                                                     >  2:187348/1‑149 (MQ=255)
                   gCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGAAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCAAGTGCTc                                                                                                                     >  2:98958/1‑149 (MQ=255)
                     aaCCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCgg                                                                                                                   <  7:274060/149‑1 (MQ=255)
                       ccTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCgggg                                                                                                                 <  2:47106/149‑1 (MQ=255)
                              ccAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTc                                                                                                            >  2:25499/1‑147 (MQ=255)
                                    agtcagGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTgg                                                                                                    >  3:285777/1‑149 (MQ=255)
                                               aggTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGc                                                                                         <  5:347682/149‑1 (MQ=255)
                                                 gTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGa                                                                                       >  8:208442/1‑149 (MQ=255)
                                                     aaTCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGCCTCGAGGCTGGCCTGGAAGGGCGACAgg                                                                                   >  6:342957/1‑149 (MQ=255)
                                                     aaTCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAg                                                                                    <  1:187348/148‑1 (MQ=255)
                                                       tCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGcc                                                                                 >  7:113849/1‑149 (MQ=255)
                                                              aGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACt                                                                           >  4:134677/1‑148 (MQ=255)
                                                               ggCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCa                                                                         <  5:1710/149‑1 (MQ=255)
                                                               ggCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCa                                                                         <  1:25499/149‑1 (MQ=255)
                                                                 ccAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAgg                                                                       <  5:253032/149‑1 (MQ=255)
                                                                  cAGCCTTAGCGTGGAAGCTTCACCGCGGGCAAAATTGACCTTATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCCCAGCCAGTGCTCGGGTGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAgg                                                                       <  1:98958/149‑1 (MQ=255)
                                                                  cAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGc                                                                      >  6:275253/1‑149 (MQ=255)
                                                                        gAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATgc                                                                >  1:341617/1‑149 (MQ=255)
                                                                               gCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTg                                                         >  6:357395/1‑149 (MQ=255)
                                                                                cAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGc                                                        <  7:208442/149‑1 (MQ=255)
                                                                                               cAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTg                                         >  3:260144/1‑149 (MQ=255)
                                                                                                  aaTTGACGTGATAGCGTCCATCGTGCTTGGCGATGTAGAGCGCGTTGTCTGCCTGTTCGAGCCAGTTCTCGTGTGACTCGAGGCTGGCCTGGAAGGGCGACAGGGCGATACTCATTCTGATGCGCAGGTGCTGTAACTGTGGGTTGCgg                                      <  4:285777/149‑1 (MQ=255)
                                                                                                   aTTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCgg                                      >  8:158818/1‑148 (MQ=255)
                                                                                                    ttGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTa                                    <  5:357395/149‑1 (MQ=255)
                                                                                                            aTCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGaaa                            <  8:113849/149‑1 (MQ=255)
                                                                                                             tCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAAc                           >  1:140029/1‑149 (MQ=255)
                                                                                                                 ggCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGc                       <  7:285356/149‑1 (MQ=255)
                                                                                                                          gCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGa              >  1:80396/1‑149 (MQ=255)
                                                                                                                           cTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGAc             >  1:174823/1‑149 (MQ=255)
                                                                                                                               gCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGACGCTc         >  4:228534/1‑149 (MQ=255)
                                                                                                                                       cAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGACGCTCCATGGCt  <  5:275253/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
GCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGACGCTCCATGGCT  >  NC_002947/1322731‑1323013

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: