Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,282,780 T→G 18.1% 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,7800TG18.1% 65.0 / 7.6 32intergenic (‑129/+100)PP_2010/PP_2011cytochrome b561/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (10/16);  new base G (2/4);  total (12/20)
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.86e-01

GTGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACAC  >  NC_002947/2282632‑2282924
                                                                                                                                                    |                                                                                                                                                
gtggGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTg                                                                                                                                                 <  1:643130/150‑1 (MQ=255)
 tggGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTgg                                                                                                                                                >  2:128149/1‑150 (MQ=255)
     gggggATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACg                                                                                                                                            >  1:504381/1‑150 (MQ=255)
          aTGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGc                                                                                                                                       >  2:614232/1‑150 (MQ=255)
               aCCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCca                                                                                                                                  <  1:614232/150‑1 (MQ=255)
               aCCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCca                                                                                                                                  >  2:113095/1‑150 (MQ=255)
                 cATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCcaca                                                                                                                                >  1:1073226/1‑150 (MQ=255)
                         gTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACt                                                                                                                        <  1:352883/150‑1 (MQ=255)
                         gTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGATCACGCCCGCTCCCACAGTCCAACt                                                                                                                        <  1:1068996/150‑1 (MQ=255)
                            cTCGTCATCAGGTCTGAATGGGGCTCTTGCTGGGGCTCCGGAGTCATTTTATTTGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCCACGGCATATTCGCTGGCACGCCCGCTCCTACAGTCCATCTGca                                                                                                                     <  2:676846/150‑1 (MQ=255)
                                   tcaGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCtt                                                                                                              <  1:128149/150‑1 (MQ=255)
                                   tcaGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAATGTGTTGCCAGCCACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttc                                                                                                             <  1:335046/150‑2 (MQ=38)
                                    caGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTg                                                                                                             <  2:1073226/150‑1 (MQ=255)
                                     aGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                                            <  1:295609/150‑3 (MQ=255)
                                          tGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttccaagac                                                                                                       >  1:1064742/1‑143 (MQ=255)
                                                   ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                              <  1:369382/150‑1 (MQ=255)
                                                   ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                              <  2:504381/150‑1 (MQ=255)
                                                   ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGtt                                                                                              <  2:1067943/150‑1 (MQ=255)
                                                   ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt                                                                                              <  2:1064742/150‑1 (MQ=37)
                                                   ctctTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGtt                                                                                              <  1:113095/150‑1 (MQ=25)
                                                            gcgcTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCTTCCAAGACTGTGCAGTTGGACTGTgg                                                                                     <  1:1145521/150‑1 (MQ=255)
                                                                     aGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCg                                                                            >  2:449266/1‑150 (MQ=255)
                                                                                  gCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAggg                                                               <  1:449266/150‑1 (MQ=255)
                                                                                       cGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                                            <  1:415980/100‑3 (MQ=255)
                                                                                       cGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGGCttcc                                                                                                            >  2:415980/1‑98 (MQ=255)
                                                                                                                  tgAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGa                               >  2:884940/1‑150 (MQ=255)
                                                                                                                      aaaTGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCaa                                                                                                                          >  1:1066862/1‑55 (MQ=255)
                                                                                                                      aaaTGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCaa                                                                                                                          <  2:1066862/55‑1 (MQ=255)
                                                                                                                                        cGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAg                                                                                                >  2:844067/1‑63 (MQ=11)
                                                                                                                                        cGGCCCCTTCGCGGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAg                                                                                                <  1:844067/63‑1 (MQ=14)
                                                                                                                                         ggCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAg        >  2:1014267/1‑150 (MQ=255)
                                                                                                                                               ttCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAacac  <  1:884940/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTGGGGGGGGATGAAACCATGATTGGTCCTCGTCATCAGGTCTGAATGGCGCTCTTGCTGGCGCTCAGGAGTCATTTTAATGGCTACCGGAAACAGATCTTTGCGCCAGGTCTGTGAAAAATGTGTTGCCAGCCACCGGCCTCTTCGCTGGCACGCCCGCTCCCACAGTCCAACTGCACAGTCTTGGAAGCCTGTGCAGTTGGACTGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGGAGCAGGCGACATCAGTTACCGGGCAGATAGGGCCAGTCATCACCCACCAGAAACAC  >  NC_002947/2282632‑2282924

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