Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,188,466 A→C 44.8% G131G (GGT→GGG nanA ← N‑acetylneuraminate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,188,4660AC44.8% ‑1.5 / 117.8 103G131G (GGT→GGGnanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base A (23/34);  new base C (13/33);  total (36/67)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-02

Only 100 of 106 total aligned reads displayed.

TTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAG  >  minE/2188401‑2188532
                                                                 |                                                                  
ttGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAcc                                                                  <  1:2277797/68‑1 (MQ=255)
          ccAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                                        >  1:485107/1‑68 (MQ=255)
          ccAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                                        >  1:759570/1‑68 (MQ=255)
          ccAGGGTCAGTTTTACCCCACTCAGGGATGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                                        >  1:1241352/1‑68 (MQ=255)
               gTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                                  >  1:3207905/1‑69 (MQ=255)
               gTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                                  >  1:2130791/1‑69 (MQ=255)
               gTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                                  >  1:3348664/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGccc                                           >  1:2246997/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGccc                                           >  1:1462338/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGccc                                           >  1:3026298/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGccc                                           >  1:1864123/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGccc                                           >  1:515358/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGccc                                           >  1:524341/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGccc                                           >  1:698378/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGccc                                           >  1:1092132/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATACGCCGAATCAATAATTGccc                                           >  1:53421/1‑69 (MQ=255)
                      ttACCCCGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATACGCCGAATCAATAATTGccc                                           >  1:1442237/1‑69 (MQ=255)
                      ttACCCAGCTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGccc                                           >  1:3467/1‑69 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:993123/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:2457971/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:1778001/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:2106412/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:2794730/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:1727271/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:1571536/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:3227598/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:1407332/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:550657/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:764262/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:863671/68‑1 (MQ=255)
                          cccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:977157/68‑1 (MQ=255)
                           ccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                                        <  1:164996/67‑1 (MQ=255)
                                cAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGAt                                  >  1:3373712/1‑68 (MQ=255)
                                cAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGAt                                  >  1:239471/1‑68 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:3202843/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:2999319/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:3389771/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:2930626/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:824288/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:1036347/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:1051878/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:1070882/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:1174105/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:135499/62‑1 (MQ=255)
                                    gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:1419431/62‑1 (MQ=255)
                                     cTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                                    <  1:2430030/61‑1 (MQ=255)
                                     cTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt                            >  1:959374/1‑69 (MQ=255)
                                     cTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt                            >  1:2875188/1‑69 (MQ=255)
                                     cTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt                            >  1:2506399/1‑69 (MQ=255)
                                      tGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtg                           <  1:2735774/69‑1 (MQ=255)
                                      tGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtg                           <  1:2982448/69‑1 (MQ=255)
                                       ggAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                                                   >  1:807515/1‑44 (MQ=255)
                                       ggAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                                                   >  1:251797/1‑44 (MQ=255)
                                       ggAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                                                   >  1:731071/1‑44 (MQ=255)
                                       ggAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                                                   >  1:686740/1‑44 (MQ=255)
                                               gtACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc                               >  1:106665/1‑56 (MQ=255)
                                               gtACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc                               >  1:1638568/1‑56 (MQ=255)
                                               gtACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc                               >  1:315074/1‑56 (MQ=255)
                                               gtACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc                               >  1:305040/1‑56 (MQ=255)
                                               gtACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCaa                   <  1:1470787/68‑1 (MQ=255)
                                                     caccaTCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGat            <  1:569110/69‑2 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:1218913/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:13022/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:424502/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2463029/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:423784/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:1158154/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:653706/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:67243/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:1334725/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:744736/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:1153779/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:1145309/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2117226/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:1025464/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2697104/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2985054/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:1803718/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:1760007/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:3057669/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:3066472/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:374298/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:3205166/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2946734/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:3337915/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2931248/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2841355/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:3355748/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2113810/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:3390538/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2544724/69‑1 (MQ=255)
                                                          tCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAACCTGAAAGGAt       <  1:604466/69‑1 (MQ=255)
                                                          tCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:677499/69‑1 (MQ=255)
                                                           cGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       <  1:2358475/68‑1 (MQ=255)
                                                           cGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtg                           <  1:36086/48‑1 (MQ=255)
                                                            ggCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       >  1:1983668/1‑67 (MQ=255)
                                                            ggCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAt       >  1:1214085/1‑67 (MQ=255)
                                                               aaaCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAtagtag  <  1:1468880/69‑1 (MQ=255)
                                                               aaaCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAtagtag  <  1:1452565/69‑1 (MQ=255)
                                                               aaaCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGAtagtag  <  1:1271043/69‑1 (MQ=255)
                                                                 |                                                                  
TTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAG  >  minE/2188401‑2188532

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