Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,020,261 A→G 30.9% S94S (AGT→AGC uidC ← predicted outer membrane porin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,020,2610AG30.9% 48.7 / 26.0 58S94S (AGT→AGCuidCpredicted outer membrane porin protein
Reads supporting (aligned to +/- strand):  ref base A (25/15);  new base G (9/9);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTTTTATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCAATAAT  >  minE/1020195‑1020323
                                                                  |                                                              
gggTTTTATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACt                                                              >  1:985651/1‑69 (MQ=255)
gggTTTTATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACt                                                              >  1:1160408/1‑69 (MQ=255)
gggTTTTATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACt                                                              >  1:2033096/1‑69 (MQ=255)
gggTTTTATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACt                                                              >  1:16896/1‑69 (MQ=255)
gggTTTTATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACt                                                              >  1:1365413/1‑69 (MQ=255)
gggTTTTATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACt                                                              >  1:297705/1‑69 (MQ=255)
   ttttATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCat                                                           <  1:2130253/69‑1 (MQ=255)
      tATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATc                                                        >  1:3183658/1‑69 (MQ=255)
      tATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATc                                                        >  1:2201451/1‑69 (MQ=255)
      tATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATc                                                        >  1:1614621/1‑69 (MQ=255)
      tATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATc                                                        >  1:738760/1‑69 (MQ=255)
          tttACTGCCCCCAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCatc                                                     >  1:82960/1‑66 (MQ=255)
          tttACTGCCCCCAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAGc                                                     >  1:1432956/1‑68 (MQ=255)
          tttACTGCCCCCAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAGc                                                     >  1:2380314/1‑68 (MQ=255)
          tttACTGCCCCCAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAGc                                                     >  1:2588953/1‑68 (MQ=255)
          tttACTGCCCCCAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAGc                                                     >  1:1573133/1‑68 (MQ=255)
          tttACTGCCCCCAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAGc                                                     >  1:2116953/1‑68 (MQ=255)
          tttACTGCCCCCAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAGc                                                     >  1:3240520/1‑68 (MQ=255)
          tttACTGCCCCCAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAGc                                                     >  1:1934397/1‑68 (MQ=255)
          tttACTGCCCCCAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAGc                                                     >  1:2861695/1‑68 (MQ=255)
                    ccAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAg                                                      <  1:2256537/57‑1 (MQ=255)
                    ccAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAg                                                      <  1:1154719/57‑1 (MQ=255)
                    ccAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAg                                                      <  1:2608729/57‑1 (MQ=255)
                    ccAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAg                                                      <  1:2841783/57‑1 (MQ=255)
                    ccAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAg                                                      <  1:1892261/57‑1 (MQ=255)
                    ccAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAg                                                      <  1:1565768/57‑1 (MQ=255)
                    ccAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAg                                                      <  1:1522804/57‑1 (MQ=255)
                    ccAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAg                                                      <  1:318584/57‑1 (MQ=255)
                    ccAGGGTAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTGCTCATATCAg                                                      <  1:1453888/57‑1 (MQ=255)
                         gCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTta                                     <  1:726813/69‑1 (MQ=255)
                         gCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTta                                     <  1:1125977/69‑1 (MQ=255)
                         gCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTta                                     <  1:869918/69‑1 (MQ=255)
                         gCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTta                                     <  1:1382282/69‑1 (MQ=255)
                          cAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTta                                     >  1:2193722/1‑68 (MQ=255)
                          cAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTta                                     >  1:2125344/1‑68 (MQ=255)
                          cAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTta                                     >  1:1305486/1‑68 (MQ=255)
                          cAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTta                                     >  1:1370217/1‑68 (MQ=255)
                          cAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTta                                     >  1:1948810/1‑68 (MQ=255)
                                cAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACat                              <  1:2201732/69‑1 (MQ=255)
                                cAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACat                              <  1:1053518/69‑1 (MQ=255)
                                cAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACat                              <  1:887314/69‑1 (MQ=255)
                                cAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACat                              <  1:534259/69‑1 (MQ=255)
                                cAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACat                              <  1:1604207/69‑1 (MQ=255)
                                cAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACat                              <  1:2493959/69‑1 (MQ=255)
                                                aTCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTc              <  1:1274752/69‑1 (MQ=255)
                                                aTCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAAAATAATACGCCCCGCCTTc              <  1:1974320/69‑1 (MQ=255)
                                                   aaGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAAc           >  1:123821/1‑69 (MQ=255)
                                                        accaccGGGTACTCATATCAGCGCGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGcc        >  1:1607277/1‑67 (MQ=255)
                                                        accaccGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGcc        >  1:303955/1‑67 (MQ=255)
                                                        accaccGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGcc        >  1:1504400/1‑67 (MQ=255)
                                                        accaccGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCa       >  1:2932684/1‑68 (MQ=255)
                                                        accaccGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCa       >  1:2633659/1‑68 (MQ=255)
                                                        accaccGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCa       >  1:88751/1‑68 (MQ=255)
                                                           accGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCaataa   >  1:2773831/1‑69 (MQ=255)
                                                           accGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCaataa   >  1:208491/1‑69 (MQ=255)
                                                           accGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCaataa   >  1:841065/1‑69 (MQ=255)
                                                            ccGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCaataat  <  1:3166426/69‑1 (MQ=255)
                                                            ccGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCaataat  <  1:2213568/69‑1 (MQ=255)
                                                                  |                                                              
GGGTTTTATTTTTACTGCCCCCAGGGCAAAGCCAAAACTTTTATCACCATCAAGATACCACCGGGTACTCATATCAGCACGAGCACCTAATTTATAAACATAATACGCCCCGCCTTCAACGCCAATAAT  >  minE/1020195‑1020323

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: