Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,282,737 T→G 16.3% intergenic (‑129/+100) 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,7370TG16.3% 79.6 / 9.4 36intergenic (‑129/+100)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/17);  new base G (2/4);  total (15/21)
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.57e-01

ATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACAC  >  NC_002947/2282599‑2282881
                                                                                                                                          |                                                                                                                                                
aTGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTg                                                                                                                                                 >  4:305940/1‑140 (MQ=255)
aTGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTg                                                                                                                                                 <  3:305940/140‑1 (MQ=255)
aTGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGccc                                                                                                                                         >  4:222284/1‑148 (MQ=255)
    aaCCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACAGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTcc                                                                                                                                    >  2:60528/1‑149 (MQ=255)
       cATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCca                                                                                                                                  >  5:32408/1‑148 (MQ=255)
                  cTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTg                                                                                                                       <  1:188171/148‑1 (MQ=255)
                  cTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGca                                                                                                                     <  5:353756/149‑1 (MQ=37)
                             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)
                                          tctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTg                                                                                                                                                 >  1:17987/1‑98 (MQ=255)
                                          tctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTg                                                                                                                                                 <  2:17987/98‑1 (MQ=255)
                                            ttGCTGGCGTTCAGGAGTCAATTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                              <  1:137123/147‑1 (MQ=255)
                                            ttGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                              <  5:382843/147‑1 (MQ=255)
                                                tGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGGAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                                                <  5:77252/140‑1 (MQ=17)
                                                tGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGGAAAATGTGTCGCCTGCTCCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAg                                                                                                >  6:77252/1‑140 (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                                                                                                >  7:62406/1‑110 (MQ=255)
                                                                               gAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAg                                                                                                <  8:62406/110‑1 (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)
                                                                                                                                          |                                                                                                                                                
ATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACAC  >  NC_002947/2282599‑2282881

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