Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,115,907 C→T *407* (TAG→TAA wcaL ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,115,9070CT96.7% 87.8 / ‑6.0 30*407* (TAG→TAAwcaLputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (13/16);  total (13/17)
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.98e-01

GGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCG  >  NC_000913/2115772‑2116041
                                                                                                                                       |                                                                                                                                      
gggTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTtat                                                                                                                                     <  2:674671/139‑1 (MQ=255)
 ggTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTtata                                                                                                                                    >  2:968799/1‑139 (MQ=255)
    tGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAg                                                                                                                                 >  2:119225/1‑139 (MQ=255)
           gTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAATGCAAGCGCCGAGCGTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgcag                                                                                                                          <  2:767842/139‑1 (MQ=255)
               aTGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAg                                                                                                                      <  1:433270/139‑1 (MQ=255)
                tGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGc                                                                                                                     >  2:345509/1‑139 (MQ=255)
                    aaCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGc                                                                                                                 <  2:804075/139‑1 (MQ=255)
                     aCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGc                                                                                                                 >  2:75226/1‑138 (MQ=255)
                     aCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGc                                                                                                                 <  1:75226/138‑1 (MQ=255)
                     aCGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCg                                                                                                                >  2:163948/1‑139 (MQ=255)
                          ttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCTGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgc                                                                                                                            <  2:926655/122‑1 (MQ=255)
                          ttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgc                                                                                                                            >  1:926655/1‑122 (MQ=255)
                          ttGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTtc                                                                                                           >  1:1300163/1‑139 (MQ=255)
                            gTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTtctc                                                                                                         <  1:163948/139‑1 (MQ=255)
                                         cgcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            <  1:968799/139‑1 (MQ=255)
                                         cgcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            <  1:345509/139‑1 (MQ=255)
                                         cgcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            <  2:926716/139‑1 (MQ=255)
                                         cgcgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCAGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctg                                                                                            <  2:941422/139‑1 (MQ=255)
                                           cgcgGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgct                                                                                          >  1:1211827/1‑139 (MQ=255)
                                                        ggTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCgtg                                                                             <  2:700899/139‑1 (MQ=255)
                                                          tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTtata                                                                                                                                    >  1:1252469/1‑82 (MQ=255)
                                                          tATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTtata                                                                                                                                    <  2:1252469/82‑1 (MQ=255)
                                                                           gcCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTcgcgcgcg                                                          <  2:1300163/139‑1 (MQ=255)
                                                                                           ggAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAg                                          <  2:1211827/139‑1 (MQ=255)
                                                                                                gTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgct                                                                                          <  2:841542/86‑1 (MQ=255)
                                                                                                gTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACctgct                                                                                          >  1:841542/1‑86 (MQ=255)
                                                                                                  gCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgcagca                                                                                                                        >  2:1274463/1‑54 (MQ=255)
                                                                                                  gCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTgcagca                                                                                                                        <  1:1274463/54‑1 (MQ=255)
                                                                                                   cGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTTTATAAAGCCTGCAGCGAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGTGCGCGCGTTTGACGACCGGAGCCAATtcg                                  >  1:887798/1‑139 (MQ=255)
                                                                                                                            tGCGTTCTCCTTTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAAcgccg         >  1:543468/1‑139 (MQ=255)
                                                                                                                                   tccttTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCg  >  1:748865/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GGGTTGTAGTCGTTAATGTTAACGCTTTGTCGGGCGGGAAGCGCGCGGGCGAAAGGGGTATGGAGGAAGGCAAGCGCCGAGCTTGCCGTCAGGAACGTGCGTCGGGAGAGTTTTTTAAATGGCATGCGTTCTCCTCTATAAAGCCTGCAGCAAGCTGGCGAGTTCTCGATTGATCACCTGCTGGTTAAAATCGTGTTCAACTTTTTCGCGCGCGCGTTTGACGACCGGAGCCAATTCGTCGGTGTCCAGTTGGCTAAACGCCGCCAGTCG  >  NC_000913/2115772‑2116041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: