Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,584 T→C 23.2% S1014S (TCT→TCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,5840TC23.2% 72.3 / 19.6 43S1014S (TCT→TCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (21/12);  new base C (7/3);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCG  >  minE/834514‑834652
                                                                      |                                                                    
aagaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtct                                                                      >  1:3823411/1‑71 (MQ=255)
   aCGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctggtccggt                                                                   <  1:4856072/71‑1 (MQ=255)
   aCGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctggtccggt                                                                   <  1:1970290/71‑1 (MQ=255)
   aCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctggt                                                                   >  1:3877035/1‑71 (MQ=255)
    cGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCtggtctggtc                                                                  <  1:3241924/71‑1 (MQ=255)
    cGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCtggtctggtc                                                                  <  1:237431/71‑1 (MQ=255)
     gagaAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctgg                                                                    <  1:2699864/68‑1 (MQ=255)
         agaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGtctggtctgg                                                               <  1:1168819/69‑1 (MQ=255)
          gaTTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGaaa                                                            <  1:3602679/71‑2 (MQ=255)
           aTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGtctcgtctgg                                                               <  1:728429/67‑1 (MQ=255)
            tCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt                                                          <  1:2317952/71‑1 (MQ=255)
             tcGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGaaag                                                           >  1:445746/2‑67 (MQ=255)
             tcGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGaaag                                                           >  1:4281304/2‑67 (MQ=255)
             tcGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGaaag                                                           >  1:4075744/2‑67 (MQ=255)
                    cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACg                                                   <  1:1001419/70‑1 (MQ=255)
                    cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGt                                                  >  1:4414918/1‑71 (MQ=255)
                    cgcgATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGt                                                  >  1:4343050/1‑71 (MQ=255)
                       gATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTc                                               <  1:1515092/71‑1 (MQ=255)
                        atatTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCt                                              <  1:4720274/71‑1 (MQ=255)
                           ttCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTa                                            >  1:3980592/1‑70 (MQ=255)
                             cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCtggtctggtc                                                                  >  1:562999/1‑46 (MQ=255)
                             cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACaa                                         >  1:913106/1‑71 (MQ=255)
                             cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                          <  1:19916/70‑1 (MQ=255)
                             cAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                          <  1:3276936/70‑1 (MQ=255)
                              aGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTggtctggtct                                                                 >  1:4521034/1‑46 (MQ=255)
                              aGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGgtctggtctg                                                                >  1:4440270/1‑47 (MQ=255)
                                gCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt                                                          >  1:3828071/1‑51 (MQ=255)
                                 cacaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGcc                                     >  1:5485810/1‑71 (MQ=255)
                                  gcaGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGcgg                                    <  1:1426495/70‑3 (MQ=255)
                                  acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                                    >  1:2810182/1‑71 (MQ=255)
                                  acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                                    >  1:775879/1‑71 (MQ=255)
                                  acaGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCg                                    >  1:3431265/1‑71 (MQ=255)
                                    aGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGt                                  <  1:4438733/71‑1 (MQ=255)
                                    aGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGt                                  <  1:3199465/71‑1 (MQ=255)
                                      ccGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa                                                       >  1:314518/1‑48 (MQ=255)
                                      ccGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa                                                       >  1:1951803/1‑48 (MQ=255)
                                          gcAAAATTATCTCCAGCCCGCCCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac                             >  1:1316838/1‑70 (MQ=255)
                                          gcAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTacac                             >  1:4136311/1‑70 (MQ=255)
                                           taaaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACAcc                            >  1:2529333/2‑70 (MQ=255)
                                           cAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACa                                          >  1:1576026/1‑56 (MQ=255)
                                            aaaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCggg                                   >  1:2979045/1‑62 (MQ=255)
                                             aaaTTATCTCCAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAAc                         >  1:692219/1‑71 (MQ=255)
                                              aaTTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACg                        >  1:5147618/1‑71 (MQ=255)
                                                ttATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGtt                      <  1:4566101/71‑1 (MQ=255)
                                                   tctcCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCAc                   >  1:422689/1‑71 (MQ=255)
                                                   tctcCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCAc                   >  1:3109015/1‑71 (MQ=255)
                                                       cAGCCCGACCTGGTCTGGTCTGGAAGATGAACATGTTTCTTACAACGCCGGTTACACCAACGTTCAc                   <  1:3874263/67‑1 (MQ=255)
                                                       cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAga               <  1:3325402/71‑2 (MQ=255)
                                                       cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc               >  1:2279472/1‑71 (MQ=255)
                                                       cAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGc               >  1:3771642/1‑71 (MQ=255)
                                                       cAGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAgt               >  1:3245122/1‑70 (MQ=255)
                                                        aGCCCCACCTGGTCCGGTCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAgtt              <  1:3154329/71‑3 (MQ=255)
                                                            cGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGa                 >  1:5564691/1‑64 (MQ=255)
                                                              aCCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATcc         >  1:3053275/1‑70 (MQ=255)
                                                               ccTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa       <  1:2859978/71‑1 (MQ=255)
                                                                ctggtccggtCTGGAAAGCTATCACGTCTCCTATAACGCGGGATACACCAACGtt                      <  1:1773061/55‑1 (MQ=255)
                                                                ctggtccggtCTGGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGtt                      <  1:4670662/55‑1 (MQ=255)
                                                                 tggtctggtcTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATg     >  1:4345503/1‑71 (MQ=255)
                                                                   gtctggtctgGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTg             <  1:4221275/61‑1 (MQ=255)
                                                                   gtctggtctgGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGc   >  1:1666562/1‑71 (MQ=255)
                                                                   gtccggtctgGAAAGCGATCACGTCTCCTATAACGCGGGATACACCAACGTTCATGAGTTAATTCCGTgg    <  1:4443603/70‑1 (MQ=255)
                                                                     ctggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCg  >  1:1296049/1‑70 (MQ=255)
                                                                      tggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa       <  1:997608/64‑1 (MQ=255)
                                                                      tggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa       <  1:1256267/64‑1 (MQ=255)
                                                                      tggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa       <  1:1249383/64‑1 (MQ=255)
                                                                      tggtctggAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCa       <  1:1082990/64‑1 (MQ=255)
                                                                      |                                                                    
AAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACACGTTTCTTACAACGCCGGTTACACCAACGTTCACGAGCTGATCCCATGGCG  >  minE/834514‑834652

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: