Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,059,436 (G)6→7 coding (181/2547 nt) torA → trimethylamine N‑oxide reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,059,4301.G92.6% 85.7 / 0.5 27coding (175/2547 nt)torAtrimethylamine N‑oxide reductase 1
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (13/12);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCG  >  NC_000913/1059315‑1059557
                                                                                                                     |                                                                                                                                  
cTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcgc                                                                                                                 >  2:151370/1‑133 (MQ=255)
 ttAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcgcg                                                                                                                <  2:89879/133‑1 (MQ=255)
  tgaCCGTCGCTGGGATGCTGGGGGCGGCCTTGGTAACGCCGAGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGGCGCACTGGGGGGGCTATCCGCGCGa                                                                                                               <  2:291695/131‑1 (MQ=255)
          gCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCGTACCGGGTCGCACTGGGGGGGCTATCCGCGTGACGGTGAAg                                                                                                       >  2:253596/1‑133 (MQ=255)
                        gCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTtgt                                                                                         <  2:16794/133‑1 (MQ=255)
                         ccGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTtgtg                                                                                        >  1:66644/1‑133 (MQ=255)
                         ccGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTtgtg                                                                                        <  1:33765/133‑1 (MQ=255)
                           gTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTggc                                                                                      >  2:497780/1‑133 (MQ=255)
                             cATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTggcgg                                                                                    <  1:151370/133‑1 (MQ=255)
                             cATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGCGGCTATCCGCGCGACGGTGAAGGATTGTCGGTTTGTggcgg                                                                                    >  2:161903/1‑133 (MQ=255)
                                         cgcgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg                                                                        >  1:395797/1‑133 (MQ=255)
                                         cgcgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg                                                                        >  2:366624/1‑133 (MQ=255)
                                         cgcgACGTACGACTGAGGTGCAAGCGGCGAATGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg                                                                        <  1:161903/133‑1 (MQ=255)
                                           cgACGTGCCACTGCGTCGCAAGTGGCGACTCACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGaa                                                                      <  1:173386/133‑1 (MQ=255)
                                                tGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGGCGCGTTGGGGCGGCAAGACCGGTCGAACTGGa                                                                 >  2:173386/1‑133 (MQ=255)
                                                                     gACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgat                                            >  2:524501/1‑133 (MQ=255)
                                                                     gACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgat                                            <  1:313952/133‑1 (MQ=255)
                                                                            gCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgattgccgga                                     >  1:440087/1‑133 (MQ=255)
                                                                                tCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgattgccggattgc                                 >  1:14243/1‑133 (MQ=255)
                                                                                   tctcGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                              <  1:514829/133‑1 (MQ=255)
                                                                                   tctcGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                              <  1:366624/133‑1 (MQ=255)
                                                                                   tctcGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                              <  2:66644/133‑1 (MQ=255)
                                                                                   tctcGAAAGAGGGCATTCTTACCGGGTCGCAATGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                              <  2:440087/133‑1 (MQ=255)
                                                                                     tcGAAAGAGGACATTTTTACCGGGCCGCCCTGTGGGGGCTAACCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATtgccggattgccggat                            <  2:185129/133‑1 (MQ=255)
                                                                                        aaaGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCAc                         >  1:314642/1‑133 (MQ=255)
                                                                                               gCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTAcaca                  >  2:112585/1‑133 (MQ=255)
                                                                                                  ttCTTACCGGGTCGCACTTGGGCGGCTATCGGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAAcg               <  1:212278/133‑1 (MQ=255)
                                                                                                          gggTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGt       >  1:355416/1‑133 (MQ=255)
                                                                                                               gCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCg  <  2:314642/133‑1 (MQ=255)
                                                                                                               gCACTGGGCGGGCTATCCACGCCACGGTGACGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCg  <  1:175703/133‑1 (MQ=255)
                                                                                                                     |                                                                                                                                  
CTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCG  >  NC_000913/1059315‑1059557

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