Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,317 G→T 32.0% P345P (CCG→CCT narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,3170GT32.0% 91.1 / 34.2 53P345P (CCG→CCTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (20/16);  new base T (12/5);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

AAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGACGCG  >  minE/836250‑836386
                                                                   |                                                                     
aaGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGaa                                                                    >  1:45822/1‑71 (MQ=255)
   ccGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCAGGAatat                                                                 <  1:2303493/71‑1 (MQ=255)
     gCACAGCAGTCGCCGGTTTGTAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCg                                                               <  1:3517377/71‑1 (MQ=255)
       acaGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGca                                                             >  1:3843525/1‑71 (MQ=255)
       acaGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGca                                                             >  1:3552912/1‑71 (MQ=255)
       acaGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGca                                                             >  1:2114912/1‑71 (MQ=255)
              tCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGcc                                                      <  1:177900/71‑1 (MQ=255)
                gCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCg                                                     <  1:1710379/70‑1 (MQ=255)
                   ggTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCTGAATATCGCacac                                                          >  1:1663398/1‑62 (MQ=255)
                     tttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTc                                               <  1:2415971/71‑1 (MQ=255)
                       tataAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTg                                             >  1:1321571/1‑71 (MQ=255)
                       tacaaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTg                                             >  1:1724754/4‑71 (MQ=255)
                         taAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGt                                                <  1:3414898/66‑1 (MQ=255)
                         caaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGaa                                                                    >  1:2648507/2‑46 (MQ=12)
                         caaaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCt                                              <  1:3323974/67‑1 (MQ=255)
                          aaaaTGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCt                                              <  1:760994/67‑1 (MQ=255)
                           aaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGaa                                                                    >  1:768524/1‑44 (MQ=12)
                           aaaTGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGaa                                                                    >  1:648273/1‑44 (MQ=12)
                             aTGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGt                                       <  1:3284045/71‑1 (MQ=255)
                                  gatggactggaAACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTtcct                                   <  1:2973165/69‑2 (MQ=255)
                                  aatggaatggAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGAttg                                                 <  1:1455029/56‑3 (MQ=255)
                                   atggaatggaAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGcc                                                      >  1:294193/1‑50 (MQ=255)
                                    tggaatggaaACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGcc                                                      >  1:3773178/1‑49 (MQ=255)
                                    tggaatggaaACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCct                                <  1:1179151/71‑1 (MQ=255)
                                    tggaatggaaACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCct                                >  1:2348671/1‑71 (MQ=255)
                                     ggactggaaACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATgg                                                 <  1:124990/53‑1 (MQ=16)
                                      gaatggaaACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCctct                              >  1:1978771/1‑71 (MQ=255)
                                       aatggaaACTGGCGCTGCCGCTGCATCCGGAATATc                                                                >  1:929345/1‑36 (MQ=255)
                                       aatggaaACTGGCGCTGCCGCTGCATCCGGAATATc                                                                >  1:1488579/1‑36 (MQ=255)
                                       aatggaaACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTg                                             >  1:702912/1‑55 (MQ=255)
                                           gaaACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGtcac                         >  1:1016794/1‑69 (MQ=255)
                                           gaaACTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGtcac                         >  1:2646398/1‑69 (MQ=255)
                                            aaaCTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGtcac                         <  1:351759/70‑3 (MQ=255)
                                            aaaCTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGTCAcc                        <  1:1189360/71‑1 (MQ=255)
                                             aaCTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTcc                        >  1:2183290/1‑70 (MQ=255)
                                             aaCTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCg                       <  1:468770/71‑1 (MQ=255)
                                             aaCTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTTCGTTCCTCCGCTGTCAcc                        >  1:1117825/1‑70 (MQ=255)
                                             aaCTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGTCAcc                        >  1:2949734/1‑70 (MQ=255)
                                             aaCTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGTCAcc                        >  1:193198/1‑70 (MQ=255)
                                               cTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTg                                             >  1:605674/1‑47 (MQ=255)
                                               cTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTg                                             >  1:1158093/1‑47 (MQ=255)
                                               cTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGAt                     >  1:1096156/1‑71 (MQ=255)
                                               cTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGTCACCGAt                     >  1:2563832/1‑71 (MQ=255)
                                                  gctctgccgctgcATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCa                  <  1:1037688/68‑1 (MQ=255)
                                                  gcgcTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCa                  <  1:938035/71‑1 (MQ=255)
                                                     ctgccgctgcATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTc               <  1:2020396/71‑1 (MQ=255)
                                                      taccgttgcACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGTCACCGATTCAGTcc              >  1:2225747/3‑70 (MQ=255)
                                                        ccgctgcATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGAtt                    >  1:3817749/1‑63 (MQ=255)
                                                              cATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTgcagca        <  1:2173391/69‑1 (MQ=255)
                                                                 ccGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGAcgc   >  1:3561121/1‑71 (MQ=255)
                                                                 ccGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGAcgc   >  1:2260819/1‑71 (MQ=255)
                                                                  cGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGAcgcg  >  1:1220806/1‑71 (MQ=255)
                                                                  cGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGAcgc   <  1:1796500/70‑1 (MQ=255)
                                                                   |                                                                     
AAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGACGCG  >  minE/836250‑836386

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: