Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,901 T→C 28.1% L787L (TTG→CTG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,9010TC28.1% 84.2 / 39.1 57L787L (TTG→CTG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base C (9/7);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

AGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAC  >  minE/833831‑833966
                                                                      |                                                                 
aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt                                                                   <  1:2738845/71‑1 (MQ=255)
aGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAttt                                                                   <  1:2103152/71‑1 (MQ=255)
 ggCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCAtttt                                                                  <  1:3019309/71‑1 (MQ=255)
  gCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTg                                                                 >  1:465022/1‑71 (MQ=255)
      gCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGa                                                             >  1:996399/1‑71 (MQ=255)
       cTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGAc                                                            <  1:2875369/71‑1 (MQ=255)
            tCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGa                                                       >  1:1051907/1‑71 (MQ=255)
             cTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc                                                      >  1:1257131/1‑71 (MQ=255)
              tggtggTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGAc                                                            >  1:2466527/1‑64 (MQ=255)
                gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt                                                    <  1:3694872/70‑1 (MQ=255)
                gtggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg                                                   >  1:1699652/1‑71 (MQ=255)
                 tggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg                                                  <  1:2478699/71‑1 (MQ=255)
                  ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg                                                  >  1:2571116/1‑70 (MQ=255)
                  ggtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTgg                                                  >  1:2964215/1‑70 (MQ=255)
                   gtTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGt                                                 <  1:1800876/70‑1 (MQ=255)
                      aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg                                                   >  1:2002814/1‑65 (MQ=255)
                      aCGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTg                                                   >  1:1969392/1‑65 (MQ=255)
                       cGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaa                                            <  1:1824569/71‑1 (MQ=255)
                        gCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaa                                           <  1:1861217/71‑1 (MQ=255)
                        gCTCGACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa                                           <  1:2948002/71‑1 (MQ=255)
                          tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa                                           >  1:3614129/3‑69 (MQ=255)
                          tcgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaa                                            >  1:1667589/3‑68 (MQ=255)
                          tGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGAc                                                      >  1:1817918/1‑58 (MQ=255)
                           cgACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGACACCTGGTACGaaa                                           >  1:2271443/2‑68 (MQ=255)
                            gACTTCCGTCTGTCGAGCACGTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaat                                        <  1:1768595/70‑2 (MQ=255)
                            gACTTCCGCATGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAga                                       >  1:1895780/1‑71 (MQ=255)
                               ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt                                                    <  1:2947024/55‑1 (MQ=255)
                               ttCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCt                                                    <  1:1625578/55‑1 (MQ=255)
                                      tGTCCAGTACCTGCCTGTTCTCCGATATCGTTCTGCCCACCGCCACCTGGTACGaaa                                           <  1:2073247/57‑1 (MQ=255)
                                        tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGaa                              <  1:3617365/68‑1 (MQ=255)
                                        tCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATAc                           >  1:1785461/1‑71 (MQ=255)
                                          gAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACtt                         >  1:1507433/1‑71 (MQ=255)
                                           aGCACCTGTCTCTATTCCGACATCATTTTGCCGACgg                                                          >  1:2317658/1‑37 (MQ=255)
                                                 tGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATg                                <  1:2630184/57‑1 (MQ=255)
                                                  gTCTCTATTCCGACATCATTTTGCCGACGGCGAcc                                                     <  1:1438584/35‑1 (MQ=255)
                                                   actctATTCCGACATCATTTTGCCGACGGCTACCTg                                                   >  1:3879545/2‑36 (MQ=38)
                                                    ctctATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCa               >  1:1652325/1‑71 (MQ=255)
                                                       ttctCCGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCAt              <  1:2153056/66‑1 (MQ=255)
                                                         ttCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGaa                              <  1:601654/51‑1 (MQ=255)
                                                         ttCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCg           <  1:2145552/70‑1 (MQ=255)
                                                          tCCGACATCATTTTGCCGACGGCGACCTGGTACGaaaaa                                         >  1:1879485/1‑39 (MQ=255)
                                                          tCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGtt         >  1:729285/1‑71 (MQ=255)
                                                           ccGACATCATTTTGCCGACGGCGACCTGGTACGaa                                            >  1:1191007/1‑35 (MQ=255)
                                                           ccGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATc             >  1:406587/1‑66 (MQ=255)
                                                            cGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGa                     >  1:2581712/1‑57 (MQ=255)
                                                            cGATATCGTTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTa       <  1:1145252/71‑1 (MQ=255)
                                                            cGATATCGTTCTGCCCACCGCCACCTGTACGAAAAAGACGATATGAACACCTCGGATATGc                >  1:684424/1‑61 (MQ=255)
                                                             gACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGt          <  1:3459817/67‑1 (MQ=255)
                                                             gACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGt          <  1:2051805/67‑1 (MQ=255)
                                                              atatcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCAt              >  1:3230386/3‑62 (MQ=255)
                                                               tatcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTc    <  1:3029061/70‑1 (MQ=255)
                                                               tatcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTc    <  1:2742127/70‑1 (MQ=255)
                                                               tatcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGAACACCTCGGATATGCATCCGTTTATTc    <  1:2565881/70‑1 (MQ=255)
                                                               catcatTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTc    <  1:720934/71‑1 (MQ=255)
                                                               catcatTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTc    <  1:169395/71‑1 (MQ=255)
                                                                atcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGaa                              >  1:2818811/1‑44 (MQ=16)
                                                                atcgtTCTGCCCACCGCCACCTGGTACGAAAAAGACGATATGaa                              >  1:2478943/1‑44 (MQ=16)
                                                                atcatTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCa   <  1:171742/71‑1 (MQ=255)
                                                                  catTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGc                >  1:2837136/1‑56 (MQ=255)
                                                                  catTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAc  >  1:1547707/1‑70 (MQ=255)
                                                                    ttttGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTAtt     <  1:3586881/65‑1 (MQ=255)
                                                                      |                                                                 
AGGCAAGCTGGATCTGGTGGTTACGCTGGACTTCCGTCTGTCGAGCACCTGTCTCTATTCCGACATCATTTTGCCGACGGCGACCTGGTACGAAAAAGACGACATGAATACTTCGGATATGCATCCGTTTATTCAC  >  minE/833831‑833966

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: