Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,842,896 Δ1 bp 17.6% 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.17.6% 90.8 / 15.9 34intergenic (+40/‑257)ccoP‑II/PP_4259cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (12/16);  new base . (2/4);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.03e-01

AGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAC  >  NC_002947/4842753‑4843041
                                                                                                                                               |                                                                                                                                                       
aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGAt                                                                                                                                                    >  7:97963/1‑149 (MQ=255)
 gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATcc                                                                                                                                                  >  7:219881/1‑149 (MQ=255)
  aTGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATcc                                                                                                                                                  <  8:198432/149‑1 (MQ=255)
  aTGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCt                                                                                                                                                 >  2:284252/1‑149 (MQ=255)
   tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATcc                                                                                                                                                  <  5:135445/148‑1 (MQ=255)
      cGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGACGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGt                                                                                                                                             <  8:216719/149‑1 (MQ=255)
        gCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                                                                                           <  4:49400/149‑1 (MQ=255)
         cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAgg                                                                                                                                          <  8:219881/149‑1 (MQ=255)
          ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                                         <  8:252588/149‑1 (MQ=255)
            aGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCACGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                                       <  8:358116/149‑1 (MQ=255)
             ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTTGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACATGATCCTTGTAGGAGCg                                                                                                                                      <  4:291986/149‑1 (MQ=255)
             ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCg                                                                                                                                      <  1:284252/149‑1 (MQ=255)
             ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCg                                                                                                                                      <  8:44432/149‑1 (MQ=255)
              gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                     >  2:43021/1‑149 (MQ=255)
               aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                                    <  5:133405/149‑1 (MQ=255)
                  aacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGtt                                                                                                                                  >  1:321169/1‑148 (MQ=255)
                    cTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTg                                                                                                                                              <  2:361690/134‑1 (MQ=255)
                    cTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTg                                                                                                                                              >  1:361690/1‑134 (MQ=255)
                     tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                               >  2:52441/1‑132 (MQ=255)
                     tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                               <  1:52441/132‑1 (MQ=255)
                                gACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGc                                                                                                              <  4:260916/149‑1 (MQ=255)
                                   acggTGCACTTGATGGCGGCTTACGTCTACAGCCTGCCGCACCAGGCTGAGCCGGCCAAGGCTGTGTACGGCATGCGTGGCCTCTTCGCGGGTAACCCCGCTCCTACACTGCATCCTTGTTGGAGGGGTTTAACCCGC‑‑‑‑‑GAAGAGGCCgg                                                                                                           <  2:321169/147‑1 (MQ=255)
                                       tGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACt                                                                                                       >  2:115022/1‑149 (MQ=255)
                                              gCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGcaca                                                                                                 <  1:43021/148‑1 (MQ=255)
                                                      cTTACTTCTACCGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGTTTACGGCAGGCCGGGCCTCTTCGCGGTTAACCCCGCTCCTACACAGCATCCTTGTCGGCGCGGTTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTtctc                                                                                        <  2:126887/149‑1 (MQ=37)
                                                               aCAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAg                                                                               >  7:100891/1‑149 (MQ=255)
                                                                   ccTGTCGCACCAGGCTGAGCCGGCCAAGGCTTAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGAc                                                                           <  7:281936/149‑1 (MQ=255)
                                                                   ccTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGAc                                                                           <  8:97963/149‑1 (MQ=255)
                                                                                                                      ctTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCACTTGACACCCGCCATAACACGGTTCATACCACGAACCCACcgc                        >  2:329316/1‑149 (MQ=255)
                                                                                                                          gcgGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGAc                     >  3:143987/1‑148 (MQ=255)
                                                                                                                              gTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAg                >  5:390753/1‑149 (MQ=255)
                                                                                                                                  aCCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCCCCCATAGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTaaa                 <  1:329316/149‑1 (MQ=255)
                                                                                                                                    ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCg            >  4:59657/1‑147 (MQ=255)
                                                                                                                                           tACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa   >  8:2675/1‑149 (MQ=255)
                                                                                                                                            acacAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAc  >  1:176478/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
AGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGC‑‑‑‑‑GAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAC  >  NC_002947/4842753‑4843041

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: