Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,731,892 C→T 11.6% A464A (GCG→GCA PP_3299 ← lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,731,8920CT11.6% 108.0 / 5.6 43A464A (GCG→GCAPP_3299lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (18/20);  new base T (3/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

TACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTT  >  NC_002947/3731750‑3732040
                                                                                                                                              |                                                                                                                                                     
tACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCGCGTgctgcta                                                                                                                                                <  1:347560/150‑2 (MQ=255)
    cgGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGt                                                                                                                                            <  2:664385/150‑1 (MQ=255)
                  tCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTc                                                                                                                              <  1:422917/150‑1 (MQ=255)
                     tCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGa                                                                                                                           >  2:542939/1‑150 (MQ=255)
                       aaCGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGaaa                                                                                                                        >  1:38010/1‑150 (MQ=255)
                          gAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACgcgc                                                                                                                                                     <  1:705912/119‑1 (MQ=255)
                          gAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACgcgc                                                                                                                                                     >  2:705912/1‑119 (MQ=255)
                             ggTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTc                                                                                                                              >  2:501660/1‑139 (MQ=255)
                             ggTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTc                                                                                                                              <  1:501660/139‑1 (MQ=255)
                                      gTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTgaga                                                                                                         >  2:525706/1‑150 (MQ=255)
                                       tcttctGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGaa                                                                                                        <  1:525706/150‑1 (MQ=255)
                                                 cAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGAt                                                                                                                  >  1:824330/1‑130 (MQ=255)
                                                 cAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGAt                                                                                                                  <  2:824330/130‑1 (MQ=255)
                                                  aaGGCATAAACGCAGTCCCACCATCGATATTCAAAATCTGAGTCAGTTGACCACCCGACACCCAGGGCGAGGTAGTCAGCAATAGCAGCACGCGCTGCTTCAGTTTGAGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCt                                                                                             <  2:158919/150‑1 (MQ=255)
                                                    ggCATTTACGCAGTCCCATCATCGCTGTGCCACATCCGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCGTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGc                                                                                           <  1:107503/150‑1 (MQ=255)
                                                             gCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGGGTAGCTGCTGATCTGCGTTCAGAAc                                                                                  <  1:542939/150‑1 (MQ=255)
                                                                gTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTAGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTc                                                                               <  1:452962/150‑1 (MQ=255)
                                                                  cccATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCCATAGCAGCACGCGCTGCTTCCGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTa                                                                             >  1:841697/1‑150 (MQ=255)
                                                                        atcGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCt                                                                                                             <  2:87260/112‑1 (MQ=255)
                                                                        atcGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCt                                                                                                             >  1:87260/1‑112 (MQ=255)
                                                                        atcGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAAc                                                                       >  1:848763/1‑150 (MQ=255)
                                                                         tcGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACt                                                                      <  2:38010/150‑1 (MQ=255)
                                                                                 ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTacac                                                              >  2:648958/1‑150 (MQ=255)
                                                                                 ccACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTAGCACCAACGGCCTTTCAGGAATGTGCTGATCTGGGTTCAGAACGTCTATCAAACTGGCTacac                                                              >  2:239850/1‑150 (MQ=255)
                                                                                  cacaTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGt                                                                                <  2:516369/131‑1 (MQ=255)
                                                                                  cacaTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGt                                                                                >  1:516369/1‑131 (MQ=255)
                                                                                     aTCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGTCTACACCGtt                                                          >  1:498151/1‑150 (MQ=255)
                                                                                       cTGAGTCAGTTGACCACCCGCCAGCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATATCCTAAGAAGCTGCGGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCt                                                        <  2:292395/150‑1 (MQ=255)
                                                                                       cTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCt                                                        <  2:51075/150‑1 (MQ=255)
                                                                                                    ccaccCGCCACCCAGGGCGACGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGaa                                           <  1:648958/150‑1 (MQ=255)
                                                                                                        ccGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGAc                                       >  2:829157/1‑150 (MQ=255)
                                                                                                            cACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAc                                   <  1:725156/150‑1 (MQ=255)
                                                                                                             aCCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcg                                  >  2:811089/1‑150 (MQ=255)
                                                                                                              cccAGGGCGAGGTAGGCAGCAATAGGCATGTTCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAAcgt                                 >  1:663291/1‑150 (MQ=255)
                                                                                                                      gAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATg                                                  <  2:15075/125‑1 (MQ=255)
                                                                                                                      gAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATg                                                  >  1:15075/1‑125 (MQ=255)
                                                                                                                            ggCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAg                   >  1:528174/1‑150 (MQ=255)
                                                                                                                                actgaagcagcGCGTGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGAt               >  1:412854/6‑150 (MQ=255)
                                                                                                                                 ctgaagcagcGCGTGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTaca                                                               >  2:313412/5‑101 (MQ=255)
                                                                                                                                 ctgaagcagcGCGTGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTaca                                                               <  1:313412/97‑1 (MQ=255)
                                                                                                                                 ctgaagcagcGCGTGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATa              >  1:197762/5‑150 (MQ=255)
                                                                                                                                 caATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATa              <  2:528174/150‑1 (MQ=255)
                                                                                                                                  aaTAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGTCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATAc             <  2:498151/150‑1 (MQ=255)
                                                                                                                                           aCGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATAc             >  2:319207/1‑141 (MQ=255)
                                                                                                                                           aCGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATAc             <  1:319207/141‑1 (MQ=255)
                                                                                                                                             gcgcTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTACACAGCTTCAGATCTTCTTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGtt  <  2:663291/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
TACGCGGTATGAGCCCAATCCTCAACGAAGGTGAGTACGTCTTCTGCACCAAGGCATTTACGCAGTCCCATCATCGCTGTGCCACATCTGAGTCAGTTGACCACCCGCCACCCAGGGCGAGGTAGGCAGCAATAGCAGCACGCGCTGCTTCAGTTTGCGCAACTGCTCGGAAATCCGATGTCCTGAGAAGCTGCTGATCTGCGTTCAGAACGTCTATCAAACTGGCTACACCGTTCTGGAATGCTACGAAGGACGAACGTCGTGCTCGCTCCAGGGATACCTCACCGCGTT  >  NC_002947/3731750‑3732040

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: