Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,953,017 C→A 69.7% V364V (GTG→GTT pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,953,0170CA69.7% 54.1 / 36.2 56V364V (GTG→GTTpepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base C (3/14);  new base A (12/27);  total (15/41)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGC  >  minE/1952958‑1953079
                                                           |                                                              
tctACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATg                                                        <  1:993092/68‑1 (MQ=255)
            aCCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCt                                           <  1:1825313/69‑1 (MQ=255)
                       ttCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCt                                <  1:2720012/69‑1 (MQ=255)
                       ttCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCt                                <  1:2819991/69‑1 (MQ=255)
                            gAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGc                           <  1:3289229/69‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:2353903/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:2723597/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:2415683/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:2455479/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:2863611/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:2902930/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:2936452/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:3130585/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:3230472/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:444971/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:469322/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:499868/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:510156/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:659020/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:768255/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:933463/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:1430992/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:110443/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:1177762/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:1227694/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:1309713/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:1347343/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:1630512/53‑1 (MQ=255)
                                       cgaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:2204999/53‑1 (MQ=255)
                                        gaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:214027/52‑1 (MQ=255)
                                        gaTCCTGACCATAAACACCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:3329761/52‑1 (MQ=255)
                                          tCCTGATCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:437537/1‑39 (MQ=37)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTcc                                            >  1:1419222/1‑38 (MQ=38)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:131156/1‑39 (MQ=38)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:2577899/1‑39 (MQ=38)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:35529/1‑39 (MQ=38)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:2330428/1‑39 (MQ=38)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:525719/1‑39 (MQ=38)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:3085831/1‑39 (MQ=38)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:667395/1‑39 (MQ=38)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:2928680/1‑39 (MQ=38)
                                          tCCTGACCATAAACACCAACGTCATGAACATCCAGTCCt                                           >  1:2062719/1‑39 (MQ=38)
                                          tCCTGACCATAAACACCAAAGTCATGAACATCCAGTCCt                                           >  1:3153358/1‑39 (MQ=25)
                                              gACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGAc         <  1:714860/69‑1 (MQ=255)
                                              gACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGAc         <  1:2450630/69‑1 (MQ=255)
                                              gACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGAc         <  1:1490562/69‑1 (MQ=255)
                                              gACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGAc         <  1:983971/69‑1 (MQ=255)
                                              gACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGAc         <  1:1031935/69‑1 (MQ=255)
                                               aCCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCa                     <  1:2715988/56‑1 (MQ=255)
                                               aCCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCa                     <  1:1348932/56‑1 (MQ=255)
                                               aCCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCa                     <  1:1291699/56‑1 (MQ=255)
                                               aCCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCa                     <  1:1720019/56‑1 (MQ=255)
                                                     aaCACCCACGTCATGGACATCCAGTCCTAGCCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGc  >  1:756845/1‑69 (MQ=255)
                                                     aaCACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGc  >  1:2305649/1‑69 (MQ=255)
                                                     aaCACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGc  >  1:295500/1‑69 (MQ=255)
                                                       cacCAACGTCATGAACATCCAGTCCTAACCAGTGGCt                                <  1:2088515/37‑1 (MQ=38)
                                                           |                                                              
TCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGC  >  minE/1952958‑1953079

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