Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,790,475 T→C Y257C (TAC→TGC)  ydiU ← UPF0061 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,790,4750TC98.1% 174.8 / ‑6.8 52Y257C (TAC→TGC) ydiUUPF0061 family protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (32/19);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCA  >  NC_000913/1790329‑1790621
                                                                                                                                                  |                                                                                                                                                  
caacaaATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATc                                                                                                                                                >  2:107389/1‑151 (MQ=255)
      aTGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGt                                                                                                                                          <  1:617255/151‑1 (MQ=255)
           gACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGcc                                                                                                                                     <  1:344479/151‑1 (MQ=255)
              aaTGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGcccc                                                                                                                                   >  2:170640/1‑150 (MQ=255)
                tGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCcagc                                                                                                                                >  2:158448/1‑151 (MQ=255)
                     gcgcCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGa                                                                                                                           <  1:561720/151‑1 (MQ=255)
                      cgcCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGAc                                                                                                                          >  1:149401/1‑151 (MQ=255)
                       gcCAGACGCTGTAAATTCCACAACGCGACGGAAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGAc                                                                                                                          >  2:658885/1‑150 (MQ=255)
                        ccAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACAt                                                                                                                        >  2:50162/1‑151 (MQ=255)
                        ccAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACAt                                                                                                                        >  1:38967/1‑151 (MQ=255)
                                   aaaTTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCGATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTAt                                                                                                             >  2:31004/1‑151 (MQ=255)
                                      ttCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATtc                                                                                                           >  1:407562/1‑150 (MQ=255)
                                           caACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCAcc                                                                                                     >  2:228455/1‑151 (MQ=255)
                                            aaCGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCAccc                                                                                                    >  2:784828/1‑151 (MQ=255)
                                                cgACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATg                                                                                                <  1:799675/151‑1 (MQ=255)
                                                 gACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATg                                                                                                >  2:826581/1‑150 (MQ=255)
                                                    ggCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCa                                                                                            <  2:559874/151‑1 (MQ=255)
                                                      cAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCaaa                                                                                          <  1:585844/151‑1 (MQ=255)
                                                       aGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAg                                                                                         <  1:728461/151‑1 (MQ=255)
                                                         gTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGcc                                                                                       >  2:790245/1‑151 (MQ=255)
                                                            gATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGAc                                                                                    >  2:211610/1‑151 (MQ=255)
                                                            gATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGAc                                                                                    >  2:779340/1‑151 (MQ=255)
                                                            gATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGAc                                                                                    >  1:226660/1‑151 (MQ=255)
                                                                  cAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCt                                                                               <  1:89936/150‑1 (MQ=255)
                                                                  cAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTg                                                                              <  2:168782/151‑1 (MQ=255)
                                                                     aGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAATCGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCa                                                                           <  1:45044/151‑1 (MQ=255)
                                                                           aaCGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTggg                                                                      >  2:572298/1‑150 (MQ=255)
                                                                             agCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCa                                                                    >  1:232757/2‑150 (MQ=255)
                                                                                ccTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAAt                                                                  >  1:432171/1‑149 (MQ=255)
                                                                                ccTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTa                                                                <  1:826523/151‑1 (MQ=255)
                                                                                       gatCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAacgacg                                                          <  1:524389/150‑1 (MQ=255)
                                                                                       gatCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGAcgc                                                         <  1:696528/151‑1 (MQ=255)
                                                                                        atCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGAcgcg                                                        >  2:613266/1‑151 (MQ=255)
                                                                                          ccGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCgg                                                       <  1:431723/150‑1 (MQ=255)
                                                                                          ccGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGt                                                      >  2:121393/1‑151 (MQ=255)
                                                                                              gTGATTGCTAATAAAACCGGGTTCGTAATCATCAAGAAAACCAACCGCCCCGTAATCATGCGACCGCACCAGCAGCGACATTATGTCGGTATTCATCACCCCATGGGCAACGCCGACCGGCTGCCATGGGGCAATAAACGAGGCGGGACg                                                   >  2:668712/1‑150 (MQ=255)
                                                                                              gTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGt                                                  >  1:607859/1‑151 (MQ=255)
                                                                                              gTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGt                                                  >  2:39714/1‑151 (MQ=255)
                                                                                                        aTAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACAt                                        >  2:155577/1‑151 (MQ=255)
                                                                                                            aaCCGGGTTCGTAATCATCAAGAAAACCAAACGCCCCGCAATCACGCGTCAGTCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACCCGGTACGTGCGACAACATCGCt                                    >  1:557750/1‑151 (MQ=255)
                                                                                                                     cGTAATCATCAAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCaga                            <  2:198671/150‑1 (MQ=255)
                                                                                                                           catcaAGAAAACCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTa                     <  2:231499/151‑1 (MQ=255)
                                                                                                                               aaGAAAACCAAACGGCCCGCAATCAAGTGACAGCCCCAGCAGCGACATGATGACGGTATTCATCACCCCATGAGCAAAGCCGACCGTCAGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTCAACCAGAGACGGTAttt                  >  1:753150/1‑150 (MQ=255)
                                                                                                                                    aaCCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCtc            >  1:159863/1‑151 (MQ=255)
                                                                                                                                     aCCAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCtc            >  2:366754/1‑150 (MQ=255)
                                                                                                                                      ccAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCtcat          >  2:571442/1‑151 (MQ=255)
                                                                                                                                      ccAAACGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCtcat          <  2:490968/151‑1 (MQ=255)
                                                                                                                                          aCGGCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCatcatc      <  2:853660/151‑1 (MQ=255)
                                                                                                                                            ggCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTg    <  2:449043/151‑1 (MQ=255)
                                                                                                                                             gCCCGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGc   >  2:364208/1‑151 (MQ=255)
                                                                                                                                               ccGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCa  >  2:672187/1‑150 (MQ=255)
                                                                                                                                               ccGCAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCa  >  1:340018/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAACAAATGGCGACAATGTCTGCGCCAGACGCTGTAAATTCCACAACGCGACGGCAGGTTGATTATCAAAGCTGTAACGCCCTTGATGATCCGAGTGATTACAAATAAAACCGGGTTCGTAATCATCAAGAAAACCAAACGGCCCGTAATCAAGCGTCAGCCCCAGCAGCGACATGTTGTCGGTATTCATCACCCCATGAGCAAAGCCGACCGTCTGCCATTGGGCAATTAACGACGCGGTACGTGCGACAACATCGCTAAACCAGAGACGGTATTTGTCCTCATCATCTGCA  >  NC_000913/1790329‑1790621

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