Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,622,598 G→A 8.3% V211I (GTT→ATT)  U712_08340 → putative L‑serine dehydratase, alpha chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,622,5980GA 8.3% 150.0 / 2.2 48V211I (GTT→ATT) U712_08340putative L‑serine dehydratase, alpha chain
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base A (2/2);  total (25/23)
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

TGTCAGGCAGAGGTTGGATCAGCTTCAGGCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCGTATTCCATGTGATGAGGTTATCGATGCCATGTATAA  >  CP006881/1622463‑1622734
                                                                                                                                       |                                                                                                                                        
tGTCAGGCAGAGGTTGGATCAGCTTCAGGCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTg                                                                                                                                       >  5:303096/1‑139 (MQ=255)
  gccggCAGAGGTTGGATCAGCTTCAGGCATGGCGGCGGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGc                                                                                                                                      <  2:46389/135‑1 (MQ=255)
        agagGTTGGATCAGCTTCAGGCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCtt                                                                                                                               >  7:315857/1‑139 (MQ=255)
            gTTGGATCAGCTTCAGGCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCttgtt                                                                                                                            >  8:36992/1‑138 (MQ=255)
                   cAGCTTCAGGCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGcc                                                                                                                    <  4:17779/139‑1 (MQ=255)
                      cTTCAGGCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACCCCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTg                                                                                                                 >  2:312919/1‑139 (MQ=255)
                      cTTCAGGCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTg                                                                                                                 >  6:56471/1‑139 (MQ=255)
                           ggCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGa                                                                                                            <  5:56471/139‑1 (MQ=255)
                            gCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGaa                                                                                                           <  1:137159/139‑1 (MQ=255)
                            gCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGaa                                                                                                           <  2:199026/139‑1 (MQ=255)
                                ggcggcTGCGGCTATTGTTGAAATGGCAGGAGGACCACCCGAACAAAGCGCGGCAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTAGAGGGGCCTTGCGTGAAACg                                                                                                        >  6:283382/1‑138 (MQ=255)
                                ggcggcTGCGGCTATTGTTGAAATGGCAGGAGGAACCCCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACg                                                                                                        >  4:99431/1‑138 (MQ=255)
                                ggcggcTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACg                                                                                                        >  3:168468/1‑138 (MQ=255)
                                ggcggcTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACg                                                                                                        >  1:214312/1‑138 (MQ=255)
                                 gcggctgcggctATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACgg                                                                                                       >  2:124891/1‑138 (MQ=255)
                                               tgttgAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGcgcg                                                                                        >  5:331680/1‑139 (MQ=255)
                                                gttgAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGCCCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGt                                                                                       >  4:101899/1‑139 (MQ=255)
                                                    aaaTGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCaaa                                                                                   >  3:124777/1‑139 (MQ=255)
                                                     aaTGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGCCCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCTCGTCAAAt                                                                                  >  6:285226/1‑139 (MQ=255)
                                                               ggaACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATtg                                                                        <  6:303096/139‑1 (MQ=255)
                                                                gaACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATtgc                                                                       <  8:315857/139‑1 (MQ=255)
                                                                  ccacCCGAACAAAGCGCGGACGCCATGGCAATTACATTAAAAAATATGCTAGGCCTCGGGTGCGACCCTGTTGCAGGGGTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATtgctg                                                                     <  2:249503/138‑1 (MQ=255)
                                                                        gAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCa                                                                                     >  6:116661/1‑117 (MQ=255)
                                                                        gAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCa                                                                                     <  5:116661/117‑1 (MQ=255)
                                                                                cgcgGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTg                                                       >  6:349043/1‑139 (MQ=255)
                                                                                 gcgGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTgg                                                      <  1:312919/139‑1 (MQ=255)
                                                                                    gAAGCCATGGCCATTACATTAAAAAATATGCTAGTCCTCGTGTGCGACCGTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCgg                                                   <  5:349043/139‑1 (MQ=255)
                                                                                           tGGCAATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTAc                                            <  1:124891/139‑1 (MQ=255)
                                                                                               cattacattaAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGc                                        <  7:179430/139‑1 (MQ=255)
                                                                                                 ttacattaAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCg                                      >  5:195396/1‑139 (MQ=255)
                                                                                                  tacattaAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCTGGGCTTGTTGAGGTGCCTTGCGTGAAACGGATTGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACATGCCGt                                     >  2:158282/1‑139 (MQ=255)
                                                                                                       taAAAAATATGCTAGGCCTCGTGTGCGCCCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCGTATTc                                 >  3:173038/1‑138 (MQ=255)
                                                                                                       taAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGcc                                                                                                                    <  8:139214/55‑1 (MQ=255)
                                                                                                       taAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGcc                                                                                                                    >  7:139214/1‑55 (MQ=255)
                                                                                                        aaaaaaTATGCTAGGGCTCGTGGGCGACCCTGTTGCAGGGCTTGTAGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCGTATTcc                                <  1:229459/138‑1 (MQ=255)
                                                                                                        aaaaaaTATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGGGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCGTATTcc                                <  8:22762/138‑1 (MQ=255)
                                                                                                              tatGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGCTATGGCATTGGCGGGCATTACAAGCCGTATTCCATGtgat                         >  2:163350/1‑139 (MQ=255)
                                                                                                              tatGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCGTATTCCATGtgat                         <  4:124777/139‑1 (MQ=255)
                                                                                                              tatGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCGTATTCCATGtgat                         <  7:36992/139‑1 (MQ=255)
                                                                                                                           gtgtGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGata                                                               <  2:252272/88‑1 (MQ=255)
                                                                                                                           gtgtGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGata                                                               >  1:252272/1‑88 (MQ=255)
                                                                                                                             gtgCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGTATTGGCGGGCATTACAAGCCGTATTCCATGTGATGAGGTTATCGATGcc          <  2:214312/139‑1 (MQ=255)
                                                                                                                              tgCGACCCTGTTGCAGGGGTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCGTATTCCATGTGATGAGGTTATCGATGCCa         <  3:99431/139‑1 (MQ=255)
                                                                                                                                  aCCCTATTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACgg                                                                                                       <  4:36382/41‑1 (MQ=255)
                                                                                                                                  aCCCTATTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACgg                                                                                                       >  3:31581/1‑41 (MQ=255)
                                                                                                                                  aCCCTATTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACgg                                                                                                       >  3:36382/1‑41 (MQ=255)
                                                                                                                                  aCCCTATTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACgg                                                                                                       <  4:31581/41‑1 (MQ=255)
                                                                                                                                      tGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCGTATTCCATGTGATGAGGTTATCGATGCCATGTATaa  >  1:137534/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                                        
TGTCAGGCAGAGGTTGGATCAGCTTCAGGCATGGCGGCTGCGGCTATTGTTGAAATGGCAGGAGGAACACCCGAACAAAGCGCGGAAGCCATGGCCATTACATTAAAAAATATGCTAGGCCTCGTGTGCGACCCTGTTGCAGGGCTTGTTGAGGTGCCTTGCGTGAAACGGAATGCGATGGGCGCGTCAAATGCGATGATTGCTGCTGATATGGCATTGGCGGGCATTACAAGCCGTATTCCATGTGATGAGGTTATCGATGCCATGTATAA  >  CP006881/1622463‑1622734

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