Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,008,406 G→A 6.5% E44K (GAA→AAA)  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,4060GA 6.5% 192.5 / 2.7 61E44K (GAA→AAA) yedEUPF0394 family sulphur transport domain‑containing inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (2/2);  total (30/31)
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.44e-01

GTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTT  >  NC_000913/2008260‑2008551
                                                                                                                                                  |                                                                                                                                                 
gTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGa                                                                                                                                                  <  2:515017/148‑1 (MQ=255)
                            gcaATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTa                                                                                                                                             >  2:237867/1‑125 (MQ=255)
                            gcaATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTa                                                                                                                                             <  1:237867/125‑1 (MQ=255)
                                ttCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTcc                                                                                                                        <  3:356049/142‑1 (MQ=255)
                                ttCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTcc                                                                                                                        >  4:356049/1‑142 (MQ=255)
                                               ttGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAAttt                                                                                                                                              <  4:280573/105‑1 (MQ=255)
                                               ttGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAAttt                                                                                                                                              >  3:280573/1‑105 (MQ=255)
                                                   ttAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGa                                                                                               >  4:538657/1‑148 (MQ=255)
                                                     aaaTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAag                                                                                             <  3:538657/148‑1 (MQ=255)
                                                       aTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTa                                                                                                                                             >  3:1138781/1‑98 (MQ=255)
                                                       aTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTa                                                                                                                                             <  4:1138781/98‑1 (MQ=255)
                                                       aTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGggc                                                                                                                                  <  1:631878/109‑1 (MQ=255)
                                                       aTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGggc                                                                                                                                  >  2:631878/1‑109 (MQ=255)
                                                       aTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAag                                                                                             >  4:136263/1‑146 (MQ=255)
                                                       aTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAag                                                                                             <  3:136263/146‑1 (MQ=255)
                                                            gggCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTgggg                                                                                      >  1:279108/1‑148 (MQ=255)
                                                                    aTCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTa                                                                              <  2:279108/148‑1 (MQ=255)
                                                                                  cgcgGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAagagg                                                                                          >  1:162237/1‑121 (MQ=255)
                                                                                  cgcgGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGt                                                                                          <  1:860822/122‑1 (MQ=255)
                                                                                  cgcgGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGt                                                                                          <  2:162237/122‑1 (MQ=255)
                                                                                  cgcgGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGt                                                                                          >  2:860822/1‑122 (MQ=255)
                                                                                  cgcgGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTg                                                                >  3:794128/1‑148 (MQ=255)
                                                                                               ctctctACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGt                                                                                          >  3:196171/1‑109 (MQ=255)
                                                                                               ctctctACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGt                                                                                          <  4:196171/109‑1 (MQ=255)
                                                                                               ctctctACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTa                                                                                   >  3:954786/1‑116 (MQ=255)
                                                                                               ctctctACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTa                                                                                   <  4:954786/116‑1 (MQ=255)
                                                                                                 ctctACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTAc                                                                                  <  2:135668/115‑1 (MQ=255)
                                                                                                 ctctACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTa                                                 >  1:539586/1‑148 (MQ=255)
                                                                                                 ctctACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTCCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTAc                                                                                  >  1:135668/1‑115 (MQ=255)
                                                                                                    tacttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAAcc                                              <  4:794128/148‑1 (MQ=255)
                                                                                                     acttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTc                                                                                                       <  4:349540/90‑1 (MQ=255)
                                                                                                     acttactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTc                                                                                                       >  3:349540/1‑90 (MQ=255)
                                                                                                        tactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATc                                                                  <  1:711924/124‑1 (MQ=255)
                                                                                                        tactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATc                                                                  >  2:711924/1‑124 (MQ=255)
                                                                                                        tactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTa                                                 <  2:852503/141‑1 (MQ=255)
                                                                                                        tactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTa                                                 >  1:852503/1‑141 (MQ=255)
                                                                                                        tactATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCa                                          >  1:906345/1‑148 (MQ=255)
                                                                                                          ctATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATc                                        >  2:1101109/1‑148 (MQ=255)
                                                                                                            aTTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTaa                                                                             <  3:203024/109‑1 (MQ=255)
                                                                                                            aTTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTaa                                                                             >  4:203024/1‑109 (MQ=255)
                                                                                                                        cTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGgatgatgat                          <  2:906345/148‑1 (MQ=255)
                                                                                                                             aCCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGt                                                                                     >  2:709282/1‑84 (MQ=255)
                                                                                                                             aCCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGt                                                                                     <  1:709282/84‑1 (MQ=255)
                                                                                                                              ccTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCa                                                    <  1:438426/116‑1 (MQ=255)
                                                                                                                              ccTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCa                                                    >  2:438426/1‑116 (MQ=255)
                                                                                                                               cTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGt                   >  4:524152/1‑148 (MQ=255)
                                                                                                                               cTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGt                   <  1:1101109/148‑1 (MQ=255)
                                                                                                                               cTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGt                   <  3:524152/148‑1 (MQ=255)
                                                                                                                                 tttGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTgg                                                                                        <  1:998954/77‑1 (MQ=255)
                                                                                                                                 tttGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTgg                                                                                        >  2:998954/1‑77 (MQ=255)
                                                                                                                                 tttGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTAt                 <  2:539586/148‑1 (MQ=255)
                                                                                                                                  ttGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATg                >  3:939475/1‑148 (MQ=255)
                                                                                                                                       cTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGAAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGa                                      >  2:644041/1‑121 (MQ=255)
                                                                                                                                       cTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGAAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTATAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGa                                      <  1:644041/121‑1 (MQ=255)
                                                                                                                                            aCGGGTAAATTCACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGAt                                                         <  4:384435/97‑1 (MQ=255)
                                                                                                                                            aCGGGTAAATTCACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCa                                                    >  3:520648/1‑102 (MQ=255)
                                                                                                                                            aCGGGTAAATTCACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCa                                                    <  4:520648/102‑1 (MQ=255)
                                                                                                                                            aCGGGTAAATTCACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAAGTATCCATCTGGAAGGAt                                                         >  3:384435/1‑97 (MQ=255)
                                                                                                                                             cGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTAc                                                                                  >  3:604751/1‑71 (MQ=255)
                                                                                                                                             cGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTAc                                                                                  <  4:604751/71‑1 (MQ=255)
                                                                                                                                                gTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCtt  <  4:939475/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GTTCAGAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGCTCCTGCAACTGTTCGGCGTCCATGCTGAAGAGTGGGGTTACTTTAAAATTATCCATCTGGAAGGATCGCCATTAACCCGCATCGACGGGATGATGATCCTCGGTATGTTTGGCGGCTGCTT  >  NC_000913/2008260‑2008551

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