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,7050GA97.3% 114.9 / ‑6.3 37intergenic (+1373/+6659)mtlD/glmMmannitol‑1‑phosphate 5‑dehydrogenase/phosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/20);  minor base C (0/1);  total (16/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.94e-01

TTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTGCTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCATTTGAATTTGC  >  CP000730/2279583‑2279838
                                                                                                                          |                                                                                                                                     
ttGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTg                                                                                                                     <  1:290579/141‑1 (MQ=255)
   aTTCTACTTCTGGTTTACTAGTGACATCGTTATCCATTGTCAGACTGTTTGTTGCTGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCg                                                                                                                  <  2:93101/141‑1 (MQ=255)
        aCTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGTCTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCa                                                                                                             <  2:21299/141‑1 (MQ=255)
        aCTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCa                                                                                                             <  2:299585/141‑1 (MQ=255)
           tCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACta                                                                                                          >  1:144721/1‑141 (MQ=255)
                   aCTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCCCCTGTTGCCGCGTCAATATTATTTgg                                                                                                  >  2:197149/1‑141 (MQ=255)
                   aCTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTgg                                                                                                  <  2:112537/141‑1 (MQ=255)
                       gTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGtgtt                                                                                              <  2:283569/141‑1 (MQ=255)
                       gTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTCTTGTTTCACCTGTTGCCGCATCACTATTATTTGGtgtt                                                                                              <  1:161233/141‑1 (MQ=255)
                        ttACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGtgttg                                                                                             >  1:66068/1‑141 (MQ=255)
                        ttACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGtgttg                                                                                             >  1:255395/1‑141 (MQ=255)
                            aTCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCgg                                                                                         <  1:166336/141‑1 (MQ=255)
                               gTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGgaga                                                                                      <  1:202239/141‑1 (MQ=255)
                                         gTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCt                                                                            >  2:246717/1‑141 (MQ=255)
                                               ttgtttgttTATGCATCTTCACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCTCCTGTTGCCGTATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCa                                                                      <  1:10176/140‑1 (MQ=255)
                                                       tgATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGCTGTTCCTGTTGTTTCTCCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTgt                                                              <  1:12643/141‑1 (MQ=255)
                                                            cATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGtct                                                         <  1:30456/141‑1 (MQ=255)
                                                                    cTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACg                                                 <  2:37249/141‑1 (MQ=255)
                                                                     tAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGt                                                <  1:231781/141‑1 (MQ=255)
                                                                      aGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGtt                                               >  1:94609/1‑141 (MQ=255)
                                                                        aaTTGTTATTAGCTTGCGGTTTATCAGTTGCAGCAGCAGTTGCTGATGCTACGGCTGCTGCACCTGTTTCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAg                                             <  1:197149/141‑1 (MQ=255)
                                                                               attaGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTcc                                      <  2:144721/141‑1 (MQ=255)
                                                                                ttaGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCa                                     >  1:9311/1‑141 (MQ=255)
                                                                                ttaGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCa                                     >  1:244743/1‑141 (MQ=255)
                                                                                       gCGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTAtttt                              <  2:253312/141‑1 (MQ=255)
                                                                                        cGGTTTATCATTTGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTc                             <  2:93597/141‑1 (MQ=255)
                                                                                                  tttGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCa                   <  2:4964/141‑1 (MQ=255)
                                                                                                  tttGCATCATCAGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCa                   >  1:249385/1‑141 (MQ=255)
                                                                                                  tttGCATCATCAGTTGCTGACGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCa                   >  1:184573/1‑141 (MQ=255)
                                                                                                         catcaGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCAt            <  1:145053/141‑1 (MQ=255)
                                                                                                          atcaGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCAtt           >  2:150603/1‑141 (MQ=255)
                                                                                                          atcaGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCAtt           >  1:222138/1‑141 (MQ=255)
                                                                                                           tcaGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCTCATCACTATTATTTGGTGTTGTCGGAGACGCGTCTGCTTTGCCATTAGCTGTCGTCTCGGCTACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCAGTAGCAttt          >  2:183996/1‑141 (MQ=255)
                                                                                                            caGTTGCTGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCATTTg         >  1:276580/1‑141 (MQ=255)
                                                                                                                   tGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCATTTGAATTTGc  >  2:124211/1‑141 (MQ=255)
                                                                                                                   tGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCATTTGAATTTGc  >  2:169560/1‑141 (MQ=255)
                                                                                                                   tGATGTTACTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCATTTGAATTTGc  <  2:263929/141‑1 (MQ=255)
                                                                                                                          |                                                                                                                                     
TTGATTCTACTTCTGGTTTACTAGTTACATCGTTATCCATTGTCGGACTGTTTGTTGATGCATCTACACTAGAATTGTTATTAGCTTGCGGTTTATCATTTGCATCATCAGTTGCTGATGTTGCTGTTGTTTCACCTGTTGCCGCATCACTATTATTTGGTGTTGTCGGAGAAGCGTCTGCTTTGCCATTAGCTGTCGTCTCAGATACGTTAGGTTGTCCAGTATTTTCTGGTGTTGCATTAGCATTTGAATTTGC  >  CP000730/2279583‑2279838

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: