Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 441,066 A→G 70.7% I167I (ATT→ATC potE ← putrescine/proton symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE441,0660AG70.7% 48.5 / 20.0 51I167I (ATT→ATCpotEputrescine/proton symporter
Reads supporting (aligned to +/- strand):  ref base A (8/7);  new base G (28/8);  total (36/15)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

ACGGTGCATGATGCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACT  >  minE/441005‑441119
                                                             |                                                     
aCGGTGCATGATGCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATg                                                 >  1:1999833/1‑68 (MQ=255)
aCGGTGCATGATGCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATg                                                 >  1:2772264/1‑68 (MQ=255)
aCGGTGCATGATGCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATg                                                 >  1:162058/1‑68 (MQ=255)
aCGGTGCATGATGCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATg                                                 >  1:2827632/1‑68 (MQ=255)
aCGGTGCATGATGCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATg                                                 >  1:3006754/1‑68 (MQ=255)
          atgCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACcgac                                      >  1:7307/1‑69 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:535071/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:2478896/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:2498096/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:2754060/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:1846200/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:1123511/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:1192740/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:3166485/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:3237410/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:3340269/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:3393584/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:1252801/1‑60 (MQ=255)
                aTTCCAGGAATCAACGTAAAGTGTAGGGCTAAACCAGAACCAGCCGATAATGCACAGGcct                                        >  1:2800480/1‑60 (MQ=255)
                   ccAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGaat                             <  1:1453806/69‑1 (MQ=255)
                   ccAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGaat                             >  1:3101177/1‑69 (MQ=255)
                   ccAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGaat                             <  1:1237639/69‑1 (MQ=255)
                   ccAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGaat                             <  1:1969284/69‑1 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATg                                                 >  1:606375/1‑48 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:3001096/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:3209070/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:2779645/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:2264645/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:1069910/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:1499974/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:1454413/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:1380532/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:1940286/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:196640/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                            >  1:1097818/1‑69 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaat                             >  1:942764/1‑68 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaat                             >  1:2189682/1‑68 (MQ=255)
                    cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaa                              >  1:574001/1‑67 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:518655/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:869728/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:49383/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:1340280/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:1375092/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:1399934/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:902898/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:1487890/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:2895768/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:1559692/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:2185519/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:175374/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:1847588/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:2735355/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:2458604/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataat                                                  >  1:2470965/1‑44 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgataa                                                   >  1:3230618/1‑43 (MQ=255)
                       gAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCgat                                                     >  1:3069794/1‑41 (MQ=255)
                            aaCGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGAccc                     >  1:3325297/1‑68 (MQ=255)
                              cGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGATCGACGACCGGAATAATGAccccc                   <  1:1342493/68‑1 (MQ=255)
                                     tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt           <  1:2153518/68‑1 (MQ=255)
                                     tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt           <  1:2222550/68‑1 (MQ=255)
                                     tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt           <  1:1044461/68‑1 (MQ=255)
                                     tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt           <  1:1878285/68‑1 (MQ=255)
                                     tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt           <  1:2626151/68‑1 (MQ=255)
                                     tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt           <  1:1771044/68‑1 (MQ=255)
                                     tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt           <  1:3157311/68‑1 (MQ=255)
                                     tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt           <  1:2042939/68‑1 (MQ=255)
                                              aaaCCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACt  <  1:2475284/69‑1 (MQ=255)
                                              aaaCCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACt  <  1:1295316/69‑1 (MQ=255)
                                              aaaCCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACt  <  1:1563941/69‑1 (MQ=255)
                                                             |                                                     
ACGGTGCATGATGCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACT  >  minE/441005‑441119

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: