Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,803 C→T 9.5% L1087L (CTC→CTT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,8030CT 9.5% 224.3 / 9.0 73L1087L (CTC→CTTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (24/42);  new base T (4/3);  total (28/45)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGC  >  minE/834738‑834867
                                                                 |                                                                
aTCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACtt                                                             <  1:833509/71‑1 (MQ=255)
 tCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGTGCTCAACTTc                                                            <  1:1367901/71‑1 (MQ=255)
   gACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACtt                                                             <  1:602333/68‑1 (MQ=255)
     cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTgc                                                        >  1:2390213/1‑70 (MQ=255)
     cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGa                                                        >  1:3660710/1‑71 (MQ=255)
     cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGa                                                        >  1:3792761/1‑71 (MQ=255)
     cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGa                                                        >  1:1291459/1‑71 (MQ=255)
     cacCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGa                                                        >  1:565879/1‑71 (MQ=255)
        ccGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGAcgc                                                     <  1:503402/71‑1 (MQ=255)
         cGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt                                                          <  1:1703389/65‑1 (MQ=255)
           ttCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACt                                                              <  1:2147493/59‑1 (MQ=255)
             cGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCAc                                                >  1:2487440/1‑71 (MQ=255)
                tGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTc                                                                  <  1:1608274/50‑1 (MQ=255)
                tGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAg                                             <  1:2722209/71‑1 (MQ=255)
                tGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAg                                             <  1:2717274/71‑1 (MQ=255)
                  aaaGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACtt                                                             <  1:160227/53‑1 (MQ=255)
                  aaaGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAg                                             >  1:3056317/1‑69 (MQ=255)
                  aaaGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa                                            <  1:2131026/70‑1 (MQ=255)
                    agaagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTc                                                            >  1:2259209/1‑52 (MQ=255)
                    agaagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTc                                                            >  1:1664772/1‑52 (MQ=255)
                    agaagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTc                                                            >  1:4116999/1‑52 (MQ=255)
                      aagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGAcgccgc                                                  <  1:839932/60‑1 (MQ=255)
                      aagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGAcgc                                                     >  1:3208786/1‑57 (MQ=255)
                      aagTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGAcg                                                      >  1:2007626/1‑56 (MQ=255)
                       agTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTggg                                      <  1:186127/71‑1 (MQ=255)
                           aTAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt                                                          <  1:2478335/47‑1 (MQ=255)
                           aTAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt                                                          <  1:2482981/47‑1 (MQ=255)
                             aGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACt                                                              <  1:2603046/41‑1 (MQ=255)
                             aGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAt                                  >  1:1689175/1‑69 (MQ=255)
                              ggCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt                                                          <  1:1525185/44‑1 (MQ=255)
                               gCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtccac                              <  1:182999/71‑1 (MQ=255)
                                ccAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa                                            >  1:1519947/1‑56 (MQ=255)
                                 cAGAAATCCAACGGCAACCAGGAAAAAGCGCTCaa                                                                >  1:3092901/1‑35 (MQ=255)
                                 cAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTc                                                            <  1:912205/39‑1 (MQ=255)
                                  aGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt                                                          <  1:4116763/40‑1 (MQ=255)
                                   gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCt                                                          <  1:3454035/39‑1 (MQ=255)
                                   gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTggg                                      <  1:2273653/59‑1 (MQ=255)
                                   gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc                                 <  1:2474456/64‑1 (MQ=255)
                                   gAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc                                 <  1:3055214/64‑1 (MQ=255)
                                    aaaTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCCGAAGTGGGGTAtccactcca                          <  1:3052458/70‑1 (MQ=255)
                                     aaTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCAcc                        >  1:1555038/1‑71 (MQ=255)
                                      aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCAc                                                <  1:1179471/46‑1 (MQ=255)
                                      aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAg                                          <  1:2562332/52‑1 (MQ=255)
                                      aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAt                                  <  1:1423086/60‑1 (MQ=255)
                                      aTCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCt                       <  1:583933/71‑1 (MQ=255)
                                       tCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa                                            <  1:615640/49‑1 (MQ=255)
                                       tCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa                                            <  1:990397/49‑1 (MQ=255)
                                        caaaCGGCTTCCCGGAAAAAGCACTTAACTTCCTGACGCCGCACCAGa                                            <  1:1225982/46‑1 (MQ=19)
                                         cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGaa                                           >  1:2073744/1‑48 (MQ=255)
                                         cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTa                      >  1:2180752/1‑69 (MQ=255)
                                         cAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTAc                     <  1:1052559/70‑1 (MQ=255)
                                          aaCGGCTTCCCGGAAAAAGCACTTAACTTCCTGACGCCGCACCAGaaa                                          >  1:3225533/1‑47 (MQ=21)
                                          aaCGGCTTCCCGGAAAAAGCACTTAACTTCCTGACGCCGCACCAGaaa                                          >  1:3818476/1‑47 (MQ=21)
                                           aCGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtcc                                >  1:988776/1‑57 (MQ=255)
                                           aCGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGc                  >  1:3779071/1‑71 (MQ=255)
                                             ggCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGa                >  1:2348567/1‑71 (MQ=255)
                                              gCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAg                                          <  1:2806355/44‑1 (MQ=255)
                                              gCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAg                                          <  1:4141499/44‑1 (MQ=255)
                                               cAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa              <  1:1441309/71‑1 (MQ=255)
                                                  cccggAAAAAGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAAcc           >  1:1675368/4‑71 (MQ=255)
                                                    aGGAAAAAGCGCTCAACTTCCTGACGCCGAACCAGAAGTgg                                       >  1:2864673/1‑41 (MQ=255)
                                                     ggAAAAAGCGCTCAACTTCCTGACGCCGCACCAGa                                            <  1:1840370/35‑1 (MQ=255)
                                                      gAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct       <  1:1600553/71‑1 (MQ=255)
                                                       aaaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa              <  1:1109918/63‑1 (MQ=255)
                                                       aaaaaGCACTTAACTTCCTGACGCCGCACCAGAAATggggca                                   >  1:1583382/1‑40 (MQ=25)
                                                        aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgctg      >  1:2656707/1‑70 (MQ=255)
                                                        aaaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCtgctga     >  1:3863982/1‑71 (MQ=255)
                                                         aaaGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct       <  1:1595360/68‑1 (MQ=255)
                                                           aGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACa              <  1:1057953/59‑1 (MQ=255)
                                                           aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc  <  1:1257276/71‑1 (MQ=255)
                                                           aGCACTTAACTTCCTGACGCCGCACCAGAAATGGGGCATTCACTCAACCTACAGTGAAAACCtgctaatgc  <  1:259185/71‑1 (MQ=255)
                                                            gcgcTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACctgct       <  1:2611229/65‑1 (MQ=255)
                                                               cTCAACTTCCTGACGCCGCACCAGAAGTGGGGTAtc                                 <  1:350975/36‑1 (MQ=255)
                                                               cTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCg                 <  1:2377482/52‑1 (MQ=255)
                                                                 cAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAAc            >  1:1835006/1‑55 (MQ=255)
                                                                 |                                                                
ATCGACACCCGTTCGGTGAAAGAAGTGATAGGCCAGAAATCCAACGGCAACCAGGAAAAAGCGCTCAACTTCCTGACGCCGCACCAGAAGTGGGGTATCCACTCCACCTACAGCGACAACCTGCTGATGC  >  minE/834738‑834867

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: