Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,842,896 Δ1 bp 29.7% intergenic (+40/‑257) ccoP‑II → / → PP_4259 cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,842,8960C.29.7% 79.0 / 36.2 37intergenic (+40/‑257)ccoP‑II/PP_4259cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (11/14);  new base . (4/7);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

GGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACC  >  NC_002947/4842749‑4843042
                                                                                                                                                    |                                                                                                                                                   
gggCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                    <  5:70726/149‑1 (MQ=255)
   cAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGa                                                                                                                                                 <  4:98503/149‑1 (MQ=255)
     gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATcc                                                                                                                                              >  2:245514/1‑149 (MQ=255)
         ccGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt‑t                                                                                                                                           <  2:130764/148‑1 (MQ=255)
           ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCtgtg                                                                                                                                          >  1:62615/1‑144 (MQ=255)
            gCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAg                                                                                                                                       <  2:192505/149‑1 (MQ=255)
               cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCATACGTCTACACCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGCGTAAACCCGCTCCTAAAAAGGATCCT‑Tga                                                                                                                                         >  3:59315/1‑143 (MQ=255)
                  gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                 <  6:80650/149‑1 (MQ=255)
                  gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                 <  1:245514/149‑1 (MQ=255)
                  gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCCGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTCGGAGCgg                                                                                                                                 <  7:247605/149‑1 (MQ=38)
                   aacaacTGCATGGCAACGAAAATGTGCACTTGCTGGGGGCTTACGTCTACATCCTTGCCCACCATGCTGAGCCGGCCAAGGCTGAGTTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg                                                                                                                                 <  1:282957/149‑1 (MQ=37)
                   aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAGGAGCggg                                                                                                                                <  2:62615/149‑1 (MQ=255)
                           cAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg                                                                                                                                 >  3:55499/1‑140 (MQ=255)
                           cAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg                                                                                                                                 <  4:55499/140‑1 (MQ=255)
                           cAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTAcc                                                                                                                          <  3:144552/147‑1 (MQ=255)
                                       aaGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                            >  2:42842/1‑149 (MQ=255)
                                                  gCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAg                                                                                                 >  6:144879/1‑149 (MQ=255)
                                                           ttACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                    >  4:328644/1‑90 (MQ=255)
                                                           ttACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAcaca                                                                                                                                                    <  3:328644/90‑1 (MQ=255)
                                                              cGTCTACAGCCTGTCGCACCAGGCTGCGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCa                                                                                      <  5:100775/148‑1 (MQ=255)
                                                               gTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAgg                                                                                    >  6:327775/1‑149 (MQ=255)
                                                                   tcAGCCTGGTGCGACAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCAAAACAGGATGCA‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGGACAGCATTTATCGGAGGCAAg                                                                                >  8:39756/2‑149 (MQ=25)
                                                                   tcAGCCTGGTGCGACAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCACCGACACAGGCTCCG‑TGTATGAGCGGGTTTAACCGCGAAGAGGCCGGTACTGGCACCGCCTTTCTCGCAGGCAAg                                                                                >  8:148676/2‑149 (MQ=37)
                                                                    cAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCaa                                                                                 <  5:144879/147‑1 (MQ=255)
                                                                           tCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCt                                                                          <  5:327775/147‑1 (MQ=255)
                                                                             gCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAACCCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGa                                                                       <  8:171398/148‑1 (MQ=255)
                                                                                  aGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCTGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                 <  5:184531/149‑1 (MQ=255)
                                                                                            ggCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAgg                                                                                                                                      <  4:43624/70‑1 (MQ=17)
                                                                                            ggCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACA‑AGGATCCTGTGTAgg                                                                                                                                      >  3:43624/1‑70 (MQ=17)
                                                                                                        g‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGTAGCGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGttt                                            <  7:39756/148‑1 (MQ=255)
                                                                                                               caggcCTGGCCTCTTGGGGGGAAAACCCGCTCCTACACAGGAGCCT‑TGTAGCAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCa                                     <  7:148676/149‑1 (MQ=38)
                                                                                                                      ggcctCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTGTCGCAGGGAAGTGACCGGGATCAATTGACACCCGCCATCACACGTTTCATACCACGAAcc                               >  6:321337/1‑148 (MQ=37)
                                                                                                                       gcctCTTCGCGGGTTAACCCGCTGCGTCACAGGATGCG‑TGTAGGAGCGGGTTTACCCGCGGAGAGGCCGGTACGGGCACAGGATTTATCGCAGGGAAGTGGCCGGGATGGATTGCAACCCGCCATAACACGTTTCCTACCACGCAccc                              >  6:316568/1‑148 (MQ=19)
                                                                                                                         ctctTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAAc                           <  1:42842/149‑1 (MQ=255)
                                                                                                                                  ggTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATAcc                                      >  3:43239/1‑129 (MQ=255)
                                                                                                                                  ggTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATAcc                                      <  4:43239/129‑1 (MQ=255)
                                                                                                                                         ccGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCg             >  4:99287/1‑147 (MQ=255)
                                                                                                                                         ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGttt                                            >  3:256220/1‑116 (MQ=25)
                                                                                                                                         ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGttt                                            <  4:256220/116‑1 (MQ=25)
                                                                                                                                             tCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGc       >  8:339803/1‑149 (MQ=255)
                                                                                                                                                  cacaGGATCCT‑TGTAGGAGCTGGTTTACCCGCGAAGAGTACGGTACTGGCACAGCATTTCTCGCATGCAAGTGAACTGGATCAATTGACACAAGCCATAACACCTTTAATACCACGAACCCAACGCGACTAAAGGTCTCAGTGCGAAcc  >  4:270635/1‑149 (MQ=255)
                                                                                                                                                  cacaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa    <  7:339803/147‑1 (MQ=255)
                                                                                                                                                   acaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAt                                                                  <  8:183238/84‑1 (MQ=255)
                                                                                                                                                   acaGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAcc  >  6:112902/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGCAGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAG‑TAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACC  >  NC_002947/4842749‑4843042

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