Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 834,527 C→T 45.9% I995I (ATC→ATT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE834,5270CT45.9% 16.5 / 76.5 61I995I (ATC→ATTnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (18/15);  new base T (17/11);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

CCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAT  >  minE/834463‑834596
                                                                |                                                                     
cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCtt                                                                  <  1:2870825/70‑1 (MQ=255)
cccTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCtt                                                                  <  1:1227682/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:833000/2‑71 (MQ=255)
             caCCGGGCGCGAACATACCCATCTGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAgg                                                    >  1:1978356/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)
                  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:3069163/1‑48 (MQ=255)
                                    tGGCGCTGCACAAAGAGGACGAGAAGATTCGCTTTCGCGATATTCAgg                                                    >  1:3407723/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:1272132/1‑71 (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:3154823/1‑71 (MQ=255)
                                                     gaCGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACctggtctgg            >  1:92992/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:719255/1‑40 (MQ=255)
                                                        gagaAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCa                                        >  1:795493/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)
                                                                |                                                                     
CCCTGAGCGAATTTACCGGTCGTGACCATACGCATCTGGCGCTGAATAAAGAAGACGAGAAGATCCGCTTCCGCGATATTCAGGCACAGCCGCGCAAAATTATCTCCAGCCCGACCTGGTCTGGTCTGGAAGAT  >  minE/834463‑834596

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: