Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,346,315 C→A 14.7% V555V (GTC→GTA acsA‑II → acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,346,3150CA14.7% 84.2 / 6.4 34V555V (GTC→GTAacsA‑IIacetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base A (2/3);  total (16/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 = 5.23e-01

CTGGATCACCGGGCGCATCGACGATGTGATCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAG  >  NC_002947/5346168‑5346460
                                                                                                                                                   |                                                                                                                                                 
ctGGATCACCGGGCGCATCGACGATGTGATCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTa                                                                                                                                                 >  2:242712/1‑150 (MQ=255)
 tGGATCACCGGGCGCATCGACGATGTGATCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTAc                                                                                                                                                <  1:117465/150‑1 (MQ=255)
 tGGATCACCGGGCGCATCGACGATGTGATCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACATCAAGGGCCAAGGCGTCTAc                                                                                                                                                >  2:117465/1‑150 (MQ=255)
       aCCGGGCGCATCGACGATGTGATCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGACAGTGAACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCGACCCCCATGACCTCAAGGGCCAAGGCGTAGACGCCttg                                                                                                                                          <  2:186387/150‑2 (MQ=255)
       aCCGGGCGCATCGACGATGTGATCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTc                                                                                                                                          <  1:141894/150‑1 (MQ=255)
        ccGGGCGCATCGACGATGTGATCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCg                                                                                                                                         <  1:242712/150‑1 (MQ=255)
                         gtgATCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGc                                                                                                                        >  2:268668/1‑150 (MQ=255)
                          tgATCAACGTGTCCGGGCACCGCATCGGCAACGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGATACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCg                                                                                                                       >  1:506286/1‑150 (MQ=255)
                             tCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGa                                                                                                                    <  1:408691/150‑1 (MQ=255)
                                aCGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGacac                                                                                                                 >  1:147062/1‑150 (MQ=255)
                                   tgtCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAg                                                                                                              <  1:121251/150‑1 (MQ=255)
                                     tCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGAt                                                                                                            >  1:28912/1‑150 (MQ=255)
                                            aCCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTAGACGCCttg                                                                                                                                          <  1:302274/113‑2 (MQ=255)
                                            aCCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTAGACGCCttg                                                                                                                                          >  2:302274/1‑112 (MQ=255)
                                                  tCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAgg                                                                                               <  1:127257/150‑1 (MQ=255)
                                                     gCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCg                                                                                            >  2:461632/1‑150 (MQ=255)
                                                                             tGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTc                                                                                                                                          <  2:92129/80‑1 (MQ=255)
                                                                             tGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTc                                                                                                                                          >  1:92129/1‑80 (MQ=255)
                                                                                 gctgcATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAAt                                                                >  1:21093/1‑150 (MQ=255)
                                                                                  ctgcATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATc                                                               <  2:21093/150‑1 (MQ=255)
                                                                                                   cGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGcc                                              >  2:375686/1‑150 (MQ=255)
                                                                                                        ccGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGc                                         <  2:315990/150‑1 (MQ=255)
                                                                                                           cgGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGa                                      <  1:202082/150‑1 (MQ=255)
                                                                                                               ggtCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCa                                  <  2:350368/150‑1 (MQ=255)
                                                                                                                 tCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGt                                >  2:228121/1‑150 (MQ=255)
                                                                                                                  cGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCGTCGCCAAGCCTGAGCTGATTCAGTg                               <  1:268668/150‑1 (MQ=255)
                                                                                                                         ccccACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCa                        >  2:112807/1‑150 (MQ=255)
                                                                                                                                  cTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTAccc               <  2:366647/150‑1 (MQ=255)
                                                                                                                                   tCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGa                                      >  1:392687/1‑126 (MQ=255)
                                                                                                                                   tCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGa                                      <  2:392687/126‑1 (MQ=255)
                                                                                                                                   tCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATAAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCa              >  1:412203/1‑150 (MQ=255)
                                                                                                                                         acgaaGGCGTAGACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCaa             >  1:78409/4‑145 (MQ=255)
                                                                                                                                         acgaaGGCGTAGACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCaa             <  2:78409/142‑1 (MQ=255)
                                                                                                                                               gCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAg  <  1:464962/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTGGATCACCGGGCGCATCGACGATGTGATCAACGTGTCCGGGCACCGCATCGGCACCGCTGAAGTGGAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAG  >  NC_002947/5346168‑5346460

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: