Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,515 A→G 42.9% E991E (GAA→GAG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,5150AG42.9% 14.1 / 47.8 42E991E (GAA→GAGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (10/14);  new base G (10/8);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTC  >  minE/834447‑834583
                                                                    |                                                                    
gTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAagaaga                                                                    <  1:2445255/71‑1 (MQ=255)
   aaaGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACg                                                                  <  1:1527250/70‑1 (MQ=255)
        cTGGGCTGCCCTGAGCGAATTTACCGGTCGTGCCCATACGCATCTGGCGCTGAATAAAGAAGACgag                                                                >  1:800475/1‑67 (MQ=255)
                cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGAt                                                           >  1:2332549/1‑64 (MQ=255)
                cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCtt                                                     <  1:1227682/70‑1 (MQ=255)
                cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCtt                                                     <  1:2870825/70‑1 (MQ=255)
                cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCt                                                      <  1:3683722/69‑1 (MQ=255)
                 ccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCtt                                                     <  1:485273/69‑1 (MQ=255)
                   tGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCcg                                                  >  1:243645/1‑70 (MQ=255)
                     aGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa                                               >  1:288234/1‑71 (MQ=255)
                      gCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGa                                               >  1:210116/1‑70 (MQ=255)
                       cGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGat                                              <  1:2344688/70‑1 (MQ=255)
                            tcaCCGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAg                                        <  1:544803/69‑1 (MQ=255)
                             caCCGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAgg                                       >  1:1978356/2‑71 (MQ=255)
                             caCCGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAgg                                       >  1:833000/2‑71 (MQ=255)
                               ccGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATAtt                                           <  1:3300136/65‑1 (MQ=255)
                               ccGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCt                                                      <  1:1319235/54‑1 (MQ=255)
                               ccGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg                                     <  1:1850715/71‑2 (MQ=255)
                               ccGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg                                     <  1:525403/71‑2 (MQ=255)
                                cGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGatat                                            <  1:3829395/63‑1 (MQ=255)
                                cGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcac                                    <  1:708986/71‑1 (MQ=255)
                                  gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGAc                                                                   <  1:3301163/38‑1 (MQ=255)
                                  gTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAg                                  <  1:3711096/71‑1 (MQ=255)
                                    cGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGcc                                >  1:1505855/1‑71 (MQ=255)
                                       gaacATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgcg                             >  1:2012385/4‑71 (MQ=255)
                                        aCCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCcgcg                             >  1:3373149/1‑70 (MQ=255)
                                          cATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCa                                         >  1:1686105/1‑56 (MQ=255)
                                           aTACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCcgc                                                 >  1:2100489/1‑47 (MQ=255)
                                           aTACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCaaa                         <  1:444828/71‑1 (MQ=255)
                                           aTACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAg                                        >  1:1265997/1‑56 (MQ=255)
                                              cccATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAAt                       >  1:726047/3‑70 (MQ=255)
                                               gCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTa                     <  1:1798718/71‑1 (MQ=255)
                                               ccATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCcgc                              >  1:71946/2‑62 (MQ=255)
                                                cATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTAt                    <  1:501873/71‑1 (MQ=255)
                                                cATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTAt                    <  1:2869360/71‑1 (MQ=255)
                                                 aTCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAAtt                      >  1:2402586/1‑68 (MQ=255)
                                                 aTCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTAtc                   <  1:1834838/71‑1 (MQ=255)
                                                  tCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTcgc                                                 <  1:461865/40‑1 (MQ=255)
                                                    tGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAgg                                       >  1:3407723/1‑48 (MQ=255)
                                                    tGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAgg                                       >  1:3069163/1‑48 (MQ=255)
                                                      gcgcTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTcc                <  1:715014/69‑1 (MQ=255)
                                                       cgcTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGc             >  1:3406273/1‑71 (MQ=255)
                                                           gAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGa         >  1:1145938/1‑71 (MQ=255)
                                                                aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCAc        >  1:3447536/1‑67 (MQ=255)
                                                                aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctgg    <  1:946677/71‑1 (MQ=255)
                                                                aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctgg    >  1:2264364/1‑71 (MQ=255)
                                                                aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctgg    >  1:1272132/1‑71 (MQ=255)
                                                                aagaagACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctgg    >  1:1483700/1‑71 (MQ=255)
                                                                 agaggaCGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctggt   >  1:3512271/1‑71 (MQ=255)
                                                                 agaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTAtctc                 <  1:1289726/57‑1 (MQ=255)
                                                                  gaagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtc  >  1:1206579/1‑71 (MQ=255)
                                                                    |                                                                    
GTGAAAGCCTGGGCTGCCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTC  >  minE/834447‑834583

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: