Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,708,493 C→G 7.2% intergenic (+31/‑259) trpS → / → yhfZ tryptophanyl‑tRNA synthetase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,708,4930CG 7.2% 210.6 / 4.9 70intergenic (+31/‑259)trpS/yhfZtryptophanyl‑tRNA synthetase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (3/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

GGTGTACGAAGCGATTGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATG  >  Exported/2708426‑2708557
                                                                   |                                                                
ggTGTACGAAGCGATTGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCa                                                                 >  1:4376342/1‑69 (MQ=255)
 gtgtACGAAGCGATTGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGc                                                              <  1:1940283/71‑1 (MQ=255)
    tACGAAGCGATTGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCg                                                            >  1:1421734/1‑70 (MQ=255)
    tACGAAGCGATTGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCg                                                            <  1:3786150/70‑1 (MQ=255)
        aaGCGATTGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATgc                                                       >  1:3180937/1‑71 (MQ=255)
         aGCGATTGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATgcg                                                      >  1:1250861/1‑71 (MQ=255)
           cGATTGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATgcggc                                                    <  1:4047004/71‑1 (MQ=255)
              ttGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTg                                                 <  1:5413398/71‑1 (MQ=255)
              ttGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTg                                                 <  1:560212/71‑1 (MQ=255)
              ttGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTg                                                 <  1:1736202/71‑1 (MQ=255)
              ttGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTg                                                 <  1:4212794/71‑1 (MQ=255)
              ttGGTTTTGTGGCGAAGCCGTAATTCACTAATGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTg                                                 <  1:2319651/71‑1 (MQ=255)
               tGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTg                                                 >  1:4471182/1‑70 (MQ=255)
               tGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTg                                                 >  1:2688817/1‑70 (MQ=255)
                ggTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTc                                                                  >  1:4876917/1‑52 (MQ=255)
                    ttGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                           >  1:797493/1‑71 (MQ=255)
                     tgtgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                           <  1:5037687/70‑1 (MQ=255)
                      gtgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTg                                                 <  1:5219663/63‑1 (MQ=255)
                       tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTTCGCTTTGTTCATGCCgg                                                           <  1:363775/52‑1 (MQ=255)
                       tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCgg                                                           <  1:5222712/52‑1 (MQ=255)
                       tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCgg                                                           <  1:1200752/52‑1 (MQ=255)
                       tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                         <  1:610389/70‑1 (MQ=255)
                       tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                        <  1:1204911/71‑1 (MQ=255)
                       tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                        <  1:736673/71‑1 (MQ=255)
                       tgGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                        <  1:186081/71‑1 (MQ=255)
                       tgGAGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCgg                                                           <  1:2591235/52‑1 (MQ=255)
                          cGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGt                                                  >  1:3838837/1‑58 (MQ=255)
                          cGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                         <  1:1095803/67‑1 (MQ=255)
                            aaGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                        <  1:2277731/66‑1 (MQ=255)
                               ccGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                 >  1:4915527/1‑70 (MQ=255)
                                 gTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTg                                                 <  1:5630881/52‑1 (MQ=255)
                                 gTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                               >  1:5523793/1‑70 (MQ=255)
                                  tAATTCACTATTGCTGGCAAATTGCGCTTTGTTCa                                                                 >  1:4698067/1‑35 (MQ=255)
                                  tAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGtgtg                                                 <  1:1891457/51‑1 (MQ=255)
                                  tAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTAc                              >  1:23368/1‑70 (MQ=255)
                                    aTTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATgcggcg                                                   >  1:5307026/1‑47 (MQ=255)
                                    aTTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                               <  1:5000217/67‑1 (MQ=255)
                                     ttCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAg                          <  1:4359624/71‑1 (MQ=255)
                                      tCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                               >  1:2076264/1‑49 (MQ=255)
                                      tCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                               >  1:3076663/1‑49 (MQ=255)
                                       cACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGa                                                          >  1:1941584/1‑37 (MQ=255)
                                       cACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATgcggcg                                                   >  1:1441757/1‑44 (MQ=255)
                                           aTTGTTGGCAAAGTGTGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                       >  1:5060043/1‑52 (MQ=11)
                                           aTTGTTGGCAAAGTGTGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                       >  1:3300924/1‑52 (MQ=11)
                                           aTTGTTGGCAAAGTGTGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                       >  1:4534983/1‑52 (MQ=11)
                                            ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATgcgg                                                     >  1:3865468/1‑37 (MQ=255)
                                            ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATgcgg                                                     >  1:4583950/1‑37 (MQ=255)
                                            ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATgcgg                                                     >  1:4387071/1‑37 (MQ=255)
                                            ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATgagg                                                     >  1:4869840/1‑37 (MQ=255)
                                            ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                          >  1:5399815/1‑48 (MQ=255)
                                            ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGc                    <  1:1674458/70‑1 (MQ=255)
                                            ttGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGc                    <  1:56510/70‑1 (MQ=255)
                                            ttGCTGGCAAATTGCGCTTTGTTAATGCCGGATgcgg                                                     >  1:3586509/1‑37 (MQ=255)
                                              gCTGGCAAATTGCGCTTTGTTCATGCCGGATgcggc                                                    >  1:957150/1‑36 (MQ=255)
                                              gCTGGCAAATTGCGCTTTGTTCATGCCGGATgcggc                                                    >  1:247301/1‑36 (MQ=255)
                                              gCTGGCAAATTGCGCTTTGTTCATGCCGGATgcggc                                                    >  1:3234032/1‑36 (MQ=255)
                                              gCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGcaca                 >  1:4685801/1‑71 (MQ=255)
                                                tGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                        <  1:859954/46‑1 (MQ=255)
                                                 aacAACTTGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCgg                                   <  1:5467278/48‑1 (MQ=11)
                                                     acttGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCg                                    >  1:3127233/3‑45 (MQ=11)
                                                     acttGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCg                                    >  1:3592912/3‑45 (MQ=11)
                                                     acttGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCg                                    >  1:2113157/3‑45 (MQ=11)
                                                     aaTTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCa          <  1:1939108/71‑1 (MQ=255)
                                                      cttGCGCTTTGTTGATGCCGGATGCGGCGTGAATTCCTTATCCg                                    <  1:2178324/43‑1 (MQ=11)
                                                       ttGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAg        <  1:421498/71‑1 (MQ=255)
                                                        tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTaa                       >  1:580117/1‑55 (MQ=38)
                                                         gcgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAgg       >  1:4676054/1‑70 (MQ=255)
                                                          cgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGaa     <  1:1651142/71‑1 (MQ=255)
                                                            cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAAt   <  1:3270070/71‑1 (MQ=255)
                                                             tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATg  >  1:1982313/1‑71 (MQ=255)
                                                               tGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATg  <  1:4888478/69‑1 (MQ=255)
                                                                   |                                                                
GGTGTACGAAGCGATTGGTTTTGTGGCGAAGCCGTAATTCACTATTGCTGGCAAATTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAGTAATGCACATATTGCAAGGAAATG  >  Exported/2708426‑2708557

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