Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,842,904:1 +G 17.4% intergenic (+48/‑249) 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,9041.G17.4% 128.9 / 30.8 46intergenic (+48/‑249)ccoP‑II/PP_4259cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (18/20);  new base G (4/4);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCC  >  NC_002947/4842761‑4843051
                                                                                                                                                 |                                                                                                                                                   
gCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                                >  8:221990/1‑149 (MQ=255)
 cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAgg                                                                                                                                               <  7:262711/149‑1 (MQ=255)
   cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAg                                                                                                                                             <  7:221990/149‑1 (MQ=255)
    aGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                                            <  8:270113/149‑1 (MQ=255)
    aGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                                            <  8:180430/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                                                                                                                                                <  8:270438/123‑1 (MQ=255)
                          cAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                                                                                                >  7:270438/1‑123 (MQ=255)
                           aaGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGTAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                     <  2:219546/149‑1 (MQ=255)
                             ggTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                                    <  6:10152/116‑1 (MQ=255)
                             ggTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                                    >  5:10152/1‑116 (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                                                                                                                                         >  3:3994/1‑92 (MQ=21)
                                                                cGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                                         <  4:3994/92‑1 (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                                                                           <  3:141861/87‑1 (MQ=255)
                                                                                                                                    tACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAt                                                                           >  4:141861/1‑87 (MQ=255)
                                                                                                                                      acaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAc            >  7:106231/4‑149 (MQ=255)
                                                                                                                                        aaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGa                                >  7:296256/2‑127 (MQ=255)
                                                                                                                                        aaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGa                                <  8:296256/126‑1 (MQ=255)
                                                                                                                                           gATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATAc      >  5:62283/1‑149 (MQ=255)
                                                                                                                                               cTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATATCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGcc  <  1:70159/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCC  >  NC_002947/4842761‑4843051

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: