Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,533 C→T 43.1% F997F (TTC→TTT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,5330CT43.1% 32.1 / 75.2 65F997F (TTC→TTTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (22/15);  new base T (17/11);  total (39/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 = 8.07e-01

TGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACAC  >  minE/834466‑834602
                                                                   |                                                                     
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:833000/2‑71 (MQ=255)
          caCCGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAgg                                                          >  1:1978356/2‑71 (MQ=255)
            ccGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATAtt                                                              <  1:3300136/65‑1 (MQ=255)
            ccGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg                                                        <  1:525403/71‑2 (MQ=255)
            ccGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAGGcg                                                        <  1:1850715/71‑2 (MQ=255)
             cGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGatat                                                               <  1:3829395/63‑1 (MQ=255)
             cGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGcac                                                       <  1:708986/71‑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:1272132/1‑71 (MQ=255)
                                             aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctgg                       <  1:946677/71‑1 (MQ=255)
                                             aagaggACGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACctgg                       >  1:2264364/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)
                                                  gaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctgg                  >  1:92992/1‑71 (MQ=255)
                                                  gaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctgg                  >  1:3154823/1‑71 (MQ=255)
                                                    cGAGAAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCtggtccggtc                >  1:3442243/1‑71 (MQ=255)
                                                     gagaAGATTCGCTTTCGCGATATTCAGGCGCAGCcgc                                                 <  1:584969/37‑1 (MQ=25)
                                                     gagaAGATTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCtggtccggtc                <  1:1052559/70‑1 (MQ=255)
                                                     gagaAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCa                                              >  1:795493/1‑40 (MQ=255)
                                                     gagaAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCa                                              >  1:719255/1‑40 (MQ=255)
                                                       gaagaTTCGCTTTCGCGATATTCAGGCGCAgcggcgt                                               >  1:2152428/1‑36 (MQ=16)
                                                       gaagaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGgtctggtctg              >  1:2688955/1‑70 (MQ=255)
                                                         agaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGa            <  1:2507330/70‑1 (MQ=255)
                                                          gaTTCGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAATTATCTCCAGCCCCACCTGGTCCGGTCTGGa            >  1:2016113/1‑69 (MQ=255)
                                                          gaTCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGgtctgtctgg             <  1:1862095/67‑1 (MQ=255)
                                                            tCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGaa           >  1:272111/1‑68 (MQ=255)
                                                            tCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt        >  1:2658065/1‑71 (MQ=255)
                                                            tCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAt        >  1:2720303/1‑71 (MQ=255)
                                                             tcGCTTTCGCGATATTCAGGCGCAGCCGCGTAAAAtt                                         <  1:958024/36‑1 (MQ=25)
                                                              cGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGa      <  1:579043/71‑1 (MQ=255)
                                                               gCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTAt                                       >  1:1614774/1‑37 (MQ=255)
                                                               gCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa     <  1:597578/71‑1 (MQ=255)
                                                                cTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAac    <  1:17509/71‑1 (MQ=255)
                                                                 ttCCGCGATATTCAGGCACAGCCGCGCAAAATTAt                                       <  1:3553451/35‑1 (MQ=255)
                                                                 ttCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctgg                       >  1:3934788/1‑51 (MQ=255)
                                                                  tCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGaa     >  1:3157990/1‑68 (MQ=255)
                                                                  tCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac  >  1:162863/1‑71 (MQ=255)
                                                                   ccGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac  >  1:1261589/1‑70 (MQ=255)
                                                                   ccGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAacac  <  1:250499/70‑1 (MQ=255)
                                                                   |                                                                     
TGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGATGAACAC  >  minE/834466‑834602

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: