Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,279,705 G→A intergenic (+1373/+6659) mtlD → / ← glmM mannitol‑1‑phosphate 5‑dehydrogenase/phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,279,7050GA96.9% 97.2 / NA 32intergenic (+1373/+6659)mtlD/glmMmannitol‑1‑phosphate 5‑dehydrogenase/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/17);  minor base C (0/1);  total (14/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 = 1.00e+00

TTGTGCCGTTATTAGTTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTGCTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCATTTGAATTT  >  CP000730/2279568‑2279836
                                                                                                                                         |                                                                                                                                   
ttGTGCCGTTATTAGTTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACtg                                                                                                                                  >  1:225873/1‑141 (MQ=255)
  gtgCCGTTATTAGTTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACtgtt                                                                                                                                <  1:39744/141‑1 (MQ=255)
         tattaGTTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCAc                                                                                                                         <  2:126938/141‑1 (MQ=255)
               ttGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTg                                                                                                                   >  2:102082/1‑141 (MQ=255)
                  aTTCTACTTCTGGTTTACTAGTTACCTCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCg                                                                                                                <  1:47032/141‑1 (MQ=255)
                    tCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCa                                                                                                              <  1:246216/141‑1 (MQ=255)
                        cTTCTGGTTTACTAGTTACTTCGTTATCCATTGTCGTACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTGGTTTCCCTTGTTGCCGCATCAc                                                                                                          <  2:72200/141‑1 (MQ=255)
                        cTTCTGGTTTACTAGTTACATCGTTATCTATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGTTGTTCCTGTGGTTTCGTCTGTTGCCGCATCAc                                                                                                          <  2:72186/141‑1 (MQ=255)
                           cTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACtat                                                                                                       <  1:152782/141‑1 (MQ=255)
                                   cTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGgt                                                                                               >  1:110450/1‑141 (MQ=255)
                                             cGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGgag                                                                                     <  2:20662/141‑1 (MQ=255)
                                              gTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGgaga                                                                                    <  1:226631/141‑1 (MQ=255)
                                               ttATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCGGCATCACTATTATTTGGTGTTGTCGGAGaa                                                                                   >  1:221684/1‑141 (MQ=255)
                                                        gTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCt                                                                          <  1:25149/141‑1 (MQ=255)
                                                                    gttgATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCt                                                              <  2:248354/141‑1 (MQ=255)
                                                                     ttgATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTg                                                             >  2:90999/1‑141 (MQ=255)
                                                                          gCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTgtcgtc                                                        >  1:12875/1‑141 (MQ=255)
                                                                             tCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCa                                                     <  1:54657/141‑1 (MQ=255)
                                                                                   cTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTTACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACg                                               >  2:120385/1‑141 (MQ=255)
                                                                                      gAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTa                                            >  1:178357/1‑141 (MQ=255)
                                                                                      gAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTa                                            <  1:160700/141‑1 (MQ=255)
                                                                                       aaTTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACCTTAg                                           >  2:99252/1‑141 (MQ=255)
                                                                                           gTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTg                                       >  1:186578/1‑141 (MQ=255)
                                                                                             tattaGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTc                                     >  1:183737/1‑141 (MQ=255)
                                                                                                      gCGGTTTATCATGTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATGTGGTGTTGTCTGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTAtttt                            <  2:85375/141‑1 (MQ=255)
                                                                                                         gTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTg                         <  2:187565/141‑1 (MQ=255)
                                                                                                            tATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATGTGGTGAGGTCGGAGGAGCGTCTGCTTGGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGgtg                      >  1:100048/1‑141 (MQ=255)
                                                                                                                aTTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGc                  <  2:225873/141‑1 (MQ=255)
                                                                                                                       tcatcaGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCa           >  1:95450/1‑141 (MQ=255)
                                                                                                                         atcaGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCAtt         <  2:193472/141‑1 (MQ=255)
                                                                                                                          tcaGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCAttt        >  1:31212/1‑141 (MQ=255)
                                                                                                                                gCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGACGCAGCGTCTGCTGTGCCATTAGCTGTCGTCGCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCATTTGAAttt  <  1:11903/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                   
TTGTGCCGTTATTAGTTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTGCTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCATTTGAATTT  >  CP000730/2279568‑2279836

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: