Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,282,774 G→C 13.9% intergenic (‑166/+63) PP_2010 ← / ← PP_2011 cytochrome b561/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,282,7740GC13.9% 84.7 / 5.3 36intergenic (‑166/+63)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/15);  new base C (2/3);  total (18/18)
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.71e-01

GTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATC  >  NC_002947/2282628‑2282914
                                                                                                                                                  |                                                                                                                                            
gTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTgg                                                                                                                                             >  3:273603/1‑148 (MQ=255)
 tCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTgg                                                                                                                                             >  3:134428/1‑147 (MQ=255)
 tCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTCCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGaa                                                                                                                                           <  2:9378/149‑1 (MQ=255)
       tGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgt                                                                                                                                     <  6:32408/149‑1 (MQ=255)
         gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg                                                                                                                                    <  4:273603/148‑1 (MQ=255)
         gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGc                                                                                                                                   <  1:60528/149‑1 (MQ=255)
         gcgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACtgtg                                                                                                                                    <  3:222284/148‑1 (MQ=37)
          cgcTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCtgtg                                                                                                                                    >  6:382843/1‑147 (MQ=255)
            ctctTTCTGTCTCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGt                                                                                                                                <  1:334045/149‑1 (MQ=37)
               ttGCTGGCGTTCAGGAGTCAATTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                                                               <  1:137123/147‑1 (MQ=255)
               ttGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                                                               <  5:382843/147‑1 (MQ=255)
                   tGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGGAAAATGTGTCGCCTG‑CTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                                                                                 >  6:77252/1‑140 (MQ=17)
                   tGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGGAAAATGTGTCGCCTG‑CTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                                                                                 <  5:77252/140‑1 (MQ=17)
                        cTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAACCCTGTGCAGTTGGACTGTGGGa                                                                                                                    <  7:248862/149‑1 (MQ=255)
                               gTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGc                                                                                                              <  4:134428/148‑1 (MQ=255)
                                tCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCAAAGTACAACTGCACAGTCTTGGAAGCCGGTGCAGATTGCCTTTGGGCGCGGGCGt                                                                                                            >  7:202940/1‑149 (MQ=255)
                                   ttttAATGGCTACCGGAAACAGATCTTTTCTCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTTTGCAGTTGGACTGTGGGAGCGGGCGt                                                                                                            <  2:61540/146‑1 (MQ=255)
                                   ttttAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCgg                                                                                                                <  5:358183/142‑1 (MQ=255)
                                   ttttAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCgg                                                                                                                >  6:358183/1‑142 (MQ=255)
                                   ttttAATGGCTACCGGAAACAGATATTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCAACGGACTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGt                                                                                                            >  1:61540/1‑146 (MQ=255)
                                             tACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCCAAGTCCAACTGCACACGCGTGGAAGCCAGGGCAGTTGGACGGTGGGAGCGGGCGTGCCCGCGAAgggg                                                                                               >  5:215233/1‑149 (MQ=255)
                                                  gAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAg                                                                                                                                 <  8:62406/110‑1 (MQ=255)
                                                  gAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAg                                                                                                                                 >  7:62406/1‑110 (MQ=255)
                                                         aTCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGGGAAGGGGCCGGAGCAGGCg                                                                                   >  7:216990/1‑149 (MQ=255)
                                                                  caaGGTCTGTTAAAATTGTGTTGCCACCCACCGGCCTCTTCGCTGGAACGCCCGCTCCCACAGTCCAACTGCAGAGTCTTGGATGCCGTTGCAGTTGTAGTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGtt                                                                          <  8:50/147‑1 (MQ=25)
                                                                                                      gctTCGCGGCCACGCCCGCCCCCACAGTCCAACGCCACTGGATTGTCCGCCTGTGCAGTTGGACGGTGGGTGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGaaa                                      <  6:215233/148‑1 (MQ=255)
                                                                                                          cGCTGGCACGCCCGCTCCCACATTCCAACTGCAGAGTGTTGGATGCCTGTGCAGTGGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAacac                                   <  8:202940/148‑1 (MQ=37)
                                                                                                                     ccGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGcctc                       >  7:165087/1‑149 (MQ=255)
                                                                                                                           ccACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGc                   <  8:165087/147‑1 (MQ=255)
                                                                                                                            cacaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAg                >  2:178132/1‑149 (MQ=255)
                                                                                                                             acaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAgg               >  8:230532/1‑149 (MQ=255)
                                                                                                                             acaGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGAGTGCCCGCGAAGGGGCAGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATAACCCACCAGAAACACACGTTCTGCCTCATGCCAgg               >  2:352618/1‑149 (MQ=38)
                                                                                                                                    aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGATACACCCGTTCTGCCTCCTGCCAGGCACCATc        >  4:307230/1‑149 (MQ=255)
                                                                                                                                    aaCTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc        >  4:324397/1‑149 (MQ=255)
                                                                                                                                    aaCTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGCGCAGGCGACATCACTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATc        >  8:228147/1‑149 (MQ=37)
                                                                                                                                        gCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCa    >  4:341551/1‑149 (MQ=255)
                                                                                                                                         cacaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCCTCGGCAt   >  1:347456/1‑149 (MQ=255)
                                                                                                                                          acaGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATc  >  8:306140/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACACCCGTTCTGCCTCCTGCCAGGCACCATCGGCATC  >  NC_002947/2282628‑2282914

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: ‑