Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 441,081 G→T 52.3% V162V (GTC→GTA potE ← putrescine/proton symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE441,0810GT52.3% ‑4.9 / 49.3 42V162V (GTC→GTApotEputrescine/proton symporter
Reads supporting (aligned to +/- strand):  ref base G (13/7);  new base T (14/8);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

ATGCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATGCG  >  minE/441015‑441137
                                                                  |                                                        
atgCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACcgac                                                        >  1:7307/1‑69 (MQ=255)
         ccAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGaat                                               <  1:1453806/69‑1 (MQ=255)
         ccAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGaat                                               <  1:1969284/69‑1 (MQ=255)
         ccAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGaat                                               >  1:3101177/1‑69 (MQ=255)
         ccAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGaat                                               <  1:1237639/69‑1 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:1097818/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:2779645/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:1380532/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:1454413/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:1499974/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:3209070/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:3001096/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:1069910/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:1940286/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:196640/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaata                                              >  1:2264645/1‑69 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaat                                               >  1:2189682/1‑68 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaat                                               >  1:942764/1‑68 (MQ=255)
          cAGGAATCAACGTAAAGTGTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGaa                                                >  1:574001/1‑67 (MQ=255)
                  aaCGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGAccc                                       >  1:3325297/1‑68 (MQ=255)
                    cGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGATCGACGACCGGAATAATGAccccc                                     <  1:1342493/68‑1 (MQ=255)
                           tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt                             <  1:3157311/68‑1 (MQ=255)
                           tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt                             <  1:2626151/68‑1 (MQ=255)
                           tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt                             <  1:2042939/68‑1 (MQ=255)
                           tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt                             <  1:2222550/68‑1 (MQ=255)
                           tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt                             <  1:2153518/68‑1 (MQ=255)
                           tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt                             <  1:1878285/68‑1 (MQ=255)
                           tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt                             <  1:1771044/68‑1 (MQ=255)
                           tgTCGGGCTAAACCAGAACCAGCCGATAATGCACAGGCCTACCACCGGAATAATGACCCCCCACACGGt                             <  1:1044461/68‑1 (MQ=255)
                                    aaaCCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACt                    <  1:2475284/69‑1 (MQ=255)
                                    aaaCCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACt                    <  1:1563941/69‑1 (MQ=255)
                                    aaaCCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACt                    <  1:1295316/69‑1 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGc             >  1:289816/1‑57 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:867342/1‑68 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:992940/1‑68 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:792781/1‑68 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:959992/1‑68 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:3331839/1‑68 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:3293069/1‑68 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:2973051/1‑68 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:2623765/1‑68 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:123235/1‑68 (MQ=255)
                                                       atGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATgcg  >  1:1057348/1‑68 (MQ=255)
                                                                  |                                                        
ATGCGGATTCCAGGAATCAACGTACAGCGTCGGGCTAAACCAGAACCAGCCAATAATGCACAGACCGACGACCGGAATAATGACCCCCCACACGGTAATGCTACTGATTTGCCCGGTAATGCG  >  minE/441015‑441137

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: