Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 958,771 A→T 10.2% intergenic (+30/‑41) serC → / → aroA phosphoserine/phosphohydroxythreonine aminotransferase/3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913958,7710AT10.2% 95.9 / 2.4 39intergenic (+30/‑41)serC/aroAphosphoserine/phosphohydroxythreonine aminotransferase/3‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (15/19);  new base T (2/2);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

TGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGC  >  NC_000913/958641‑958902
                                                                                                                                  |                                                                                                                                   
tgGTCGGCGGAATGCGCGCTTCTATTTATCACGCCATGCCGCTGGAAGTCGTTACAGCGCTGACTGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGtt                                                                                                                         <  2:721168/143‑1 (MQ=255)
   tCGGTGGAATGCGCGCTTCTAGTTATAACGCTATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCGGGGTTGTGTTCGGACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGGCTGGCTGTGGGGattaa                                                                                                                      <  2:236218/143‑3 (MQ=255)
     ggcggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAt                                                                                                                    <  2:203516/143‑1 (MQ=255)
     ggcggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGATTTTAt                                                                                                                    >  1:339275/1‑143 (MQ=255)
      gcggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAtt                                                                                                                   >  1:513263/1‑143 (MQ=255)
      gcggAATGCGCGCTTCTATTTATAACGCCATGCAGCTGGAAGGCGTTAAAGCGCTGACAGACGTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAAGTGTGCTTAGTCCCCCCAGCCAGCCTGTGGGGTGTCTAtt                                                                                                                   >  1:100728/1‑143 (MQ=255)
       cggAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTAttt                                                                                                                  >  1:392329/1‑143 (MQ=255)
         gcgtGCGCGCTCCTATTTATAGCGCCAGGACGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCt                                                                                                                <  1:3708/140‑1 (MQ=255)
            tGCGCGCTTCTATTTGTGACGCCTTGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCtgtt                                                                                                             <  2:748894/143‑1 (MQ=255)
                   ttCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTagaga                                                                                                      <  1:575234/143‑1 (MQ=255)
                        tttATTGCGTCATGCCGCTGGAAGGCGATAGAGCGCTGACAGGCTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAgagttga                                                                                                 <  2:715098/143‑1 (MQ=255)
                              aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                           >  1:192896/1‑143 (MQ=255)
                              aCGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCAt                                                                                           <  1:427499/143‑1 (MQ=255)
                               cGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATg                                                                                          <  1:739501/143‑1 (MQ=255)
                                  cATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGGTTTTATTTCTGTTGTAGAGGGTTGAGTTCATGGaa                                                                                       >  2:157705/1‑143 (MQ=255)
                                     gCCGCTGGAAGGCGTTAAAGCGCTGACAGACTGCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAACTTTTGCTTAATCCCCACAGCCAGCCTGTGGGGACTTTATTTCTGTTGTAGAGGGTTGAGCTCATGTAATcc                                                                                    >  2:466995/1‑143 (MQ=255)
                                      ccGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATccc                                                                                   <  1:2024/143‑1 (MQ=255)
                                       cGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCt                                                                                  <  1:666529/143‑1 (MQ=255)
                                               ggCGTTCAAGCGCTGGCAGACTGCATCGTTGAGTTCGAACGCCGTCACGGTGAATGCCGAAATTTTGCTTAATCCGCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTAc                                                                          <  2:831960/143‑1 (MQ=255)
                                                 cgctgaaGCGGTGACAGGCTGGATGGTTGAGTTCGCACGCCGTCACGGTTAATGCCGAAATGTTGCGTAAGCCCCACAGGCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACaa                                                                        <  2:50938/138‑1 (MQ=255)
                                                  gTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATATTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAAc                                                                       >  2:702515/1‑143 (MQ=255)
                                                  gTTAAAGCGCTGACAGACTTCATGGTTAAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGCTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGATACAAc                                                                       >  2:49525/1‑143 (MQ=255)
                                                     aaaCCGCTGACAGACCTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCa                                                                    <  1:26348/143‑1 (MQ=255)
                                                          gcTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCt                                                               <  1:627149/143‑1 (MQ=255)
                                                                  aCTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGCGTTCATGGAATCCCTGACGTTACAACCCACCGCCCGACTCGa                                                       >  2:764959/1‑143 (MQ=255)
                                                                         ggTTGAGTTCGAACGCCGTCACGGTGAATGCCGAACTTTCGCTTAATCCCCACAGCCGGCCTGTGGGGTTTTTATTTCTGGTGTCGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACt                                                <  2:647036/143‑1 (MQ=255)
                                                                          gTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCAACAGCATGTACGTGTGTGGTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTa                                               <  1:203906/143‑1 (MQ=255)
                                                                                ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATc                                         >  2:258151/1‑143 (MQ=255)
                                                                                ttCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATc                                         >  1:352981/1‑143 (MQ=255)
                                                                                    aaCGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGcc                                     <  1:572850/143‑1 (MQ=255)
                                                                                     aCGCGGTCGCGGTTAATGCCGAAATTTTGCTTAATCCCCGCAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGccc                                    <  2:248350/143‑1 (MQ=255)
                                                                                        ccctgaCGGTTACTGCCGAACTTTTGCTTGATCCCCACGGCCACCCTGTGGGGTTTTTCTTTCTGTTGTAGGGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGt                                 <  2:71111/138‑1 (MQ=255)
                                                                                            cACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACATCACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCa                             >  2:98820/1‑143 (MQ=255)
                                                                                             gcGGTTAATGCCGAAATTTTGCTTAATCCCCGCAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCaa                            <  2:846290/142‑1 (MQ=255)
                                                                                               ggTTAATGCTGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTCGCGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAaga                          <  2:145659/143‑1 (MQ=255)
                                                                                                            aTTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcg             >  1:699056/1‑143 (MQ=255)
                                                                                                            aTTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACcg             >  1:571488/1‑143 (MQ=255)
                                                                                                                      aaccccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTGTATTg   >  2:516547/3‑143 (MQ=255)
                                                                                                                       accccACAGGCTGGCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGc  >  1:869711/2‑143 (MQ=255)
                                                                                                                                  |                                                                                                                                   
TGGTCGGCGGAATGCGCGCTTCTATTTATAACGCCATGCCGCTGGAAGGCGTTAAAGCGCTGACAGACTTCATGGTTGAGTTCGAACGCCGTCACGGTTAATGCCGAAATTTTGCTTAATCCCCACAGCCAGCCTGTGGGGTTTTTATTTCTGTTGTAGAGAGTTGAGTTCATGGAATCCCTGACGTTACAACCCATCGCTCGTGTCGATGGCACTATTAATCTGCCCGGTTCCAAGAGCGTTTCTAACCGCGCTTTATTGC  >  NC_000913/958641‑958902

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