Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,301,577 G→C 9.8% intergenic (+67/‑1233) PP_3774 → / → PP_5627 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,301,5770GC 9.8% 108.7 / 3.6 41intergenic (+67/‑1233)PP_3774/PP_5627hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (18/19);  new base C (2/2);  total (20/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 = 4.24e-01

CAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGCCT  >  NC_002947/4301430‑4301723
                                                                                                                                                   |                                                                                                                                                  
cAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgg                                                                                                                                                   <  6:79011/149‑1 (MQ=255)
           cacCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCCCCGCGCTGGGCttt                                                                                                                                        <  2:143303/149‑2 (MQ=255)
            acCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCtt                                                                                                                                         >  1:9071/1‑147 (MQ=255)
            acCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCtt                                                                                                                                         <  2:9071/147‑1 (MQ=255)
            acCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTc                                                                                                                                        >  5:222623/1‑148 (MQ=255)
                aCGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTAt                                                                                                                                     >  5:64274/1‑147 (MQ=255)
                    aaCCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTg                                                                                                                               >  5:54905/1‑149 (MQ=255)
                     aCCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTgg                                                                                                                              <  4:173872/149‑1 (MQ=255)
                       ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCt                                                                                                                                <  4:51187/145‑1 (MQ=255)
                       ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCt                                                                                                                                >  3:51187/1‑145 (MQ=255)
                       ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGc                                                                                                                            <  6:64274/149‑1 (MQ=255)
                       ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGc                                                                                                                            >  3:69720/1‑149 (MQ=255)
                           gccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGaa                                                                                                                        <  1:250641/149‑1 (MQ=255)
                             cgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGATGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACa                                                                                                                      <  7:171243/149‑1 (MQ=255)
                                      cTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGATGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCg                                                                                                             >  7:215155/1‑149 (MQ=255)
                                         gATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGc                                                                                                            <  7:48068/147‑1 (MQ=255)
                                                atcatcAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGcgcg                                                                                                                                                    <  2:250974/100‑1 (MQ=255)
                                                atcatcAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCACcgcg                                                                                                                                                    >  1:250974/1‑100 (MQ=255)
                                                             tctcGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTg                                                                                       >  1:4200/1‑148 (MQ=255)
                                                                 gCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGaa                                                                                    <  6:54905/147‑1 (MQ=255)
                                                                 gCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAg                                                                                   <  1:73266/148‑1 (MQ=255)
                                                                     atgcttggCTGAACGTCGATGCCCAGCCCTGCTATCAATACCCCTCAGTGCGTCGCATGAATGCAAAAAGCCCAGCGCGTGGCTGGGCTTGTATCAGCTGGGCTGAACAGTAGTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACt                                                                              <  7:93091/142‑1 (MQ=255)
                                                                             cTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTggg                                                                      >  8:99125/1‑149 (MQ=255)
                                                                                  cGTCGATGCCCATCCCTGCTATCAACACCACTCAGGGCGTCTCATGAATGCAAAAATGCCACCGCGGTGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTAGACACTGTACTGGAAGTTACTCCATTGGGCGCtt                                                                 <  4:16959/149‑1 (MQ=255)
                                                                                  cGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCtt                                                                 <  7:99125/149‑1 (MQ=255)
                                                                                    tCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTtgt                                                               <  2:4200/149‑1 (MQ=255)
                                                                                           ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCtt                                                                 <  4:5057/140‑1 (MQ=255)
                                                                                           ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCtt                                                                 >  3:5057/1‑140 (MQ=255)
                                                                                           ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAAc                                                         >  7:93827/1‑148 (MQ=255)
                                                                                           ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCAGAACAGTACTTCCGCCAGCCTATACACACTGTACTGGAAGTTACTACAGAGGGCGCTTGTGATAAc                                                         >  7:290618/1‑148 (MQ=255)
                                                                                           ccAGCCCTGCTATCAAAACCCCTCAGAGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAAc                                                         >  6:16178/1‑148 (MQ=255)
                                                                                                      aTCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCgg                                             >  4:117365/1‑149 (MQ=255)
                                                                                                       tCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCgg                                             >  2:175380/1‑148 (MQ=255)
                                                                                                            atcccTCATGGCGTAGGATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTTGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTgg                                       <  8:290618/147‑1 (MQ=255)
                                                                                                            aCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTCATGGAAGTTACTCCATTGGGCGCTTGTGATAACTATAACGCTCGGCATTgg                                       >  7:228574/1‑149 (MQ=255)
                                                                                                            aCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTgg                                       <  4:127085/149‑1 (MQ=255)
                                                                                                                        gTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATg                           <  6:222623/149‑1 (MQ=255)
                                                                                                                        gTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATg                           >  1:118610/1‑149 (MQ=255)
                                                                                                                           gCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGt                        >  2:75981/1‑149 (MQ=255)
                                                                                                                             aTGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTa                                                                                                   <  6:57473/72‑1 (MQ=255)
                                                                                                                                aaTGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCcgg                    <  5:16178/148‑1 (MQ=255)
                                                                                                                                  tgatagaaGCCCAGCCCCGCGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCcggcgg                 >  1:124428/7‑149 (MQ=255)
                                                                                                                                    atagaaGCCCAGCCCCGCGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGc               >  1:246945/5‑149 (MQ=255)
                                                                                                                                     tagaaGCCCAGCCCCGCGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGcc              <  2:246945/146‑1 (MQ=255)
                                                                                                                                                 gcgGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTctgcct  <  4:69720/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAGTACTTCCGCCAGCCTTTACACACTGTACTGGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGCCT  >  NC_002947/4301430‑4301723

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