Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,955,668 C→G 9.3% G447A (GGT→GCT)  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
*NC_0009131,955,6680CG 9.3% 134.9 / 6.7 43G447A (GGT→GCT) torZtrimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (16/22);  new base G (2/2);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

CTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCGCCGCAGACATTTCCGGCAGCACGCCGCCAACA  >  NC_000913/1955525‑1955811
                                                                                                                                               |                                                                                                                                               
ctgcGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTCCCGCCACCCGCCGACCAGATCATTTTGATATTCGGATACGGTTGTTCTTTACCGTTa                                                                                                                                            <  4:87994/149‑1 (MQ=255)
 tgcGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTAt                                                                                                                                           >  1:104217/1‑149 (MQ=255)
 tgcGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTAt                                                                                                                                           >  8:226423/1‑149 (MQ=255)
          gTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATAtt                                                                                                                                  >  2:283868/1‑149 (MQ=255)
           tAGCATTCAGACACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATAttt                                                                                                                                 <  5:128205/149‑1 (MQ=255)
           tAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATAAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATAttt                                                                                                                                 <  4:171640/149‑1 (MQ=255)
                   aGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAg                                                                                                                         >  5:319107/1‑149 (MQ=255)
                   aGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAg                                                                                                                         <  5:281820/149‑1 (MQ=255)
                     aaaCGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGa                                                                                                                       <  6:9186/149‑1 (MQ=255)
                        cgacgaTCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGAttt                                                                                                                    <  3:317975/149‑1 (MQ=255)
                             atcatcTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAg                                                                                                               <  7:226423/149‑1 (MQ=255)
                             atcatcTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAg                                                                                                               <  2:312175/149‑1 (MQ=255)
                              tcatcTCCGGTTTCTGCCATGCTTTATTCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGACCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCTCAAGGATTTTCCAgc                                                                                                              >  8:170924/1‑149 (MQ=255)
                              tcatcTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGa                                                                                                                                      <  8:407988/125‑1 (MQ=255)
                              tcatcTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGa                                                                                                                                      >  7:407988/1‑125 (MQ=255)
                               catcTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTAGCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgcg                                                                                                             >  1:283367/1‑149 (MQ=255)
                               catcTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGACCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTAGCGTTATGCTGATATTTGCCGCCAGGATTTTCCAgcg                                                                                                             >  5:22221/1‑149 (MQ=255)
                                     cGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATc                                                                                                         <  7:213099/147‑1 (MQ=255)
                                     cGGTTTATGACATGCTTTAATCAGACGATTAGTATCCTGGTGATGGTTAAAGTTACCTCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATcga                                                                                                       <  7:411537/149‑1 (MQ=255)
                                              ccATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTAGCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACggg                                                                                              <  2:283367/149‑1 (MQ=255)
                                              ccATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTAGCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACggg                                                                                              <  6:22221/149‑1 (MQ=255)
                                              ccATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACgg                                                                                               <  3:358269/148‑1 (MQ=255)
                                                 tGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCaa                                                                                           >  6:378250/1‑149 (MQ=255)
                                                   cTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATATTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAAct                                                                                         <  2:19219/149‑2 (MQ=255)
                                                     ttAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAg                                                                                        >  3:332688/1‑148 (MQ=255)
                                                         tCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATc                                                                                   <  1:283868/149‑1 (MQ=255)
                                                           aGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGc                                                                                 >  6:235065/1‑149 (MQ=255)
                                                               gATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCg                                                                                                                                               >  5:406884/1‑83 (MQ=255)
                                                               gATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCg                                                                                                                                               <  6:406884/83‑1 (MQ=255)
                                                                         cTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCca                                                                    <  1:19652/148‑1 (MQ=255)
                                                                          tGGTGATGGGTAAAGATACCGCCCCCCGCCCACCAGATCAATATGATAATCGGAAACGTATGATCTTTACCCTTATGCTGATATTAGCCGCCAGGATTTTCCAGCCCATCGACCATACCGGCAACAGGAATCGCCGTCAATCCGCcaa                                                                   >  8:427509/1‑147 (MQ=255)
                                                                             tGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCatcatc                                                               <  5:235065/149‑1 (MQ=255)
                                                                              gATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCcatcat                                                                >  1:423050/1‑147 (MQ=255)
                                                                                  ggTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCg                                                           <  5:378250/148‑1 (MQ=255)
                                                                                           aCCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTg                                                  <  6:319107/148‑1 (MQ=255)
                                                                                                   ccGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCg                                         >  5:297605/1‑149 (MQ=255)
                                                                                                           caGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATcgc                                 <  7:427509/149‑1 (MQ=255)
                                                                                                             gATCATTTTGATATTCGGATACGTTTGTTCTTTAACGTTATGCTGATATTTTCCGAAAGGATTTTACAGCGCATCGACGATAATGGAAACAGGAATCTCCGTCATTCAGCCATCATCAGCGGCATCTCTGGCATGTCAGGCAATcgcag                               >  8:318604/1‑147 (MQ=255)
                                                                                                                        tatTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCGCCGCAGACATTTcc                    >  8:267676/1‑149 (MQ=255)
                                                                                                                          ttCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGTATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCGCCGCAGACATTTCCgg                  <  7:170924/149‑1 (MQ=255)
                                                                                                                                    tttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCGCCGCAGACATTTCCGGCAGCAcgccg        >  4:297605/1‑149 (MQ=255)
                                                                                                                                     ttGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGACATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCGCCGCAGACATTTCCGGCAGCAcgccgc       >  3:116792/1‑149 (MQ=255)
                                                                                                                                             tACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCGCCGCAGACATTTCCGGCAGCACGCCGCCAaca  <  4:332688/146‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATGCTTTAATCAGACGATTAGTATCCTGGTGATGGGTAAAGTTACCGCCACCCGCCCACCAGATCATTTTGATATTCGGATACGTTTGTTCTTTACCGTTATGCTGATATTTGCCGCCAGGATTTTCCAGCGCATCGACGATACGGGCAACAGGAATCGCCGTCATTCCGCCATCATCAGCGGCCTCGCTGGCATGTCCGGCAATCGCCGCAGACATTTCCGGCAGCACGCCGCCAACA  >  NC_000913/1955525‑1955811

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