Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,040,0750CG40.5% 46.2 / 11.9 38L200V (CTG→GTG) yedYmembrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base C (4/18);  new base G (15/0);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 5.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCAC‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGG  >  NC_000913/2039937‑2040224
                                                                                                                                          |                                                                                                                                                      
gtgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGGAAAATTCtt                                                                                                                                            >  2:271602/1‑151 (MQ=255)
   ccGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATg                                                                                                                                         <  2:610489/151‑1 (MQ=255)
        gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATTTCg                                                                                                                                      >  2:378775/1‑149 (MQ=255)
             gTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGAtt                                                                                                                               <  1:499287/151‑1 (MQ=255)
             gTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGACCAGCAGGACCGCTTTATCGGCGGAGGGCTGGAATATCCTTATTTCGCAGGa                                                                                                                                 >  1:1057530/1‑149 (MQ=255)
              ttttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATACTTATGTCGAAgg                                                                                                                                  >  2:990529/1‑147 (MQ=255)
              ttttCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGGGAAATATCCCTTTTTCGAAgg                                                                                                                                  >  2:731528/1‑147 (MQ=255)
                          caAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTGGAATATTCCTATGTCGAAGGGTTGCGGCTcgac                                                                                                                     >  1:359503/1‑148 (MQ=255)
                            aaTTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCGCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGc                                                                                                                <  1:851165/151‑1 (MQ=255)
                              ttGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTGAAATATTCTTATTTCGGAGGGTTGGGTCTCGACGAAg                                                                                                                 >  1:747064/1‑148 (MQ=255)
                               tGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCaa                                                                                                              <  1:203153/150‑1 (MQ=255)
                                   ggCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATTCTTATGTCGAAGGAGTGCGTCTCGACGAAGCaa                                                                                                              >  2:1175419/1‑146 (MQ=255)
                                   ggCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCAt                                                                                                         <  1:348454/151‑1 (MQ=255)
                                             gAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCacac                                                                                               <  1:575208/151‑1 (MQ=255)
                                              aaCCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTcacaat                                                                                              <  1:737586/151‑3 (MQ=255)
                                              aaCCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGCTGAAAAATCCTTATGTCGAAGGATTTCGTCTCGACGAAGCAATGCATCCGCCcac                                                                                                 >  2:1103673/1‑148 (MQ=255)
                                                     cAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGGGAAAAATCCTTATGTCGAAGGGTTGGGGCTCGACGAAGCAATGCATCCGCTCACACt                                                                                              >  2:570073/1‑144 (MQ=255)
                                                         aaTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAg                                                                                   <  2:289608/151‑1 (MQ=255)
                                                            gCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGgt                                                                                 <  1:627368/150‑1 (MQ=255)
                                                             cGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTATGTCGAAGGATTGCGTCTCCGCGAAGCAAAGCATCCCCTCACACCGATGACCGTTgg                                                                                  >  1:1067228/1‑148 (MQ=255)
                                                                     gTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCa                                                                       <  2:891339/151‑1 (MQ=255)
                                                                       cGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAg                                                                     <  2:63447/151‑1 (MQ=255)
                                                                         cTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAAAGCATCCGCTCACACCGATGACCGGAGGGGTTTATGGCaa                                                                      >  1:223159/1‑148 (MQ=255)
                                                                              gAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATTCCTATTTTGACGGATTTGGTTTCGGCGAAGCAAAGCATCCGCTCCCACCGGTGACCGCAGGTGTTTTTGGCAAAGGGtt                                                                >  1:559536/1‑149 (MQ=255)
                                                                                aaCAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGCGAAATATCCTTTTTTCGAAGGGATGCGTCCCGACGAAGCAATGCATCCCCCCACACCGGTGACCGGAGGGGTTTTTGGCAAGGCGGTAccc                                                            >  2:708097/1‑150 (MQ=255)
                                                                                      ttATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAAACCTTTTGTCGGAGGGGTGCGGCTCGGCGGAGCAAAGCATCCGCTCCCACCGATGAACGCAGGGGTTTTTGGCAAGGCGGTACcccccc                                                        >  2:833666/1‑147 (MQ=255)
                                                                                             ccGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCCCACACTGATGACCGTAGGTGTTTATGGCAAGGCTTTACCGCCAC‑AAAATGgcg                                              <  2:869655/151‑1 (MQ=255)
                                                                                               ggAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAAACCTTTTGTCGGAGGGATGGGTCTCCGCGAAGCAATGCAATCGCCCACACCGATGACCGGAGGGGTTTATTGCAAGGCGGTAACCCCAC‑aaaaaagg                                               >  2:355216/1‑144 (MQ=255)
                                                                                               gaaaCAGATGCCAGGCAAGCAGACCCGTTTAACCGGGGGTGGTATGAAATACCTTTAGGGCAAAGAATTGCATCGCGACAAAGCAAAGCATCCCCACACACCGATGACCGTAGGTGTTTGTGGCAAGGCGTTACCGCCAC‑AAAAt                                                  <  2:1053600/143‑1 (MQ=255)
                                                                                                     gATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTTGGGTGAAAAATCCTTATGGCCTAGGGGTGGGGCCCGCCGAAAGAGAGAATTCGCTCCCACAGATGACCGAAGGTGGTTTTGGCCAGGTGTTAACGCCCCAAAAATGGCGGCCCggg                                       >  2:412537/1‑150 (MQ=255)
                                                                                                     gATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAAAATCCTTATGTTGAAGGGTTTCGTCTCCACGAAACAAAGCATCCGCCCACACTTATGGCCGGAGGGGGTTATTGGAAGGGGTTACCCCCAC‑aaaaaaggggggcgg                                        >  1:511570/1‑138 (MQ=255)
                                                                                                            ggccagAAGGCCCGTTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGGGTCTCGACGAAGCAATGCACCCGCTCACACTGATGCCGGTAGGTGTTTATGGAAAGGCTTTACCGCCAC‑AAAATGGCGCGCCGGTGCGACTGa                               <  2:583145/151‑1 (MQ=255)
                                                                                                                    ggACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCAC‑AAAATGGCGCGCCGGTGCGACTGATtgtg                          <  1:91547/148‑1 (MQ=255)
                                                                                                                                 gggggTGGGCTAAAATATCCTTAGGTCGAAGGATTGCGTCTCGACAAAGCAATGCATCCGCTCACCCTGATGACCGTAGGTGTTTAGGGAAAGGTGTTACCGCCAC‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAatat              <  2:327232/144‑1 (MQ=255)
                                                                                                                                 ggcggTGGTCAAAAATACCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCCTCCCCTCACAATGATGCCCGTCGGTGTTTGTGGCAAGGCGTTACCGCCAC‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAat                <  2:292259/145‑1 (MQ=255)
                                                                                                                                    ggggggCTGAAAAATCCTTATGTCGAAGTTTTGTGTCTGGGGAAGGCAAACCATCCTCTCACACGGATGCCCGTAGATGTTTATGGCAAGGCTTTTCCGCCAC‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCtt         <  1:994968/146‑1 (MQ=255)
                                                                                                                                     gTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGAAGGGCGTTACCGCCAC‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCttt        <  2:846162/149‑1 (MQ=255)
                                                                                                                                        ggCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCAC‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAgg   >  2:87009/1‑151 (MQ=255)
                                                                                                                                          cTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCAC‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAggg  >  1:331230/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                      
GTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCAC‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGG  >  NC_000913/2039937‑2040224

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑