Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,770,856 (T)6→5 coding (242/357 nt) ydiL → DUF1870 domain‑containing protein YdiL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,770,8510T.97.7% 171.1 / ‑1.8 43coding (237/357 nt)ydiLDUF1870 domain‑containing protein YdiL
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (17/25);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CGTCAGTGGGAAAATGGCAAGTGTGCCATTCCTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAA  >  NC_000913/1770714‑1770985
                                                                                                                                         |                                                                                                                                      
cGTCAGTGGGAAAATGGCAAGTGTGCCATTCCTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGttt                                                                                                                                      >  2:100760/1‑140 (MQ=255)
  tCAGTGGGAAAATGGCAAGTGTGCCATTCCTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGttttt                                                                                                                                    >  1:548022/1‑140 (MQ=255)
     gTGGGAAAATGGCAAGTGTGCCATTCCTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTccc                                                                                                                                 >  2:663435/1‑140 (MQ=255)
                        gCCATTCCTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCaa                                                                                                             <  1:659671/140‑1 (MQ=255)
                        gCCATTCCTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCaa                                                                                                             <  1:19322/140‑1 (MQ=255)
                         ccATTCCTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAAc                                                                                                            <  2:19652/140‑1 (MQ=255)
                         ccATTCCTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAAc                                                                                                            >  1:604746/1‑140 (MQ=255)
                               cTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCt                                                                                                      >  1:495017/1‑140 (MQ=255)
                               cTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCt                                                                                                      >  1:572759/1‑140 (MQ=255)
                                 gATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTAc                                                                                                    >  2:40598/1‑140 (MQ=255)
                                  aTCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTAcc                                                                                                   >  2:21091/1‑140 (MQ=255)
                                      tgtAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGa                                                                                               >  1:538378/1‑140 (MQ=255)
                                      tgtAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGa                                                                                               >  1:646051/1‑140 (MQ=255)
                                        tAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATg                                                                                             >  1:104207/1‑140 (MQ=255)
                                         aGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATgg                                                                                            >  1:405255/1‑140 (MQ=255)
                                         aGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATgg                                                                                            <  2:485002/140‑1 (MQ=255)
                                                 aaCTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTa                                                                                    <  2:604746/140‑1 (MQ=255)
                                                   cTGTTGGCTATGCGTCAGCAAATAAAAAAACCTCTTCATACCATCATACAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTAtt                                                                                  <  2:284995/140‑1 (MQ=255)
                                                   cTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTAtt                                                                                  <  2:439695/140‑1 (MQ=255)
                                                   cTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTAtt                                                                                  <  2:572759/140‑1 (MQ=255)
                                                   cTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTAtt                                                                                  <  1:100760/140‑1 (MQ=255)
                                                      ttGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGAc                                                                               <  2:466761/140‑1 (MQ=255)
                                                              gCGTCAGCAAAGAAAAAAACATCTTCATGGCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGCCAATTTAATTGACTAGAAGAt                                                                       >  2:77104/1‑140 (MQ=255)
                                                                  cAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTAt                                                                   <  2:538378/140‑1 (MQ=255)
                                                                   aGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTTTACCCTGGTGGCAACTTTATTAACTGGAAGATCTATc                                                                  >  1:319800/1‑140 (MQ=255)
                                                                   aGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATc                                                                  >  2:641258/1‑140 (MQ=255)
                                                                     cAAAGAAAAAAACATCTTCATGTCACCATCGAGAAAACAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCaa                                                                <  1:183490/140‑1 (MQ=255)
                                                                          aaaaaaaCATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCGTGATGGCGATTTTATTGACTGAAAGATCTATCAATCGGt                                                           >  2:183490/1‑140 (MQ=255)
                                                                          aaaaaaaCATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGt                                                           <  2:43407/140‑1 (MQ=255)
                                                                               aaCATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAgctg                                                      <  1:166073/140‑1 (MQ=255)
                                                                               aaCATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAgctg                                                      <  2:405255/140‑1 (MQ=255)
                                                                                aCATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAgctgc                                                     <  1:610761/140‑1 (MQ=255)
                                                                                 cATCTTCATGCCATCATAGAGAAAATAAATAACCGCCTCGGCAACAACCCCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAgctgct                                                    <  1:470098/140‑1 (MQ=255)
                                                                                               catAGAGAAAATAAATCGCCGCTTCGGCAACAACACCATCCGTTTTTCTCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGcaca                                      <  2:319800/140‑1 (MQ=255)
                                                                                               catAGAGAAAATAAATAACCGCATCGGCAACAACATCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCCGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGcaca                                      <  2:495017/140‑1 (MQ=255)
                                                                                               catAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGcaca                                      >  1:490602/1‑140 (MQ=255)
                                                                                               catAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGcaca                                      <  2:116519/140‑1 (MQ=255)
                                                                                                 tAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACg                                    >  1:116310/1‑140 (MQ=255)
                                                                                                      aaaaTAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATtata                                            <  2:640640/127‑1 (MQ=255)
                                                                                                      aaaaTAAATAACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATtata                                            >  1:640640/1‑127 (MQ=255)
                                                                                                        aataaataACCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGa                             <  1:663435/140‑1 (MQ=255)
                                                                                                            aataaCCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCgg                         <  1:500712/140‑1 (MQ=255)
                                                                                                                aCCGCATCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTGGCTGCTGAATTATACGCACACGATCTGGAGCGGCttt                     >  1:260576/1‑140 (MQ=255)
                                                                                                                     aTCGGCAACAACACCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTaa                <  2:548022/140‑1 (MQ=255)
                                                                                                                             caacaCCATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCca        <  2:283153/140‑1 (MQ=255)
                                                                                                                                  ccATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGaa   <  2:333180/140‑1 (MQ=255)
                                                                                                                                   cATGCGTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACAAGATCTGGAGCGGCTTTGTTAATTTTTCCACAGaaa  >  2:271029/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                      
CGTCAGTGGGAAAATGGCAAGTGTGCCATTCCTGATCGTGTAGTCGAGCAACTGTTGGCTATGCGTCAGCAAAGAAAAAAACATCTTCATGCCATCATAGAGAAAATAAATAACCGCATCGGCAACAACACCATGCGTTTTTTCCCCGACTTAACCGCATTTCAACGAGTCTACCCTGATGGCAATTTTATTGACTGGAAGATCTATCAATCGGTAGCTGCTGAATTATATGCACACGATCTGGAGCGGCTTTGTTAATTTTTCCACAGAAA  >  NC_000913/1770714‑1770985

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