Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 27,559 A→C 36.5% A77A (GCA→GCC carA → carbamoyl phosphate synthetase small subunit, glutamine amidotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE27,5590AC36.5% 73.8 / 111.7 126A77A (GCA→GCCcarAcarbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base A (50/30);  new base C (21/25);  total (71/55)
Fisher's exact test for biased strand distribution p-value = 9.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

Only 100 of 147 total aligned reads displayed.

ACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCT  >  minE/27494‑27621
                                                                 |                                                              
acttaTCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGcaca                                                              >  1:3356755/1‑68 (MQ=255)
acttaTCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACaa                                                             >  1:1024426/1‑69 (MQ=255)
         cATATTGGCAATTTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:1455524/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:810251/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:1621229/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:88791/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:264938/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:2937648/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:2607729/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:245105/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:1212293/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:2090354/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:2061845/68‑1 (MQ=255)
         cATATTGGCAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:413961/68‑1 (MQ=255)
         cATATTGGAAATGTCGGCACTAATGACGCCGATGAAGAATCTTCTCAGGTTCATGCCCAAGGTCTGGt                                                     <  1:2629200/68‑1 (MQ=255)
             ttGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGAtta                                                >  1:1291641/1‑68 (MQ=255)
             ttGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGAt                                                  >  1:1534009/1‑67 (MQ=255)
             ttGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTc                                                >  1:2681245/1‑69 (MQ=255)
             ttGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTc                                                >  1:2881136/1‑69 (MQ=255)
             ttGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTc                                                >  1:1438525/1‑69 (MQ=255)
             ttGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTc                                                >  1:889666/1‑69 (MQ=255)
             ttGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTc                                                >  1:1000881/1‑69 (MQ=255)
             ttGGCAATGTCGGCACCAATGACGCCGATGAACAATCTTCTCAGGTACATGCACAAGGTCTTGTGATTc                                                >  1:662870/1‑69 (MQ=255)
               ggCATTGTCGGCACTAATGACGCCGATGAAGAATCTTCACAGGTTCATGCCCAAGGTCTGGt                                                     <  1:782233/62‑1 (MQ=255)
                 cAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgcg                                             <  1:424802/68‑1 (MQ=255)
                 cAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgcg                                             <  1:1190458/68‑1 (MQ=255)
                 cAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgcg                                             <  1:198714/68‑1 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGcac                                                               >  1:2333408/1‑45 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGcac                                                               >  1:480373/1‑45 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGcac                                                               >  1:3112588/1‑45 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGcac                                                               >  1:740235/1‑45 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGcac                                                               >  1:1323986/1‑45 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGc                                       >  1:1883140/1‑69 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGc                                       >  1:2104704/1‑69 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGc                                       >  1:788851/1‑69 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGc                                       >  1:1237738/1‑69 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGc                                       >  1:1367647/1‑69 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGc                                       >  1:2927853/1‑69 (MQ=255)
                      tCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGc                                       >  1:1566534/1‑69 (MQ=255)
                      tCGGCACCAAAGACGCCGATGAAGAATCTTCTCAGGTACATGcac                                                               >  1:1850671/1‑45 (MQ=255)
                              aaTGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGAtt                               >  1:406178/1‑69 (MQ=255)
                              aaTGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGAtt                               >  1:3164024/1‑69 (MQ=255)
                              aaTGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGAtt                               >  1:1083177/1‑69 (MQ=255)
                              aaTGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGAtt                               >  1:1509567/1‑69 (MQ=255)
                                          gaagaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgc                                              <  1:1793335/42‑1 (MQ=255)
                                          gaagaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgc                                              <  1:199105/42‑1 (MQ=255)
                                          gaagaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgc                                              <  1:3247754/42‑1 (MQ=255)
                                          gaagaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgc                                              <  1:2713489/42‑1 (MQ=255)
                                          gaagaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgc                                              <  1:2256993/42‑1 (MQ=255)
                                          gaagaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgc                                              <  1:2291299/42‑1 (MQ=255)
                                          gaagaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTcgc                                              <  1:2999797/42‑1 (MQ=255)
                                            agaaTCTTCTCAGGTTCATGCCCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTcc                  <  1:671613/68‑1 (MQ=255)
                                            agaaTCTTCTCAGGTTCATGCCCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTcc                  <  1:3293338/68‑1 (MQ=255)
                                            agaaTCTTCTCAGGTTCATGCCCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTcc                  <  1:660370/68‑1 (MQ=255)
                                            agaaTCTTCTCAGGTTCATGCCCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTcc                  <  1:1703999/68‑1 (MQ=255)
                                            agaaTCTTCTCAGGTTCATGCCCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTcc                  <  1:1691830/68‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:2845930/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:971288/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:2679654/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:3195027/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:850958/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:2735713/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:524675/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:536920/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:2917671/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:1804755/69‑1 (MQ=255)
                                             gaaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGACAGCAACTTCCGt                <  1:513675/69‑1 (MQ=255)
                                              aaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:1373749/68‑1 (MQ=255)
                                              aaTCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGt                <  1:751258/68‑1 (MQ=255)
                                                tcttctCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAAt             >  1:85359/1‑69 (MQ=255)
                                                tcttctCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAAt             >  1:310183/1‑69 (MQ=255)
                                                  ttctCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATa            >  1:2731752/1‑68 (MQ=255)
                                                      cAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCg         <  1:2625123/67‑1 (MQ=255)
                                                       aGGTTCATGCCCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:1019236/1‑69 (MQ=255)
                                                       aGGTTCATGCCCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:3050468/1‑69 (MQ=255)
                                                       aGGTTCATGCCCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:164051/1‑69 (MQ=255)
                                                       aGGTTCATGCCCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:2049620/1‑69 (MQ=255)
                                                       aGGTTCATGACCAAGGTCTGGTGATTCGCGACCTGc                                       >  1:966532/1‑36 (MQ=25)
                                                       aGGTTCATGACCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:3381901/1‑69 (MQ=255)
                                                       aGGTTCATGACCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:3161882/1‑69 (MQ=255)
                                                       aGGTTCATGACCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:2990667/1‑69 (MQ=255)
                                                       aGGTTCATGACCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:937106/1‑69 (MQ=255)
                                                       aGGTTCATGACCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:272553/1‑69 (MQ=255)
                                                       aGGTTCATGACCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:2599810/1‑69 (MQ=255)
                                                       aGGTTCATGACCAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAg      >  1:3422321/1‑69 (MQ=255)
                                                       aGGTTCATGACCAAGATCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGa        >  1:2323676/1‑67 (MQ=255)
                                                           aCATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGAcc   >  1:813066/1‑68 (MQ=255)
                                                           aCATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGAcc   >  1:2738685/1‑68 (MQ=255)
                                                           aCATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGAc    >  1:2456634/1‑67 (MQ=255)
                                                           aCATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACct  >  1:2364396/1‑69 (MQ=255)
                                                           aCATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACct  >  1:610367/1‑69 (MQ=255)
                                                           aCATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACct  >  1:1853110/1‑69 (MQ=255)
                                                           aCATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACct  >  1:2932274/1‑69 (MQ=255)
                                                           aCATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACct  >  1:298433/1‑69 (MQ=255)
                                                                cacaAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGaa       <  1:263902/59‑1 (MQ=255)
                                                                cacaAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGaa       <  1:1872894/59‑1 (MQ=255)
                                                                cacaAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGaa       <  1:1174215/59‑1 (MQ=255)
                                                                cacaAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGaa       <  1:922883/59‑1 (MQ=255)
                                                                cacaAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGaa       <  1:3140959/59‑1 (MQ=255)
                                                                cacaAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGaa       <  1:1081495/59‑1 (MQ=255)
                                                                 |                                                              
ACTTATCCCCATATTGGCAATGTCGGCACCAATGACGCCGATGAAGAATCTTCTCAGGTACATGCACAAGGTCTGGTGATTCGCGACCTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCT  >  minE/27494‑27621

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: