Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,955,6930GA31.7% 157.0 / 99.3 104P439S (CCT→TCT) torZtrimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base G (37/34);  new base A (20/13);  total (57/47)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 111 total aligned reads displayed.

ATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCG  >  NC_000913/1955644‑1955742
                                                 |                                                 
atatTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAgg                                                   <  1:1502625/50‑1 (MQ=255)
 tatTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGa                                                  >  1:4010195/1‑50 (MQ=255)
 tatTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGa                                                  >  1:2001728/1‑50 (MQ=255)
  atTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAt                                                 >  2:4049744/1‑50 (MQ=255)
  atTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGAAt                                                 >  1:2780779/1‑50 (MQ=255)
   ttCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAtt                                                <  1:2527672/50‑1 (MQ=255)
   ttCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAt                                                 <  1:4594935/49‑1 (MQ=255)
   ttCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGAAtt                                                <  2:315573/50‑1 (MQ=255)
     cGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAtttt                                              >  1:1344735/1‑50 (MQ=255)
     cGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAttg                                               <  2:2663697/49‑2 (MQ=255)
     cGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGAAttt                                               <  1:2310347/49‑1 (MQ=255)
      ggATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTc                                             >  2:2176062/1‑50 (MQ=255)
      ggATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTc                                             <  1:3566073/50‑1 (MQ=255)
       gATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTcc                                            <  2:758132/50‑1 (MQ=255)
       gATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTc                                             >  1:2227890/1‑49 (MQ=255)
       gATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTc                                             >  2:3075793/1‑49 (MQ=255)
       aaTACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGAACTTTcc                                            <  2:5380479/49‑1 (MQ=255)
         tACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAg                                          <  2:1992619/50‑1 (MQ=255)
         tACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAg                                          >  2:1904647/1‑50 (MQ=255)
         tACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAg                                          <  2:1218647/50‑1 (MQ=255)
           cGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgcg                                        >  2:724444/1‑50 (MQ=255)
           cGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgc                                         <  1:109576/49‑1 (MQ=255)
            gTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgcgc                                       <  1:2945504/50‑1 (MQ=255)
            gTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgcg                                        >  1:5299817/1‑49 (MQ=255)
            gTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAgcgc                                       <  2:1145714/50‑1 (MQ=255)
            gTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAgcg                                        >  1:1375432/1‑49 (MQ=255)
             tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCa                                      >  2:3445664/1‑50 (MQ=255)
              ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCAt                                     >  1:4286642/1‑50 (MQ=255)
               tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                     <  1:3827520/49‑1 (MQ=255)
               tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCAt                                     <  1:1010905/49‑1 (MQ=255)
               tGTTCTTTACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATc                                    >  2:2599431/1‑50 (MQ=255)
                gTTCTTTACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATc                                    >  1:3325445/1‑49 (MQ=255)
                 ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                  <  2:4002929/50‑1 (MQ=255)
                 ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                  >  1:3487904/1‑50 (MQ=255)
                 ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                  >  2:4026658/1‑50 (MQ=255)
                 ttCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                  <  1:2628121/50‑1 (MQ=255)
                 ttCTTTACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATcga                                  >  1:2951992/1‑50 (MQ=255)
                  tCTTTACCGTTATGCTGATATTTTCCGCCAGAATTTTCCAGCGCATcgac                                 >  1:4030545/1‑50 (MQ=255)
                  tCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                 >  1:3629483/1‑50 (MQ=255)
                   cTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgac                                 >  2:1300336/1‑49 (MQ=255)
                    tttACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcgacga                               >  1:1588898/1‑50 (MQ=255)
                      tACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATcgacg                                <  2:5135698/47‑1 (MQ=255)
                      tACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATcgacg                                <  2:4234756/47‑1 (MQ=255)
                      tACCGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATa                             >  1:229587/1‑50 (MQ=255)
                       aCCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATAc                            >  1:3959134/1‑50 (MQ=255)
                       aCCGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATAc                            >  2:5249293/1‑50 (MQ=255)
                        ccGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATAc                            <  1:2124572/49‑1 (MQ=255)
                        ccGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATAc                            <  1:3513247/49‑1 (MQ=255)
                        ccGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATAc                            <  2:2963589/49‑1 (MQ=255)
                         cGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACgg                          <  2:3053783/50‑1 (MQ=255)
                         cGTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACgg                          <  2:4412905/50‑1 (MQ=255)
                          gTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACggg                         <  1:2864716/50‑1 (MQ=255)
                          gTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACggg                         <  1:4234415/50‑1 (MQ=255)
                          gTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACggg                         >  2:2245956/1‑50 (MQ=255)
                          gTTATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACggg                         <  2:2761311/50‑1 (MQ=255)
                           ttATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACggg                         >  1:2505039/1‑49 (MQ=255)
                           ttATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACggg                         >  2:4910065/1‑49 (MQ=255)
                           ttATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGc                        <  2:53104/50‑1 (MQ=255)
                           ttATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACggg                         >  1:3215288/1‑49 (MQ=255)
                            tATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGc                        >  1:3278911/1‑49 (MQ=255)
                            tATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGc                        >  1:699618/1‑49 (MQ=255)
                            tATGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACGGGc                        >  2:5113632/1‑49 (MQ=255)
                             aTGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCaa                      >  2:2665321/1‑50 (MQ=255)
                             aTGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCaa                      >  1:4865855/1‑50 (MQ=255)
                             aTGCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACGGGCaa                      >  1:3723062/1‑50 (MQ=255)
                              tGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCaa                      >  1:4137761/1‑49 (MQ=255)
                               gCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACGGGCaa                      >  2:342087/1‑48 (MQ=255)
                               gCTGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACGGGCaa                      >  2:4024134/1‑48 (MQ=255)
                                cTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAg                   <  2:108579/50‑1 (MQ=255)
                                cTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAg                   >  1:5236931/1‑50 (MQ=255)
                                 tGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAg                   >  2:3806921/1‑49 (MQ=255)
                                 tGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAg                   >  2:5276803/1‑49 (MQ=255)
                                 tGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAg                   >  1:4112959/1‑49 (MQ=255)
                                 tGATATTTGCCGCCAGAATTTTCCAGCGCATCGACGATACGGGCAACAg                   >  2:1464113/1‑49 (MQ=255)
                                  gATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAgg                  >  2:3888980/1‑49 (MQ=255)
                                   atatatGCCGACAGGATTTTCCAGCGCAGCGACGATACGGGCAACAGGaa                <  2:2465509/50‑1 (MQ=38)
                                   atatTTGCCGCCAGGATTTTCCAGCGCATCGACGATAAGGGCAACAGGaa                <  1:5301655/50‑1 (MQ=255)
                                    tatTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAAt               <  1:1230306/50‑1 (MQ=255)
                                     atTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATc              >  1:2834332/1‑50 (MQ=255)
                                      tttGCCGCCAGGATATTCCAGCGCATCGACGATACGGGCAACAGGAATCg             <  2:5010477/50‑1 (MQ=255)
                                      tttGCCGCCAGAATTTTCCAGCGCATCGACGATACGGGCAACAGGAATc              <  1:591042/49‑1 (MQ=255)
                                        tGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGc            <  1:1888650/49‑1 (MQ=255)
                                         gccgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGcc           <  2:4087164/49‑1 (MQ=255)
                                          ccgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCg          <  1:5089597/49‑1 (MQ=255)
                                          ccgccAGAATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCg          <  2:108981/49‑1 (MQ=255)
                                           cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt         >  2:1986783/1‑49 (MQ=255)
                                           cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt         >  1:1136774/1‑49 (MQ=255)
                                           cgccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTc        <  2:3694428/50‑1 (MQ=255)
                                           cgccAGAATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGt         >  1:2340665/1‑49 (MQ=255)
                                            gccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCa       >  2:1469335/1‑50 (MQ=255)
                                            gccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCa       >  1:2951741/1‑50 (MQ=255)
                                             ccAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCAt      >  2:1370324/1‑50 (MQ=255)
                                              cAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCAtt     <  1:305141/50‑1 (MQ=255)
                                              cAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCAtt     <  1:2727649/50‑1 (MQ=255)
                                               aGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTc    >  1:109443/1‑50 (MQ=255)
                                                ggATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTcc   <  1:504961/50‑1 (MQ=255)
                                                ggATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTcc   <  2:1232999/50‑1 (MQ=255)
                                                ggATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTcc   >  2:1746276/1‑50 (MQ=255)
                                                 gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCg  >  1:1579500/1‑50 (MQ=255)
                                                 gATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCg  >  2:2970099/1‑50 (MQ=255)
                                                 |                                                 
ATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCG  >  NC_000913/1955644‑1955742

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