Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,236,903 G→A 69.2% A150T (GCT→ACT)  purD → phosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,236,9030GA69.2% 150.4 / 97.6 107A150T (GCT→ACT) purDphosphoribosylglycinamide synthetase phosphoribosylamine‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base G (4/29);  new base A (42/32);  total (46/61)
Fisher's exact test for biased strand distribution p-value = 1.46e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 107 total aligned reads displayed.

AACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCG  >  minE/2236841‑2236953
                                                              |                                                  
aaCCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGc                                                >  1:936209/1‑67 (MQ=255)
aaCCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGc                                                >  1:238013/1‑67 (MQ=255)
aaCCTGCGCTGGCGTATCTGCGTGAGAAAAGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTg                                                 >  1:328364/1‑66 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:194181/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1579626/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1595536/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:836677/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:808169/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:705636/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1795508/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1803465/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1829280/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1892128/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1011036/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:642501/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:2048616/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:45833/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:2117691/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:402068/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:28715/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:388543/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:387747/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1119478/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1282660/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1279305/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1245346/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:967569/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1456173/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:106865/67‑1 (MQ=255)
                      tGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGGCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:507637/67‑1 (MQ=255)
                       gagaAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:211442/66‑1 (MQ=255)
                        agaAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:917957/65‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:364830/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1068785/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1947674/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:953572/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:63112/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:2000060/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:548252/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1005769/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:2065338/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1509467/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:2129321/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:415337/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:281640/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1608582/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1600966/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:812806/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1392398/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1402087/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1670750/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1388578/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:149342/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1350560/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1331211/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1316258/57‑1 (MQ=255)
                                            aTTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGGGGCGATGACGCTGgaa              <  1:2031262/57‑1 (MQ=255)
                                             ttAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:744752/56‑1 (MQ=255)
                                              tAAAGCTGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1664440/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          >  1:1117993/1‑43 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          >  1:1475716/1‑43 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGAc                     >  1:1961559/1‑48 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:928183/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:838248/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:436458/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:990186/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:435146/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:658526/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:775588/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:771846/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:945208/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:680346/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1757922/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1018513/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:109541/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1279573/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1377190/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1442075/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1510508/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:15651/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1634473/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1646381/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1666424/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1666476/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1681841/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:371213/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1896940/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1988061/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:1988285/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:2040155/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:2066099/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:2072987/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:2106979/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:23182/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:277188/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagcg  >  1:311848/1‑67 (MQ=255)
                                              tAAAGCGGACGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagc   >  1:235955/1‑66 (MQ=255)
                                              tAAAGCGGACGGTCAGACTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAgaagcggaagc   >  1:1894074/1‑66 (MQ=255)
                                                     gACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGgaa              <  1:1587945/48‑1 (MQ=255)
                                                      aCGGTCTGACTGCCGGGAAAGGCGTTATCGTGGCg                          <  1:1792057/35‑1 (MQ=255)
                                                              |                                                  
AACCTGCGCTGGCGTATCTGCGTGAGAAAGGCGCGCCAATCGTCATTAAAGCGGACGGTCTGGCTGCCGGGAAAGGCGTTATCGTGGCGATGACGCTGGAAGAAGCGGAAGCG  >  minE/2236841‑2236953

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: