Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,183,240 G→A 46.1% E1172E (GAG→GAA gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,183,2400GA46.1% 8.4 / 52.4 65E1172E (GAG→GAAgltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (13/22);  new base A (11/19);  total (24/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGA  >  minE/2183175‑2183304
                                                                 |                                                                
cAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGttt                                                               <  1:938810/69‑1 (MQ=255)
cAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGttt                                                               <  1:3416095/69‑1 (MQ=255)
             tGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGaaa                                                  <  1:3107955/69‑1 (MQ=255)
                      aCCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAg                                          >  1:2086867/1‑68 (MQ=255)
                      aCCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAg                                          >  1:251866/1‑68 (MQ=255)
                      aCCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAg                                          >  1:2562227/1‑68 (MQ=255)
                      aCCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAg                                          >  1:275429/1‑68 (MQ=255)
                      aCCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAg                                          >  1:1454877/1‑68 (MQ=255)
                        cACTATCACGGCCTGCCATTCAAGGTTACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCt                                        <  1:330542/68‑1 (MQ=255)
                        cACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCt                                        <  1:181779/68‑1 (MQ=255)
                        cACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCt                                        <  1:1384676/68‑1 (MQ=255)
                        cACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCt                                        <  1:2272918/68‑1 (MQ=255)
                         aCTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGa                                      >  1:1403252/1‑69 (MQ=255)
                         aCTATCACGGCCTGCCATTCAAGGTAACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGa                                      >  1:1749127/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:2964245/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:2902941/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:425398/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:831913/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:2116110/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:951117/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:1083758/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:1283896/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:1407531/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:1336597/1‑69 (MQ=255)
                          cTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGAt                                     >  1:1359752/1‑69 (MQ=255)
                                            cAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTaa                    >  1:2095017/1‑68 (MQ=255)
                                            cAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAac                   >  1:2127412/1‑69 (MQ=255)
                                            cAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAac                   >  1:858721/1‑69 (MQ=255)
                                            cAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAac                   >  1:1817550/1‑69 (MQ=255)
                                            cAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAac                   >  1:2616827/1‑69 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCTGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:2700259/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:3169364/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:1624871/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:3302912/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:2874393/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:2763931/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:2607554/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:2568197/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:1849582/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:531840/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:623867/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:2407961/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:777590/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:2272949/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:1062326/69‑1 (MQ=255)
                                                      aaTTACTTTGAATTTATCGCCCGTAAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:2114303/69‑1 (MQ=255)
                                                       aTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:3205149/68‑1 (MQ=255)
                                                       aTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:213027/68‑1 (MQ=255)
                                                       aTTACTTTGAATTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTGGGCGTAACACGTCtggtg         <  1:610816/68‑1 (MQ=255)
                                                            tttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCt    <  1:1568929/68‑1 (MQ=255)
                                                            tttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCt    <  1:2027774/68‑1 (MQ=255)
                                                            tttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCt    <  1:82702/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:1917901/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:2504694/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:1818415/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:2417586/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:695210/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:231637/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:78273/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:28719/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:2882245/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:886269/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:2125916/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTg   <  1:3114333/68‑1 (MQ=255)
                                                             ttGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGa  >  1:1962934/1‑69 (MQ=255)
                                                                 |                                                                
CAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCCGTGAAACCCGCGAGCTGATGGCACAGCTTGGCGTAACACGTCTGGTGGATCTGA  >  minE/2183175‑2183304

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: