Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,743,460 C→T 10.6% intergenic (‑237/+10) fbp ← / ← glgP class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,743,4600CT10.6% 125.9 / 5.9 47intergenic (‑237/+10)fbp/glgPclass 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (23/19);  new base T (3/2);  total (26/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.95e-01

CTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGA  >  NC_002947/5743312‑5743600
                                                                                                                                                    |                                                                                                                                            
ctctCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACg                                                                                                                                             >  1:465731/1‑150 (MQ=255)
cgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccattcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACg                                                                                                                                             >  1:61225‑M1/81‑150 (MQ=255)
      gACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCGGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctca                                                                                                                                       >  2:355557/1‑150 (MQ=255)
       aCGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctcac                                                                                                                                      <  2:293235/150‑1 (MQ=255)
           accattcgcgaatatgccaccgagatctggaaggcattggattgaggtgtgtccattcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctaacctca                                                                                                                                  >  2:562151‑M1/70‑150 (MQ=255)
               aTAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATcc                                                                                                                              <  1:545376/150‑1 (MQ=255)
                       gTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGCAGGGCGTCTTTGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCtt                                                                                                                      >  2:172712/1‑150 (MQ=255)
                          gTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCa                                                                                                                   <  1:368408/150‑1 (MQ=255)
                                  ttCGCCAGGACGGTGGAATGGTTGTGTGTCCTTTCGCGGGCGCGCTCGCTTCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCgg                                                                                                           <  2:465731/150‑1 (MQ=39)
                                      ccAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCaa                                                                                                                            <  1:148587/129‑1 (MQ=255)
                                      ccAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCaa                                                                                                                            >  2:148587/1‑129 (MQ=255)
                                                   cTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATg                                                                                          <  1:123497/150‑1 (MQ=255)
                                                     tgaggttaggccgttcgcgggcacgcccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGACACACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                                                                        >  2:562151‑M2/82‑150 (MQ=255)
                                                           gtgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTc                                                                                  <  1:355557/150‑1 (MQ=255)
                                                            tgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGACTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTcc                                                                                 >  2:444410/1‑150 (MQ=255)
                                                               ccATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                                                              >  1:122652/1‑150 (MQ=255)
                                                                cgttcgcgggcacgcccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgacg                                                                             >  1:61225‑M2/71‑150 (MQ=255)
                                                                cATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgacg                                                                             >  2:535144/1‑150 (MQ=255)
                                                                  ttCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGACACACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACgag                                                                           >  2:503762/1‑150 (MQ=255)
                                                                         gcacgcccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATg                                                                                          >  1:539277‑M2/62‑128 (MQ=255)
                                                                         gcacgcccgctcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATg                                                                                          <  2:539277‑M2/67‑1 (MQ=255)
                                                                                 GCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGCGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgt                                                            >  1:324467‑M1/1‑150 (MQ=255)
                                                                                   tcccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgAGCGCGCCCGCGAATGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCGGccgtgc                                                          >  1:485894‑M2/52‑150 (MQ=255)
                                                                                     ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTcc                                                                                 >  2:63961‑M2/50‑125 (MQ=255)
                                                                                     ccacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTcc                                                                                 <  1:63961‑M2/76‑1 (MQ=255)
                                                                                     CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                                                               >  1:144653‑M1/1‑111 (MQ=255)
                                                                                     CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                                                               <  2:144653‑M1/111‑1 (MQ=255)
                                                                                      cacaggtgtctcgcaagcttcaaggcttgcgctgtacctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccgt                                                       >  2:238108‑M2/49‑150 (MQ=255)
                                                                                         AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGACACACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt                                                    >  1:503762‑M1/1‑150 (MQ=255)
                                                                                                 cgcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCCTATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCATGCCGGGTTCGATACCg                                            >  2:108824/1‑150 (MQ=255)
                                                                                                       ccTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCTACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCg                                            <  1:423341/144‑1 (MQ=255)
                                                                                                       ccTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCg                                            >  2:423341/1‑144 (MQ=255)
                                                                                                        cTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcga                                                                               <  1:367588/108‑1 (MQ=255)
                                                                                                        cTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcga                                                                               >  2:367588/1‑108 (MQ=255)
                                                                                                         ttGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgcgc                                    >  1:183971/1‑150 (MQ=255)
                                                                                                            aggCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGccac                                 <  1:495344/148‑1 (MQ=255)
                                                                                                            aaGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAACGATGCCCTCCGATGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGccac                                 <  1:350456/150‑1 (MQ=255)
                                                                                                             aGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc                                                         >  2:221554/1‑125 (MQ=255)
                                                                                                             aGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgcc                                                         <  1:221554/125‑1 (MQ=255)
                                                                                                                    gagaCACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTg                         >  1:34623/1‑150 (MQ=255)
                                                                                                                       acacCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCgg                      <  2:350900/150‑1 (MQ=255)
                                                                                                                        tacctgtgggagcgAGCGCGCCCGCGAATGGACACACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgt                     <  1:562151‑M2/136‑1 (MQ=255)
                                                                                                                          cctgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATTCCGCCTTGCGCCACCACTCCTGCGgtgt                   <  1:108824‑M2/138‑1 (MQ=255)
                                                                                                                            tgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATAc                                              <  2:426566‑M2/111‑1 (MQ=255)
                                                                                                                            tgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATAc                                              >  1:426566‑M2/11‑121 (MQ=255)
                                                                                                                            tgtgggagcgGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGgtgtgt                 <  2:122652‑M2/140‑1 (MQ=255)
                                                                                                                                           gCCCGCGAATGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGACCCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGa  <  2:61225/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCACCACTCCTGCGGTGTGTGCCACAGATCCTCGA  >  NC_002947/5743312‑5743600

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