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

AGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCG  >  NC_002947/4842753‑4843038
                                                                                                                                                |                                                                                                                                               
aGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATc                                                                                                                                           >  2:83840/1‑149 (MQ=255)
        gCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                   >  8:221990/1‑149 (MQ=255)
         cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAgg                                                                                                                                  <  7:262711/149‑1 (MQ=255)
           cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAg                                                                                                                                <  7:221990/149‑1 (MQ=255)
            aGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                               <  8:180430/149‑1 (MQ=255)
            aGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                               <  8:270113/149‑1 (MQ=255)
            aGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGc                                                                                                                               <  1:83840/149‑1 (MQ=255)
               aacaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                            <  7:77100/149‑1 (MQ=255)
                  aacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGtt                                                                                                                          >  5:111387/1‑148 (MQ=255)
                             aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                >  3:98332/1‑147 (MQ=255)
                                  cAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                                                                                   >  7:270438/1‑123 (MQ=255)
                                  cAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                                                                                   <  8:270438/123‑1 (MQ=255)
                                   aaGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGTAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                        <  2:219546/149‑1 (MQ=255)
                                     ggTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                       >  5:10152/1‑116 (MQ=255)
                                     ggTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                       <  6:10152/116‑1 (MQ=255)
                                        gCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTg                                                                                                   >  6:146914/1‑149 (MQ=255)
                                                   cggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                            <  6:281668/113‑1 (MQ=255)
                                                   cggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                            >  5:281668/1‑113 (MQ=255)
                                                         aCGTCTACAGCCTGTCGCACCAGGCTTAGCCGGCCATGGCTTAGTATGGCTGGCGTGGCCTCTTCGCGGGTAACCCCGCTCCTACACAGCATCCTTGTAGTTGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGc                                                                                   <  6:111387/148‑1 (MQ=37)
                                                         aCGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCCCTCCAACAAAGGATACATGTAGGACCGGGTTTACCCGCCAAGACGCCCGTAATGGCACAGCAtatatcc                                                                                    >  2:92450/1‑141 (MQ=37)
                                                           gTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCGGCGTACCAGGGTACGTGTAGGAGCGGGGTTTCCCGCGAAGAGGGCGGGTCTGGCACCGCATTtctggacgg                                                                                >  7:127871/1‑143 (MQ=25)
                                                              tACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCGCCTACCCAGGATCCATGTAGGGGCGGGTTTACCCGCGCAGCCGCCGGTACGGGCACCGCATTTGTCGCAggg                                                                               >  3:21287/1‑146 (MQ=37)
                                                                   ccTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTTTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGAc                                                                        <  5:279516/149‑1 (MQ=255)
                                                                   ccTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAACCGCGCCTACACAGGATCCTTGTAAGAGCCGGTTTACACGCGAAGAGGCCGGTAAAGGCACAGCATATCTCGCAGTCAAATgc                                                                         >  3:51438/1‑147 (MQ=37)
                                                                       tCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCt                                                                      <  4:98332/147‑1 (MQ=255)
                                                                        cGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                            <  4:3994/92‑1 (MQ=21)
                                                                        cGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                            >  3:3994/1‑92 (MQ=21)
                                                                         ccACCAGGCTGAGCCGGCCAAGGCTGAGTACGGCCGGCCTGGCCTCTTCGCGGGTAAACCCCCTCCTACACCGTATCCTTGTAGCAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGAt                                                                  <  6:98696/148‑1 (MQ=37)
                                                                           aCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATc                                                                 >  7:212867/1‑148 (MQ=255)
                                                                               ggatctgCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                             <  7:241180/142‑1 (MQ=255)
                                                                                      ccGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTAAACCCCCCCGACACACGGTGCATGGAAGGCCGGGGTTAACACCGAAGAGGCGGGGGCAGGGACAGGATTTGTGGCCGGCAAGGGACCAGGGTGAAtttttaaag                                                      >  8:241180/1‑141 (MQ=255)
                                                                                                 tGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCCCAGCATTTCTCGCAGGGaa                                                                             >  7:148563/1‑114 (MQ=255)
                                                                                                        ggcaggcCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATAc                                   <  8:212867/149‑1 (MQ=255)
                                                                                                                     tctTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                                 >  4:193777/1‑42 (MQ=255)
                                                                                                                         ggcgGGTAAACCCGCTCATACTCAGTATCCTTGTTGTAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAAcgc                     <  1:92450/145‑1 (MQ=39)
                                                                                                                         cgcgGTAAAACCCGTCCATACACAGTATCCTTGTTGGAGCGGTTTAACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAAcgc                     <  7:257731/146‑1 (MQ=255)
                                                                                                                             tgTAAACCCTCTCATACACATCATCCTTTTAGTAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACATGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTaaa              <  3:233436/148‑1 (MQ=39)
                                                                                                                              gTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                                         >  2:19625/1‑57 (MQ=255)
                                                                                                                                    ccGCTCCTACACCGTATCCTTGTCGGTGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa       <  8:39641/149‑1 (MQ=255)
                                                                                                                                    ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa       <  5:146914/149‑1 (MQ=255)
                                                                                                                                      gCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGt     >  5:83869/1‑149 (MQ=255)
                                                                                                                                           cTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCg  >  2:70159/1‑148 (MQ=255)
                                                                                                                                            tACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAt                                                              >  4:141861/1‑87 (MQ=255)
                                                                                                                                            tACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAt                                                              <  3:141861/87‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AGATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCG  >  NC_002947/4842753‑4843038

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: