Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,777,217 C→G 8.9% L512V (CTG→GTG)  betA‑II → choline dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,777,2170CG 8.9% 134.6 / 2.6 43L512V (CTG→GTG) betA‑IIcholine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (21/18);  new base G (2/2);  total (23/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.51e-01

CGAAACCGCCTTCCACCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGA  >  NC_002947/5777093‑5777358
                                                                                                                            |                                                                                                                                             
cGAAACCGCCTTCCACCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAAc                                                                                                                              <  8:102428/142‑1 (MQ=255)
cGAAACCGCCTTCCACCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAAc                                                                                                                              >  7:102428/1‑142 (MQ=255)
      cGCCTTCCACCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGccacca                                                                                                                 >  7:27545/1‑149 (MQ=255)
      cGCCTTCCACCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGcc                                                                                                                     >  8:165413/1‑145 (MQ=255)
      cGCCTTCCACCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGcc                                                                                                                     <  7:165413/145‑1 (MQ=255)
       gCCTTCCACCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCcaccac                                                                                                                <  8:27545/149‑1 (MQ=255)
              aCCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCAt                                                                                                          >  6:27748/1‑148 (MQ=255)
              aCCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTACACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCAt                                                                                                          >  1:277863/1‑148 (MQ=255)
                 cGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATc                                                                                                      <  5:27748/149‑1 (MQ=255)
                           gTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCgag                                                                                                >  8:2691/1‑145 (MQ=255)
                               aaGATGGGCACCGACGACATGGCAGTGGTAGATGGTGAATGCCGCGTGCACGGCATTAAGGGCCTGCGCGTGGTCGATTCGTCGATCATGCCGGTGATCATCACCTTATACCTAAACGCCACCACGATCATGATCGCCGAGAAGATTTc                                                                                        <  5:334944/149‑1 (MQ=255)
                                aGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGcc                                                                                                   <  6:356730/137‑1 (MQ=255)
                                aGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGcc                                                                                                   >  5:356730/1‑137 (MQ=255)
                                 gATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGAGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCg                                                                                       >  1:74440/1‑148 (MQ=255)
                                                             gATGGTTAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATc                                                                                                            >  6:214394/1‑99 (MQ=255)
                                                             gATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGcc                                                          >  7:210140/1‑149 (MQ=255)
                                                                   gAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCgcag                                                    <  7:187233/149‑1 (MQ=255)
                                                                   gAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCgcag                                                    >  8:130147/1‑149 (MQ=255)
                                                                     aGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCg                                                                                                  <  4:290725/101‑1 (MQ=255)
                                                                     aGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCg                                                                                                  >  3:290725/1‑101 (MQ=255)
                                                                      ggCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCAc                                                 <  2:141877/149‑1 (MQ=255)
                                                                         cgcgTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGc                                              <  7:130147/149‑1 (MQ=255)
                                                                                     aTGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGc                                                                     >  4:110385/1‑114 (MQ=255)
                                                                                          gggCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCgg                             <  8:210140/149‑1 (MQ=255)
                                                                                          gggCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGAGCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCACGTAGTACTTGGCgg                             <  2:277863/149‑1 (MQ=255)
                                                                                                   cgTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCAc                    <  5:226565/149‑1 (MQ=255)
                                                                                                        tCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAGGTACTACGAGGCGGGGGATGCACCAGt                >  5:43999/1‑148 (MQ=255)
                                                                                                           aTGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGc                                                                     <  7:183370/92‑1 (MQ=255)
                                                                                                           aTGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGc                                                                     >  8:183370/1‑92 (MQ=255)
                                                                                                              ccTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGt                <  6:319537/142‑1 (MQ=255)
                                                                                                              ccTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGt                >  5:319537/1‑142 (MQ=255)
                                                                                                               cTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCa                  >  8:199877/1‑139 (MQ=255)
                                                                                                               cTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCa                  <  7:199877/139‑1 (MQ=255)
                                                                                                               cTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGc         >  3:222985/1‑148 (MQ=255)
                                                                                                               cTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGACGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGc         >  1:204558/1‑148 (MQ=255)
                                                                                                                tCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGt                <  6:230239/140‑1 (MQ=255)
                                                                                                                tCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGt                >  5:230239/1‑140 (MQ=255)
                                                                                                                  gATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAGGTACTACGCGGCGGGGGATGCACCAGTGAAAGGCaa       <  6:43999/147‑1 (MQ=255)
                                                                                                                   aTCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGcc    >  3:257394/1‑149 (MQ=255)
                                                                                                                   aTCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGcc    <  4:222985/149‑1 (MQ=255)
                                                                                                                     gttGCCGGTGATGATCAGCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGa  >  3:320200/3‑149 (MQ=255)
                                                                                                                      ttGCCGGTGATGATCAGCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCAcc                                                >  4:345212/2‑102 (MQ=255)
                                                                                                                      ttGCCGGTGATGATCAGCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCAcc                                                <  3:345212/101‑1 (MQ=255)
                                                                                                                            |                                                                                                                                             
CGAAACCGCCTTCCACCCGTCCTGCTCGTGCAAGATGGGCACCGACGACATGGCAGTGGTCGATGGTGAAGGCCGCGTGCACGGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGA  >  NC_002947/5777093‑5777358

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