Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,008,558 A→T 8.5% A94A (GCA→GCT)  yedE → UPF0394 family sulphur transport domain‑containing inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,008,5580AT 8.5% 152.5 / 2.8 47A94A (GCA→GCT) yedEUPF0394 family sulphur transport domain‑containing inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (23/20);  new base T (2/2);  total (25/22)
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.96e-01

CCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCA  >  NC_000913/2008414‑2008704
                                                                                                                                                |                                                                                                                                                  
ccGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAgcgc                                                                                                                                                >  2:274030/1‑149 (MQ=255)
 cGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCt                                                                                                                                               <  2:141252/149‑1 (MQ=255)
   ttGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCtgt                                                                                                                                             >  3:110687/1‑149 (MQ=255)
     ggggCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTCGGCGGCTGCTTTGCCGCAGCGCTGTgg                                                                                                                                           <  3:379365/149‑1 (MQ=255)
              aGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATTGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAAt                                                                                                                                  <  3:9792/149‑1 (MQ=255)
              aGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCACTGTGGGCCAACAAt                                                                                                                                  <  3:149676/149‑1 (MQ=255)
              aGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTAAGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAAt                                                                                                                                  <  4:327875/149‑1 (MQ=255)
                   cTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCTGCGCTGTGGGCCAACAATGTCaa                                                                                                                             >  4:69264/1‑149 (MQ=255)
                   cTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCTGCGCTGTGGGCCAACAATGTCaa                                                                                                                             >  2:493461/1‑149 (MQ=255)
                    tGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCaaa                                                                                                                            >  1:436757/1‑149 (MQ=255)
                       aaCTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGAATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTg                                                                                                                         >  2:477649/1‑149 (MQ=255)
                          tGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGa                                                                                                                      <  1:86439/149‑1 (MQ=255)
                            ttCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGAGGGCCAACAATGTCAAACTGCGAAt                                                                                                                    >  3:362936/1‑149 (MQ=255)
                                           gAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCTGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTAt                                                                                                     <  1:493461/149‑1 (MQ=255)
                                           gAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCTGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTAt                                                                                                     <  3:69264/149‑1 (MQ=255)
                                            aaGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATc                                                                                                    <  3:309682/149‑1 (MQ=255)
                                            aaGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCA‑TAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATcc                                                                                                   >  3:435771/1‑149 (MQ=255)
                                                gTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCa                                                                                                <  2:51193/149‑1 (MQ=255)
                                                         cTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAgg                                                                                       >  4:24565/1‑149 (MQ=255)
                                                          tttAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGc                                                                                      <  4:110528/149‑1 (MQ=255)
                                                          tttAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGc                                                                                      >  1:352830/1‑149 (MQ=255)
                                                             aaaaTTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCcat                                                                                   >  3:331019/1‑149 (MQ=255)
                                                               aaTTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCcatca                                                                                 <  4:310100/149‑1 (MQ=255)
                                                               aaTTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGATGCTTTGCCGCAGCGCTGTGGGCCAACCATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCCGGCcatca                                                                                 >  4:92146/1‑149 (MQ=255)
                                                                  tATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAATCTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTg                                                                              <  2:495915/149‑1 (MQ=255)
                                                                   aTCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTgg                                                                             >  3:70848/1‑149 (MQ=255)
                                                                      cATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTggcgg                                                                          <  3:161300/149‑1 (MQ=255)
                                                                      cATCGGGAAGGATCGCCATTAACCCGAATCGCCGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTggcgg                                                                          >  2:373060/1‑149 (MQ=255)
                                                                       aTCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTggcggc                                                                         >  3:526639/1‑149 (MQ=255)
                                                                           ggAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCTGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTggcgg                                                                          >  4:492456/1‑144 (MQ=255)
                                                                           ggAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCcaa                                                                                                                                      <  4:314672/84‑1 (MQ=255)
                                                                           ggAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCcaa                                                                                                                                      >  3:314672/1‑84 (MQ=255)
                                                                           ggAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTggcgg                                                                          <  3:492456/144‑1 (MQ=255)
                                                                            gAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTAt                                                                    >  1:126040/1‑149 (MQ=255)
                                                                                       aTTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTg                                                         <  2:436757/149‑1 (MQ=255)
                                                                                                         gatgatgatCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATgg                                       >  1:494052/1‑149 (MQ=255)
                                                                                                         gatgatgatCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTAACGCCGGTTTTGGCGCGCGTCTGGCAATgg                                       >  3:302677/1‑149 (MQ=255)
                                                                                                          atgatgatCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGTATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATggg                                      <  1:477649/149‑1 (MQ=255)
                                                                                                                 tCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAAc                               >  1:227262/1‑149 (MQ=255)
                                                                                                                  ccTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAAcc                              >  3:545399/1‑149 (MQ=255)
                                                                                                                        tATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCCTTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTggcgg                        >  4:257320/1‑149 (MQ=255)
                                                                                                                              tGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGTGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCt                  <  4:445771/149‑1 (MQ=255)
                                                                                                                               ggcggcTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCtt                 <  3:200192/149‑1 (MQ=255)
                                                                                                                                      gctTTGCCGCAGCGCTGTGGGCTAACAATGTCAAAATGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGt          >  2:202047/1‑149 (MQ=255)
                                                                                                                                        tttGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTAt        >  1:74922/1‑149 (MQ=255)
                                                                                                                                         ttGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTAtt       <  2:227262/149‑1 (MQ=255)
                                                                                                                                              gCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCa  <  4:435771/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTTTGCCGCAGCGCTGTGGGCCAACAATGTCAAACTGCGAATGCCGCGCAGCCGTATCCGCATTATGCAGGCCATCATTGGCGGCATTATCGCCGGTTTTGGCGCGCGTCTGGCAATGGGCTGTAACCTGGCGGCGTTCTTTACCGGTATTCCTCA  >  NC_000913/2008414‑2008704

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