Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,842,949:1 +G 41.9% 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.G41.9% 31.8 / 55.0 31intergenic (+48/‑249)ccoP‑II/PP_4259cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (9/9);  new base G (7/6);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAA  >  NC_002947/4842811‑4843085
                                                                                                                                           |                                                                                                                                        
ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                               <  2:331176/150‑1 (MQ=255)
ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                               <  2:601946/150‑1 (MQ=255)
   acaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAgg                                                                                                                                    >  2:508146/1‑142 (MQ=255)
   acaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAgg                                                                                                                                    <  1:508146/142‑1 (MQ=255)
         gCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                                   <  2:453453/137‑1 (MQ=255)
         gCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                                   >  1:453453/1‑137 (MQ=255)
                   gACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                     <  2:28209/125‑1 (MQ=255)
                   gACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAg                                                                                                                                     >  1:28209/1‑125 (MQ=255)
                                                   cAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                              <  2:611298/100‑1 (MQ=21)
                                                   cAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                              >  1:611298/1‑100 (MQ=21)
                                                      ccTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGAcc                                                                         <  2:95767/150‑1 (MQ=255)
                                                      ccTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGAcc                                                                         >  1:29638/1‑150 (MQ=255)
                                                            gCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                                                                                                   <  2:22377/86‑1 (MQ=21)
                                                            gCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                                                                                                   >  1:22377/1‑86 (MQ=21)
                                                                              cAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATaa                                                 >  2:284230/1‑150 (MQ=255)
                                                                                       gTAAGGCAGGCCTGGGCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCa                                        <  2:429092/150‑1 (MQ=255)
                                                                                            gcaggcCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCa                                   >  1:338637/1‑150 (MQ=255)
                                                                                                gcctggcctCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGaa                               <  1:284230/150‑1 (MQ=255)
                                                                                                   cggcctCTTCGCGGGTACACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                               <  1:449508/114‑1 (MQ=18)
                                                                                                   cggcctCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                               >  2:449508/2‑115 (MQ=18)
                                                                                                          ttCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGa                     >  2:603981/1‑150 (MQ=255)
                                                                                                                  ctctcccGCCCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGCCGGCCCCGGGACCACCACCTCATCACTCGAAGACAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGt             <  1:43931/146‑1 (MQ=255)
                                                                                                                      cccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa         >  2:383599/1‑150 (MQ=255)
                                                                                                                       ccGCTCCTACACAGGATCCTTGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATGTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATAGCACGAACGCAACGCGACTAAAgggcagct        >  2:43931/1‑144 (MQ=39)
                                                                                                                       ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGt                                                                             >  1:183221/1‑81 (MQ=18)
                                                                                                                       ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGt                                                                             <  2:183221/81‑1 (MQ=18)
                                                                                                                              ctacaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa  >  2:476988/6‑150 (MQ=255)
                                                                                                                                    ggATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                             >  2:222317/1‑36 (MQ=40)
                                                                                                                                    ggATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                             <  1:222317/36‑1 (MQ=255)
                                                                                                                                      aTCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa  <  2:172521/141‑1 (MQ=255)
                                                                                                                                      aTCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa  >  1:172521/1‑141 (MQ=255)
                                                                                                                                          ttGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                               <  2:372436/76‑1 (MQ=255)
                                                                                                                                          ttGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                               >  1:372436/1‑76 (MQ=255)
                                                                                                                                           |                                                                                                                                        
GGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAA  >  NC_002947/4842811‑4843085

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