Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,148,320 C→T 8.7% A74T (GCC→ACC)  PP_5695 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,148,3200CT 8.7% 128.1 / 3.0 46A74T (GCC→ACC) PP_5695hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base T (2/2);  total (25/21)
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.63e-01

GAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGGCGCTCTGGCAAGGGATGCCGGCCT  >  NC_002947/5148177‑5148454
                                                                                                                                               |                                                                                                                                      
gaagaaATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCaa                                                                                                                                     <  6:340654/147‑1 (MQ=255)
 aagaaATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGTCAAgg                                                                                                                                   >  5:346428/1‑148 (MQ=255)
                 cTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCg                                                                                                                    >  3:63731/1‑147 (MQ=255)
                         gCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACagag                                                                                                                            >  4:69295/1‑131 (MQ=255)
                         gCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACagag                                                                                                                            <  3:69295/131‑1 (MQ=255)
                                   ccTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGt                                                                                                <  4:63731/149‑1 (MQ=255)
                                           ccccTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACag                                                                                                                              >  8:153239/1‑111 (MQ=255)
                                           ccccTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACag                                                                                                                              <  7:153239/111‑1 (MQ=255)
                                               ttACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCAt                                                                                     >  6:5464/1‑148 (MQ=255)
                                                  cGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAAt                                                                                 <  5:5464/149‑1 (MQ=255)
                                                   gAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTACATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATc                                                                                >  6:72763/1‑149 (MQ=255)
                                                      ccAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCt                                                                                                     <  4:22470/125‑1 (MQ=255)
                                                      ccAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCt                                                                                                     >  3:22470/1‑125 (MQ=255)
                                                      ccAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGt                                                                             >  1:97735/1‑149 (MQ=255)
                                                        aaGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTc                                                                            >  4:88240/1‑148 (MQ=255)
                                                        aaGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTc                                                                            >  7:604/1‑148 (MQ=255)
                                                           gCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAg                                                                         <  1:198621/148‑1 (MQ=255)
                                                           gCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAg                                                                         >  2:198621/1‑148 (MQ=255)
                                                             cGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGTCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAggg                                                                       >  8:130424/1‑148 (MQ=255)
                                                              gCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGcc                                                                                                                                       <  6:237309/83‑1 (MQ=255)
                                                                    cAAAGGAGCGGCGGGATGATGGTGTTTTAGTGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCa                                                               <  6:22954/149‑1 (MQ=255)
                                                                            cggcggGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCg                                                         >  5:3856/1‑147 (MQ=255)
                                                                                  gatgatGGTGTTTAAGGGGTTTACATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCAt                                                 <  5:72763/149‑1 (MQ=255)
                                                                                      atgGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGTCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCaa                                              <  7:130424/148‑1 (MQ=255)
                                                                                      atgGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGTCAAGGCACAGAGGTAGGGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCaa                                              <  6:346428/148‑1 (MQ=255)
                                                                                      atgGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCaa                                              <  3:88240/148‑1 (MQ=255)
                                                                                           gtTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCg                                         >  6:225411/1‑148 (MQ=255)
                                                                                                aGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGt                                                            <  7:136256/124‑1 (MQ=255)
                                                                                                aGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGt                                                            >  8:136256/1‑124 (MQ=255)
                                                                                                 ggggTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTg                                  >  3:339915/1‑149 (MQ=255)
                                                                                                   ggTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGtt                                >  4:297862/1‑149 (MQ=255)
                                                                                                         catcatGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCgg                           <  5:225411/148‑1 (MQ=255)
                                                                                                          atcatGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGgc                         <  8:604/149‑1 (MQ=255)
                                                                                                           tcatGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCa                                                                                      <  3:228225/87‑1 (MQ=255)
                                                                                                           tcatGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCa                                                                                      >  4:228225/1‑87 (MQ=255)
                                                                                                           tcatGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGgcg                        >  4:173599/1‑149 (MQ=255)
                                                                                                              tGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAAt                                                                                 <  3:326008/89‑1 (MQ=255)
                                                                                                              tGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAAt                                                                                 >  4:326008/1‑89 (MQ=255)
                                                                                                              tGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGGCGctc                     >  7:86274/1‑149 (MQ=255)
                                                                                                                cccAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGGCGCTCTg                   >  5:50692/1‑149 (MQ=255)
                                                                                                                   aaTTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCAGTCCAGGGTCATGTCACGTGCGTAGCCCATCCACTCCGAAGCTTGTTCGCGGGCGCTCTGGc                 >  1:324685/1‑148 (MQ=255)
                                                                                                                        aGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGGCGCTCTGGCAAggg            <  6:50692/148‑1 (MQ=255)
                                                                                                                        aGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGGCGCTCTGGCAAggg            <  3:297862/148‑1 (MQ=255)
                                                                                                                           atcatcCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGGCGCTCTGGCAAGGGATGc        <  8:86274/149‑1 (MQ=255)
                                                                                                                                ccatcAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGACCAGATTAGCTGAGGTTTGTCTTACATCAATACGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGGCGCTCTGGCAAGGGATGCCggac   >  7:69723/1‑147 (MQ=255)
                                                                                                                                 catcAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGGCGCTCTGGCAAGGGATGCCGGCCt  <  6:3856/149‑1 (MQ=255)
                                                                                                                                      gcgcCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCAt                                                 >  4:325590/1‑97 (MQ=255)
                                                                                                                                               |                                                                                                                                      
GAAGAAATGAGCGAGGCCTGAGATTGCCGGGGCTGCCTTGCAGCCCCTTACGAGCCAAGGCCGCTCTGCAAAGGAGCGGCGGGATGATGGTGTTTAAGGGGTTTTCATCATGCCCAATTCAGCATCATCCATCAGCGCCTTGGCCAAGGCACAGAGGTAGTGCGAAGCCCAGATTAGCTGAGGTTTGTCTTCCATCAATCCGTCCAGGGTCATGTCACGTGCGTAGCCCATCAACTCCGAAGCTTGTTCGCGGGCGCTCTGGCAAGGGATGCCGGCCT  >  NC_002947/5148177‑5148454

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