Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,259,631 A→G 7.2% R431R (CGT→CGC torZ ← trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,259,6310AG 7.2% 193.7 / 4.0 69R431R (CGT→CGCtorZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base A (36/28);  new base G (3/2);  total (39/30)
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

TTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCG  >  minE/1259563‑1259695
                                                                    |                                                                
ttCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATAc                                                                 <  1:3431187/70‑1 (MQ=255)
   ggATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGc                                                             <  1:750469/71‑1 (MQ=255)
          tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                                  <  1:1396129/59‑1 (MQ=255)
          tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                                  <  1:1846398/59‑1 (MQ=255)
           ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGaa                                                     >  1:2301660/1‑71 (MQ=255)
            tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCaa                                                           >  1:4416331/1‑64 (MQ=255)
            tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:5262237/1‑71 (MQ=255)
            tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:850539/1‑71 (MQ=255)
            tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:567388/1‑71 (MQ=255)
             gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCATCAGGAATc                                                   >  1:168942/1‑71 (MQ=255)
             gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    <  1:4522051/70‑1 (MQ=255)
             gTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATc                                                   >  1:3700420/1‑71 (MQ=255)
              ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                                      <  1:5199668/67‑1 (MQ=255)
                cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:5448335/1‑67 (MQ=255)
                cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:429945/1‑67 (MQ=255)
                cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGcc                                                >  1:2998410/1‑71 (MQ=255)
                  ttACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              >  1:128213/1‑71 (MQ=255)
                   tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:4964833/70‑1 (MQ=255)
                   tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:4261229/70‑1 (MQ=255)
                   tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              >  1:1514098/1‑70 (MQ=255)
                   tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:340776/70‑1 (MQ=255)
                     ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAggg                                                      >  1:3145899/1‑59 (MQ=255)
                     ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAggg                                                      >  1:5087717/1‑59 (MQ=255)
                     ccGTTGTGTTGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAttt                                                  >  1:1695375/1‑62 (MQ=255)
                          aTGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCg                                       <  1:5195118/70‑1 (MQ=255)
                           tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                                      >  1:4152849/1‑70 (MQ=255)
                           tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                                      >  1:1056701/1‑70 (MQ=255)
                             ttGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAt                                                    <  1:5642868/53‑1 (MQ=255)
                              tGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATa                                                                  <  1:3528062/39‑1 (MQ=255)
                              tGATATGCGCCACCAGGGTTTTCCAGTGCTTCAACAATGCGGGCAACAGGGAt                                                    <  1:5274536/53‑1 (MQ=255)
                                atatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatca                                 <  1:5301244/70‑1 (MQ=255)
                                 tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                                         >  1:3113091/1‑61 (MQ=255)
                                 tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                                         >  1:4898602/1‑61 (MQ=255)
                                 tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc                                         >  1:4371361/1‑61 (MQ=255)
                                 tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                >  1:489677/1‑70 (MQ=255)
                                   tttGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGaa                                                     >  1:4554316/1‑47 (MQ=255)
                                   tttGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                                      <  1:1910830/46‑1 (MQ=255)
                                     tGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCg                           >  1:332939/1‑71 (MQ=255)
                                       ccgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGc                                      <  1:476519/58‑1 (MQ=255)
                                        cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                >  1:1364440/1‑63 (MQ=255)
                                        cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc                        <  1:4989511/71‑1 (MQ=255)
                                        cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc                        <  1:1653200/71‑1 (MQ=255)
                                        cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGcc                        >  1:82557/1‑71 (MQ=255)
                                            aGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGcc                                                <  1:4784928/43‑1 (MQ=255)
                                            aGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGc                    >  1:2855123/1‑71 (MQ=255)
                                              gATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                                      <  1:3863570/35‑1 (MQ=255)
                                              gATTTTCCAGCGCATCGACGATACGGGCAACAGGa                                                      <  1:1339194/35‑1 (MQ=255)
                                              gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:4889252/43‑1 (MQ=255)
                                              gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:5176533/43‑1 (MQ=255)
                                              gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:5182770/43‑1 (MQ=255)
                                              gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt                                              <  1:1315152/43‑1 (MQ=255)
                                              gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCatcatc                               <  1:4849533/58‑1 (MQ=255)
                                               aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTc                                             <  1:2801507/43‑1 (MQ=255)
                                               aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTc                                             <  1:782481/43‑1 (MQ=255)
                                               aTTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTgg                 <  1:210047/71‑1 (MQ=255)
                                                ttttCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:4680944/1‑35 (MQ=255)
                                                ttttCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:1220249/1‑35 (MQ=255)
                                                ttttCCAGCGCATCGACGATACGGGCAACAGGAAt                                                    >  1:1115151/1‑35 (MQ=255)
                                                ttttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAg                             >  1:4939809/1‑58 (MQ=255)
                                                ttttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAg                             >  1:1573795/1‑58 (MQ=255)
                                                 tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                <  1:634346/54‑1 (MQ=255)
                                                 tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                <  1:895565/54‑1 (MQ=255)
                                                 tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                <  1:3976185/54‑1 (MQ=255)
                                                 tttCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                <  1:1938859/54‑1 (MQ=255)
                                                   tCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCa                                            >  1:3915666/1‑40 (MQ=255)
                                                   tCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCa                                            >  1:1238459/1‑40 (MQ=255)
                                                      aGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTcc          >  1:1461784/1‑71 (MQ=255)
                                                      aGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTcc          >  1:1176889/1‑71 (MQ=255)
                                                      aGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTc           >  1:4970305/1‑70 (MQ=255)
                                                              gacgaTACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATcg  >  1:4138658/1‑71 (MQ=255)
                                                               acgaTACGGGCAACAGGAATCGCCGTCATTCCGCCAtcatca                              >  1:5117048/1‑42 (MQ=255)
                                                               acgaTACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAAt    >  1:4219019/1‑68 (MQ=255)
                                                                    |                                                                
TTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCG  >  minE/1259563‑1259695

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