Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 836,102 T→C 11.7% L274L (TTG→CTG) ‡ narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE836,1020TC11.7% 171.9 / 12.0 68L274L (TTG→CTG) ‡narHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base T (40/20);  new base C (2/6);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 4.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCA  >  minE/836033‑836169
                                                                     |                                                                   
aTTGAAGCGGGTCAGACGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgtt                                                                    >  1:2292953/1‑71 (MQ=255)
 ttGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                                   >  1:5677282/1‑71 (MQ=255)
  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                                   >  1:866778/1‑70 (MQ=255)
  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                                   <  1:2321273/70‑1 (MQ=255)
  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                                   >  1:3607389/1‑70 (MQ=255)
  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                                   >  1:3838078/1‑70 (MQ=255)
  tGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCtgttg                                                                   >  1:4447281/1‑70 (MQ=255)
   gAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGAGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa                                                                 >  1:269440/1‑71 (MQ=255)
    aaGCGTGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa                                                                 <  1:3936445/70‑1 (MQ=255)
     aGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACg                                                               <  1:4309061/71‑1 (MQ=255)
      tccggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGa                                                              >  1:276075/4‑71 (MQ=255)
      gCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                              >  1:4667558/1‑71 (MQ=255)
        cggTCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGac                                                             <  1:812112/69‑2 (MQ=255)
          gTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTAc                                                                <  1:2662269/65‑1 (MQ=255)
           tCAACCGACCGTGTGCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGACGCCg                                                         <  1:5551407/71‑1 (MQ=255)
           tCAACCGACCGTGTGCTCAGAAACCTGCATGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTACGACGcc                                                          >  1:2781159/1‑70 (MQ=255)
            cAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGa                                                        >  1:1801949/1‑71 (MQ=255)
             aGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGAc                                                       >  1:4856313/1‑71 (MQ=255)
             aGCAGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGAc                                                       >  1:4834711/1‑71 (MQ=255)
              gCCGACCGTGTGATCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGcc                                                          <  1:1669589/67‑1 (MQ=255)
               ccGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTAc                                                                <  1:1369153/60‑1 (MQ=255)
                cGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                                                     >  1:125299/1‑70 (MQ=255)
                cGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                                                     <  1:3257921/70‑1 (MQ=255)
                    cGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                              >  1:1920537/1‑57 (MQ=255)
                       gtgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAAc                                             >  1:1712571/1‑71 (MQ=255)
                       gtgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                                                <  1:1753677/52‑1 (MQ=255)
                       gtgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                                                <  1:1364331/52‑1 (MQ=255)
                       gtgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                                                <  1:4624859/52‑1 (MQ=255)
                        tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                              >  1:1093989/1‑53 (MQ=255)
                        tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                              >  1:5142558/1‑53 (MQ=255)
                        tgCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGa                                                              >  1:1306591/1‑53 (MQ=255)
                        tgCTCAGAAACCTGCGTGGGTCGCATCCGGTATCTGGGCGTGCTGCTTTAc                                                                <  1:1447620/51‑1 (MQ=255)
                          cTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTg                                          <  1:5285365/71‑1 (MQ=255)
                             aGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTa                                                                 >  1:5124531/1‑45 (MQ=255)
                             aGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGc                                         <  1:447159/69‑1 (MQ=255)
                                aaCCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagc                                      >  1:1542849/1‑69 (MQ=255)
                                 aCCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagccag                                   <  1:1057236/71‑1 (MQ=255)
                                 aCCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagccag                                   <  1:4551761/71‑1 (MQ=255)
                                   cTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGCCGCTATTGAACGTGCAGCCAGCa                                 >  1:3687299/1‑71 (MQ=255)
                                   cTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCa                                 >  1:5308503/1‑71 (MQ=255)
                                    tgtgTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCAc                                >  1:516212/1‑71 (MQ=255)
                                     gtgtCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGcagc                                      >  1:1822240/1‑64 (MQ=255)
                                      tgtCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCAcc                               >  1:1109144/1‑70 (MQ=255)
                                           gtcgTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                                                     <  1:4986321/43‑1 (MQ=255)
                                           gtcgTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGc                                                     <  1:4479607/43‑1 (MQ=255)
                                             cgTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac                         >  1:4635675/1‑69 (MQ=255)
                                             cgTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac                         >  1:1731161/1‑69 (MQ=255)
                                              gTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac                         >  1:4330242/1‑68 (MQ=255)
                                               tATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaacgaga                     <  1:5545138/71‑1 (MQ=255)
                                                aTCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAAc                                             >  1:4102071/1‑46 (MQ=255)
                                                aTCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaacgaga                     >  1:3000351/1‑70 (MQ=255)
                                                 tCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaacgaga                     >  1:3121/1‑69 (MQ=255)
                                                 tCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac                         >  1:84392/1‑65 (MQ=255)
                                                 tCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgagaac                         >  1:5270669/1‑65 (MQ=255)
                                                   cGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACcgag                            >  1:2697642/1‑60 (MQ=255)
                                                     ttATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATc               >  1:1709170/1‑71 (MQ=255)
                                                     ttATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATc               >  1:1152525/1‑71 (MQ=255)
                                                     ttATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATc               >  1:3305631/1‑71 (MQ=255)
                                                       aTCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGt                                           <  1:183513/41‑1 (MQ=255)
                                                       aTCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGc                                         >  1:4633019/1‑43 (MQ=255)
                                                       aTCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGc                                         >  1:3194964/1‑43 (MQ=255)
                                                       aTCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCtt             >  1:5082793/1‑71 (MQ=255)
                                                       aTCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCtt             >  1:4232619/1‑71 (MQ=255)
                                                       aTCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCtt             >  1:4687490/1‑71 (MQ=255)
                                                         cTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTa           >  1:5420955/1‑71 (MQ=255)
                                                          ttGGCGTGCTGTTGTACGATGCCGACGCTATTGAACg                                            <  1:471528/37‑1 (MQ=255)
                                                               gTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCtt             >  1:1987960/1‑63 (MQ=255)
                                                                 gCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGAt                >  1:4725820/1‑58 (MQ=255)
                                                                 gCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGAt                >  1:1334109/1‑58 (MQ=255)
                                                                 gCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGAt                >  1:2100409/1‑58 (MQ=255)
                                                                  cTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGt    <  1:39227/69‑1 (MQ=255)
                                                                  cTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGt    <  1:3024941/69‑1 (MQ=255)
                                                                  cTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGt    <  1:4242686/69‑1 (MQ=255)
                                                                  cTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGt    <  1:968287/69‑1 (MQ=255)
                                                                  cTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCa  <  1:5328740/71‑1 (MQ=255)
                                                                    gttgtACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCa  >  1:1292840/1‑69 (MQ=255)
                                                                     ttgtACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCttt            <  1:5327897/58‑1 (MQ=255)
                                                                     ttgtACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCtt             <  1:2444341/57‑1 (MQ=255)
                                                                     ttgtACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCtt             <  1:5422309/57‑1 (MQ=255)
                                                                     |                                                                   
ATTGAAGCGGGTCAGCCGACCGTGTGCTCAGAAACCTGTGTCGGTCGTATCCGTTATCTTGGCGTGCTGTTGTACGATGCCGACGCTATTGAACGTGCAGCCAGCACCGAGAACGAGAAAGATCTTTACCAGCGTCA  >  minE/836033‑836169

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: