Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,842,949:1 +G 18.2% 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,9491.G18.2% 88.6 / 21.1 33intergenic (+48/‑249)ccoP‑II/PP_4259cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (12/15);  new base G (4/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

CCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAT  >  NC_002947/4842803‑4843098
                                                                                                                                                   |                                                                                                                                                     
ccGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgtgt                                                                                                                                                    <  1:508331/150‑5 (MQ=255)
     ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAg                                                                                                                                               <  1:126636/150‑1 (MQ=38)
      cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCCAGGTTCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                                              <  1:8457/150‑1 (MQ=39)
        ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                            <  1:631674/150‑1 (MQ=255)
        ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                            <  1:20644/150‑1 (MQ=255)
                              aaGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt                                                                                                                                                        <  1:426553/117‑1 (MQ=255)
                              aaGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCt                                                                                                                                                        >  2:426553/1‑117 (MQ=255)
                                           tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAt                                                                                                         <  2:775250/150‑1 (MQ=255)
                                                            aGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGc                                                                                                                <  1:223394/126‑1 (MQ=255)
                                                            aGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGc                                                                                                                >  2:223394/1‑126 (MQ=255)
                                                                 gTCGCACCAGGCTGAGCCGGTCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                        >  2:818353/1‑113 (MQ=255)
                                                                 gTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                        <  1:818353/113‑1 (MQ=255)
                                                                      aCCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCATGGCCTCTTCGAGGGTAAACCCGCTCCTACAGACGAGCCTTGTATGAGCGGGTTTAACCGCGAAGAAGCCGGTAGAGGAACAGCATTTATCGCAGGCAAGTTACCTGGATCaa                                                                              >  1:311458/1‑150 (MQ=255)
                                                                        cAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            <  1:14167/150‑1 (MQ=255)
                                                                          ggCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGa                                                                          <  1:407849/150‑1 (MQ=255)
                                                                               aGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAATGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACAccc                                                                     >  1:817713/1‑150 (MQ=255)
                                                                                gCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            <  1:445988/142‑1 (MQ=255)
                                                                                gCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            >  2:445988/1‑142 (MQ=255)
                                                                                              aGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTc                                                      <  1:13412/150‑1 (MQ=255)
                                                                                                        gcctggcctCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTAATGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGaa                                            <  2:817713/150‑1 (MQ=255)
                                                                                                                         gTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAgg                           >  2:624273/1‑150 (MQ=255)
                                                                                                                              cccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCTGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa                      <  2:311458/150‑1 (MQ=255)
                                                                                                                              cccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGa                                  <  1:763852/138‑1 (MQ=255)
                                                                                                                              cccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGa                                  >  2:763852/1‑138 (MQ=255)
                                                                                                                              cccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa                      >  2:789627/1‑150 (MQ=255)
                                                                                                                               ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                     >  2:826777/1‑150 (MQ=255)
                                                                                                                               ccGCTCCTACAAGGTTCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCGGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                     >  2:806480/1‑150 (MQ=38)
                                                                                                                                 gCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            >  2:644578/1‑93 (MQ=255)
                                                                                                                                 gCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            <  1:644578/93‑1 (MQ=255)
                                                                                                                                 gCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg                   >  1:515152/1‑150 (MQ=255)
                                                                                                                                  cTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGAAGTGc                  >  2:28019/1‑150 (MQ=255)
                                                                                                                                   tCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGaa                                            <  2:139884/123‑1 (MQ=255)
                                                                                                                                   tCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGaa                                            >  1:139884/1‑123 (MQ=255)
                                                                                                                                     cTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa               <  1:826777/150‑1 (MQ=255)
                                                                                                                                      ctacaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa               >  1:589146/6‑150 (MQ=255)
                                                                                                                                           aGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATa         >  2:666778/1‑150 (MQ=255)
                                                                                                                                                  ttGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAt  >  1:229243/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAT  >  NC_002947/4842803‑4843098

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