Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 3,951,159 T→C intergenic (‑123/+72) PP_3486 ← / ← PP_3487 cytochrome c/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,951,1590TC97.1% 114.3 / NA 35intergenic (‑123/+72)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (21/13);  minor base G (0/1);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCG  >  NC_002947/3951017‑3951298
                                                                                                                                              |                                                                                                                                           
gcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcgc                                                                                                                                       <  5:246712/149‑1 (MQ=255)
   gAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGAc                                                                                                                                    >  4:85817/1‑149 (MQ=255)
    aGGTCGAAGTTGCCATGCCAAGCACGCAGTGATGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCGCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGTGAGCGCTGCGCGCTCCTTTCGCGAc                                                                                                                                    <  5:269967/148‑1 (MQ=38)
     ggTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGc                                                                                                                                  <  3:85817/149‑1 (MQ=255)
         gAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCATGACTCGCCCAAGGGCAGCGCACCCGCCCCCCCAAGGGCAGCACAAGCGCAGAAAGAGCGGAGCGCTTCGCCCTCCATTAGCGACGCacag                                                                                                                              >  8:100785/1‑146 (MQ=255)
            gTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCg                                                                                                                           <  2:32706/149‑1 (MQ=255)
                                    aGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcgcg                                                                                                                                      >  6:221270/1‑114 (MQ=255)
                                    aGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcgcg                                                                                                                                      <  5:221270/114‑1 (MQ=255)
                                      ccATATCGCAAGCCCTGTGGGAGCGGCCATTCGTCGCGAAAGGGCTGAGCCCCACCCCCGGCAGTGGCATCACAAACCCAGAAACAGGGGAGCCCTTACCGCCACAATAGCGACGCAAGGCCGAACCCACAAAAACCGCTGCCACacag                                                                                                 >  6:193177/1‑149 (MQ=255)
                                                                        cgcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAATGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                                                                  >  6:351544/1‑146 (MQ=255)
                                                                        cgcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAATGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                                                                  <  5:351544/146‑1 (MQ=255)
                                                                         gcgAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCaccgac                                                              >  5:301438/1‑149 (MQ=255)
                                                                            aatgggCGGGGCCGCAGCCCGGGACGTGGAAGCCCAAGCGCATAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCTACGCATGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCgg                                                           <  7:100785/146‑1 (MQ=255)
                                                                              ggggCTGCTGCGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTAGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGc                                                          >  1:186990/2‑148 (MQ=255)
                                                                                gggggcgcTGCAGCCCCGGCAGTGGCAGCTCAAGCGCACAAAGAGGGGAGCGGGGCGCGCTCCTTTCGCTACGAAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg                                                       <  6:70183/145‑1 (MQ=255)
                                                                                    gcgcAGCAGCCCCGGCTGTGGCATCTCACCCGCATAATCAGGGGAGCGCTGCGCGCGCGTTCCGTCACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGc                                                   <  5:193177/149‑1 (MQ=255)
                                                                                          cagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTcc                                                                                                                       >  5:268299/1‑75 (MQ=255)
                                                                                          cagcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTcc                                                                                                                       <  6:268299/75‑1 (MQ=255)
                                                                                            gcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                           >  1:182205/1‑149 (MQ=255)
                                                                                            gcCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                           >  5:90038/1‑149 (MQ=255)
                                                                                                  gCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAAc                                          <  1:354284/144‑1 (MQ=255)
                                                                                                  gCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAAc                                          >  2:354284/1‑144 (MQ=255)
                                                                                                          aGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCt                             >  4:303600/1‑149 (MQ=255)
                                                                                                            cTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCCTGCGCTg                            >  8:8897/1‑148 (MQ=255)
                                                                                                               aaGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                              >  3:138988/1‑127 (MQ=255)
                                                                                                               aaGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                              <  4:138988/127‑1 (MQ=255)
                                                                                                               aaGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGc                         >  4:53993/1‑148 (MQ=255)
                                                                                                                aGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGTGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCaa                       >  2:12805/1‑149 (MQ=255)
                                                                                                                         agagGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTg              >  8:219723/1‑149 (MQ=255)
                                                                                                                          gaaaggAGCGCGCAGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTgg             >  3:144751/5‑149 (MQ=255)
                                                                                                                             gggAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCTATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCgg          >  1:299770/1‑149 (MQ=255)
                                                                                                                               gAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCCGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTg        >  7:73486/1‑149 (MQ=255)
                                                                                                                                  cgcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCAt     >  3:323433/1‑149 (MQ=255)
                                                                                                                                   gcTGCGCGGTCCTTGCGCTAGGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATc    <  8:77774/149‑1 (MQ=255)
                                                                                                                                       cgcgcTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCg  <  3:303600/147‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
GCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAGCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCCG  >  NC_002947/3951017‑3951298

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: