Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,631,597 C→T 10.2% intergenic (+1371/‑64) PP_4094 → / → PP_4096 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,631,5970CT10.2% 102.3 / 3.1 39intergenic (+1371/‑64)PP_4094/PP_4096hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (16/19);  new base T (2/2);  total (18/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.49e-01

CAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAGGAGCAACAACATGCGCAAACGGAT  >  NC_002947/4631455‑4631731
                                                                                                                                              |                                                                                                                                      
cAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATa                                                                                                                                  <  7:129587/149‑1 (MQ=255)
 aaaCCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCa                                                                                                                                    >  7:18878/1‑146 (MQ=255)
 aaaCCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATa                                                                                                                                  <  8:18878/148‑1 (MQ=255)
 aaaCCGCCGCGGCGCCACTTCTATCGAAAGCCTTCAACCCCACAGCGTCCATCGTTGAAACAAACGACTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTCGCACTGGCAAGGACGGAACCTGAACAGCAAGACAGCGTCATAg                                                                                                                                 >  1:290920/1‑149 (MQ=255)
        cgcgGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATc                                                                                                                          <  6:241183/149‑1 (MQ=255)
                           aaaGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGt                                                                                                                                      >  1:151095/1‑118 (MQ=255)
                           aaaGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGt                                                                                                                                      <  2:151095/118‑1 (MQ=255)
                              gCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGAc                                                                                                    <  3:223331/149‑1 (MQ=255)
                                    aaCCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCGTAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATg                                                                                              >  7:290537/1‑149 (MQ=255)
                                      ccccACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATgcc                                                                                            >  1:279470/1‑149 (MQ=255)
                                       cccACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATgccg                                                                                           <  6:176475/149‑1 (MQ=255)
                                                      gTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACt                                                                             >  4:119522/1‑148 (MQ=255)
                                                      gTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTc                                                                            <  3:119522/149‑1 (MQ=255)
                                                       ttGAAACTAAGTAGTAGCCGGTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCa                                                                           <  2:152664/149‑1 (MQ=255)
                                                            aCAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGTGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATg                                                                      >  3:308145/1‑149 (MQ=255)
                                                              aaaGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTg                                                                                                       <  5:171419/114‑1 (MQ=255)
                                                              aaaGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTg                                                                                                       >  6:171419/1‑114 (MQ=255)
                                                                       gCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCa                                                           >  3:178956/1‑149 (MQ=255)
                                                                          gttgttGCGCAGCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCGAGCAAGTCGGCATGAGTTGACTCGATGCCGCCCTCCGCATGACTCTGGATGGAAACGGCGCACCg                                                        <  2:290920/149‑1 (MQ=255)
                                                                             gttCCGCACCGTGCAAACCGTTAATATCTGGCACTGGCAAGGACGGAACCTTAAGATCAAGCCAGTGTCATAGCATAATCCAGCAACTAGGAATCAGTTGACTCTATGCCGACCTCCGAATGACTCAGGATGGAAACGGCTCACCGGt                                                      <  4:308145/148‑1 (MQ=255)
                                                                                    aCCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGcc                                              <  4:178956/149‑1 (MQ=255)
                                                                                             cccGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCt                                     <  8:290537/149‑1 (MQ=255)
                                                                                                    atatCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTg                                                                                                       >  5:133220/1‑76 (MQ=255)
                                                                                                    atatCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTg                                                                                                       <  6:133220/76‑1 (MQ=255)
                                                                                                       tCTGGCACTGGCAAGGAAGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGc                             >  8:151991/1‑147 (MQ=255)
                                                                                                         tggcactggcaAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCa                           <  2:55395/147‑1 (MQ=255)
                                                                                                         tggcactggcaAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCa                           >  1:55395/1‑147 (MQ=255)
                                                                                                         tggcactggcaAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAg                          >  4:167531/1‑148 (MQ=255)
                                                                                                          ggcactggcaAGGACGGAACCTGAAGAGCAAGCCAGTGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAg                          >  8:258187/1‑147 (MQ=255)
                                                                                                          ggcactggcaAGGACGGAACCTGAAGAGCAAGCCAGTGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAg                          <  7:258187/147‑1 (MQ=255)
                                                                                                           gcactggcaAGGAAGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAGGAg                       <  2:279470/149‑1 (MQ=255)
                                                                                                               tggcaAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGt                                            >  4:129332/1‑124 (MQ=255)
                                                                                                               tggcaAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGt                                            <  3:129332/124‑1 (MQ=255)
                                                                                                                     ggACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAgg                         <  2:230989/137‑1 (MQ=255)
                                                                                                                     ggACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAgg                         >  1:230989/1‑137 (MQ=255)
                                                                                                                      gACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAGGAGCAACAACATgc            <  3:167531/149‑1 (MQ=255)
                                                                                                                         ggAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAGGAGCAACAACATgcgc          >  8:224459/1‑148 (MQ=255)
                                                                                                                               tGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAGGAGCAACAACATGCGCAAACGGa   >  1:276258/1‑149 (MQ=255)
                                                                                                                                 aaGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAGGAGCAACAACATGCGCAAACGGAt  <  7:151991/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                      
CAAACCGCCGCGGCGCCACTTCTATCGAAAGCGTTCAACCCCACAGCGTCCATCGTTGAAACAAAGGAGTAGCCGTTGTTCCGCACCGTGCCACCCGTTCATATCTGGCACTGGCAAGGACGGAACCTGAAGAGCAAGCCAGCGTCATAGCCTCATCCAGCAACTCGGCATCAGTGGACTCGATGCCGCCCTCCGCATGACTCAGGATGGAAACGGCGCACCGGTTTCGAGCCGTTCGATCTGGAATTGCCAGGAGCAACAACATGCGCAAACGGAT  >  NC_002947/4631455‑4631731

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