Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,322,861 A→C 9.3% I332S (ATC→AGC)  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,8610AC 9.3% 123.7 / 4.4 43I332S (ATC→AGC) PP_1155YaiC/YhcK/ArdA family protein
Reads supporting (aligned to +/- strand):  ref base A (21/18);  new base C (2/2);  total (23/20)
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.81e-01

GCGGTCCTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGACGCTC  >  NC_002947/1322717‑1323006
                                                                                                                                                |                                                                                                                                                  
gCGGTCCTGGCGATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTa                                                                                                                                                <  7:95588/149‑1 (MQ=255)
gCGGTCCTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGt                                                                                                                                                 <  4:210953/148‑1 (MQ=255)
              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                                                                                                                                       >  2:22975/1‑128 (MQ=255)
                            ttCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGCTGTACATCGccaa                                                                                                                                       <  1:22975/130‑3 (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                                                                  <  1:25499/149‑1 (MQ=255)
                                                                             ggCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCa                                                                  <  5:1710/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)
                                                                                                                                                |                                                                                                                                                  
GCGGTCCTGGCCATGCGCGAGCGGCTCCTTCAGGCAACCTCCCGCCAGTCAGTCAGGAAGGAGGTTCAATCAGGATAGGCCAGCCTGAGCGTGGCAGCTTCACCGCGGGCAAAATTGACCTGATCGCGGCCAGCGTGCTTGGCGATGTACAGCGCCTTGTCTGCCTGTTCGAGCCAGTGCTCGGGGGACTCGAGGCTGGCCTGGAAGGGCGACAGGCCGATACTCAGGCTGATGCGCAGGTGCGGTAACTGTGGGTTGCGGTAGCTGGAAACCCGCTCGCGTAGACGCTC  >  NC_002947/1322717‑1323006

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: