Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,756,447 G→A 44.4% V345V (GTG→GTA pssA → phosphatidylserine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,756,4470GA44.4% 53.1 / 153.9 189V345V (GTG→GTApssAphosphatidylserine synthase
Reads supporting (aligned to +/- strand):  ref base G (38/67);  new base A (28/56);  total (66/123)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

Only 100 of 197 total aligned reads displayed.

GAGATCAATCTGCGTCGTTTCCTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTTGATGATAAGTGGAT  >  minE/1756385‑1756515
                                                              |                                                                    
gagaTCAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTct                                                                 <  1:939776/68‑2 (MQ=255)
gagaTCAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:2188690/68‑1 (MQ=255)
gagaTCAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:2260874/68‑1 (MQ=255)
gagaTCAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:113575/68‑1 (MQ=255)
gagaTCAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:2368325/68‑1 (MQ=255)
gagaTCAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:2619505/68‑1 (MQ=255)
gagaTCAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:3007318/68‑1 (MQ=255)
gagaTCAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:3073973/68‑1 (MQ=255)
gagaTCAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:420768/68‑1 (MQ=255)
gagaTAAATCTGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:2275045/68‑1 (MQ=255)
         cTGCGTCGTTTCCTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGaaa                                                       <  1:600332/69‑1 (MQ=255)
         cTGCGTCGTTTCCTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGATTCGGTTATGGaaa                                                       <  1:2520135/69‑1 (MQ=255)
          tGCGTCGTTTCCTGAGTCGTTTACAGTATTACGTCAACACTGACCAGTTAGTAGTTCg                                                                 <  1:1399586/58‑1 (MQ=255)
          tGCGTCGTTTCCTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGaaa                                                       <  1:2357163/68‑1 (MQ=255)
                gtTTCCTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATgac                                                 >  1:918884/1‑68 (MQ=255)
                gtTTACTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCg                                                                 >  1:2149458/1‑52 (MQ=255)
                  ttCCTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATgacgac                                              <  1:851499/69‑1 (MQ=255)
                  ttCCTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATgacgac                                              <  1:2089991/69‑1 (MQ=255)
                  ttCCTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATgacgac                                              <  1:812227/69‑1 (MQ=255)
                       gAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGacaaca                                          <  1:1152528/68‑1 (MQ=255)
                       gAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGacaaca                                          <  1:1941929/68‑1 (MQ=255)
                       gAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGacaaca                                          <  1:1895564/68‑1 (MQ=255)
                       gAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGacaaca                                          <  1:3355795/68‑1 (MQ=255)
                       gAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGacaaca                                          <  1:1249782/68‑1 (MQ=255)
                       gAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGacaaca                                          <  1:1805468/68‑1 (MQ=255)
                            gTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATc                                    <  1:3237586/69‑1 (MQ=255)
                            gTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATc                                    <  1:2818523/69‑1 (MQ=255)
                               tacAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCt                                       <  1:64711/61‑1 (MQ=255)
                               tacAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCt                                       <  1:617887/61‑1 (MQ=255)
                               tacAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCt                                       <  1:2282048/61‑1 (MQ=255)
                               tacAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCt                                       <  1:3214553/61‑1 (MQ=255)
                               tacAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCt                                       <  1:199133/61‑1 (MQ=255)
                               tacAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCt                                       <  1:1011682/61‑1 (MQ=255)
                               tGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCAc                                  >  1:2703158/1‑68 (MQ=255)
                               tGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCAc                                  >  1:2940679/1‑68 (MQ=255)
                                acAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCt                                       <  1:430604/61‑1 (MQ=255)
                                 cAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCTATCACCt                                <  1:1042082/68‑1 (MQ=255)
                                 cAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCTATCACCt                                <  1:2328344/68‑1 (MQ=255)
                                 cAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCTATCACCt                                <  1:685084/68‑1 (MQ=255)
                                 cAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCTATCACCt                                <  1:1981704/68‑1 (MQ=255)
                                 cAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCTATCACCt                                <  1:155502/68‑1 (MQ=255)
                                 cAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCTATCACCt                                <  1:144967/68‑1 (MQ=255)
                                 cAGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCTATCACCt                                <  1:1802203/68‑1 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATgac                                                 >  1:2799500/1‑50 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:1595799/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:2916145/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:3136334/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:2751207/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:1489310/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:26503/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:177498/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:2386005/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:558198/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:2350388/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:1274363/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:1127024/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:94965/1‑55 (MQ=255)
                                  aGTATTACGTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCTATCACCt                                <  1:2926369/67‑1 (MQ=255)
                                  aGTATTACGACAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGacaa                                            >  1:2856131/1‑55 (MQ=255)
                                      ttACGTCAATACTGACCAGCTATTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGATAgg                          <  1:41267/69‑1 (MQ=255)
                                      ttACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCCCCTGAAAgg                          <  1:2209399/69‑1 (MQ=255)
                                      ttACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAgg                          >  1:452965/1‑69 (MQ=255)
                                      ttACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAgg                          <  1:2345761/69‑1 (MQ=255)
                                      ttACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAgg                          <  1:1079087/69‑1 (MQ=255)
                                       tACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAgg                          <  1:271994/68‑1 (MQ=255)
                                        aCGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGa                        <  1:2964861/69‑1 (MQ=255)
                                        aCGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGa                        <  1:932620/69‑1 (MQ=255)
                                        aCGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGa                        >  1:1568515/1‑69 (MQ=255)
                                        aCGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGa                        <  1:1359223/69‑1 (MQ=255)
                                        aCGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGa                        <  1:1965700/69‑1 (MQ=255)
                                        aCGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGa                        <  1:1851993/69‑1 (MQ=255)
                                          gTCAACACTGACCAGTTAGTAGTTCGGTTATGGAAAGATGACGACAACACCTATCACCt                                <  1:1657480/59‑1 (MQ=255)
                                            cAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGAtgtg                    <  1:2237799/69‑1 (MQ=255)
                                            cAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGAtgtg                    <  1:2746133/69‑1 (MQ=255)
                                            cAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGAtgtg                    <  1:3307606/69‑1 (MQ=255)
                                            cAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGAtgtg                    <  1:1184458/69‑1 (MQ=255)
                                            cAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGAtgtg                    <  1:767440/69‑1 (MQ=255)
                                            cAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGAtgtg                    <  1:1190439/69‑1 (MQ=255)
                                             aaTACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTgg                   <  1:1066702/69‑1 (MQ=255)
                                             aaTACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTgg                   <  1:1714798/69‑1 (MQ=255)
                                             aaTACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTgg                   <  1:1735889/69‑1 (MQ=255)
                                             aaTACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTgg                   <  1:385948/69‑1 (MQ=255)
                                             aaTACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTgg                   <  1:1251077/69‑1 (MQ=255)
                                             aaTACTGACCAGCTAGTGGTCCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTgg                   <  1:203566/69‑1 (MQ=255)
                                                aCTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGAt                       >  1:2054592/1‑62 (MQ=255)
                                                aCTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGtt                >  1:121998/1‑69 (MQ=255)
                                                aCTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGtt                >  1:2597088/1‑69 (MQ=255)
                                                   gACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTgg                   <  1:549653/63‑1 (MQ=255)
                                                    aCCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTtgat             >  1:425286/1‑68 (MQ=255)
                                                    aCCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTtgat             >  1:631186/1‑68 (MQ=255)
                                                    aCCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTtgat             >  1:2576016/1‑68 (MQ=255)
                                                    aCCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTtgat             >  1:1748140/1‑68 (MQ=255)
                                                    aCCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTtgat             >  1:3169465/1‑68 (MQ=255)
                                                    aCCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTtgat             >  1:2114507/1‑68 (MQ=255)
                                                       aGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTTGATGATa         <  1:1925737/69‑1 (MQ=255)
                                                           aGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTTGATGATAAGTg     <  1:3345179/69‑1 (MQ=255)
                                                              ggTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTTGATGATa         >  1:248907/1‑62 (MQ=255)
                                                              ggTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTTGATGATAAGTGGAt  >  1:1810127/1‑69 (MQ=255)
                                                              ggTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTTGATGATAAGTGGAt  >  1:266097/1‑69 (MQ=255)
                                                              ggTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTTGATGATAAGTGGAt  >  1:3032896/1‑69 (MQ=255)
                                                              |                                                                    
GAGATCAATCTGCGTCGTTTCCTGAGCCGTTTGCAGTATTACGTCAATACTGACCAGCTAGTGGTTCGGTTATGGAAAGATGACGACAACACCTATCACCTGAAAGGGATGTGGGTTGATGATAAGTGGAT  >  minE/1756385‑1756515

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