Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 264,098 G→T L50L (CTC→CTA)  ykfH ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA264,0980GT95.9% 140.7 / NA 50L50L (CTC→CTA) ykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (22/25);  minor base C (2/0);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTT  >  NC_000913_3_pae_tpiA/263954‑264231
                                                                                                                                                |                                                                                                                                     
ggCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGGCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGa                                                                                                                                   >  2:537537/1‑149 (MQ=255)
     ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                              >  2:212601/1‑149 (MQ=255)
      cgggCGCAGTGTAACCCCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATAGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                                                                                             <  1:97508/147‑1 (MQ=255)
       aGGCGCAGTGTAATCTCCCGGTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGc                                                                                                                            <  2:386428/149‑1 (MQ=255)
         gcgcAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTGGGCGTATACCGCTAGGACTCCGGCgg                                                                                                                          >  1:58876/1‑148 (MQ=255)
              gtgtATTCCCCCGCTGCAGGGACCCGGTGGGGATATAAATAACACCTCTTCACTATTGTCAGTTAAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATcc                                                                                                                     <  1:297389/149‑1 (MQ=255)
               tgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTACGACTCCGGCGATATCCg                                                                                                                    >  2:17015/1‑149 (MQ=255)
                    tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACg                                                                                                               >  2:443803/1‑149 (MQ=255)
                     ctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                                                                                              >  2:577135/1‑149 (MQ=255)
                     ctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTCATGGTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGt                                                                                                              >  1:286233/1‑149 (MQ=255)
                      tcCCGCTGCAGGCCACTGGTAGGGTTGCTCATTCCAGTTCAGCACTAATGTCAGCTCATGATGGTGTTCGTCTGGGAGGTATAGGGCCCATTGCGGTCCTGGCTGCGTTCGGCGTATACCGCCAGGGCTCCGGCGGTATCAGGGAGGTc                                                                                                             >  2:137238/1‑149 (MQ=255)
                            tGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGCTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGcc                                                                                                       >  1:182233/1‑149 (MQ=255)
                            tGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGcc                                                                                                       >  1:162767/1‑149 (MQ=255)
                                 gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTcc                                                                                                            >  1:253011/1‑139 (MQ=255)
                                 gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTcc                                                                                                            <  2:253011/139‑1 (MQ=255)
                                   cccTGGTGGGGTTGCTCATTACGGTTCTCCGCTATTGTCAGTTCAGGGGGATGCTGATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                <  2:523520/149‑1 (MQ=255)
                                    tgtGGTAGGTTGGCTCATGACTGTTCTCGCATATGGTAAGTTCAGGGTGATGCTGAGCAGGCAGGTATAGGGCCCAGTGGGGGCCTGGCTGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                               <  2:583524/147‑1 (MQ=255)
                                          aGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGGCCTGGCGGCGTTCGACGTATACCGCTAGGGCTCCGGCGATAACAGGAACGCCCCTGCCGGGGTGAtatcgga                                                                                         >  2:245181/1‑143 (MQ=255)
                                             gTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGc                                                                                      <  1:178806/149‑1 (MQ=255)
                                               ggCTCATTACAGTTCTCCATTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTa                                                                                    <  2:249033/148‑1 (MQ=255)
                                                 cTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTAcc                                                                                  <  2:162767/149‑1 (MQ=255)
                                                          gTTCTCCACTATTGTCGGTTAAGGGTGTTGCTCAGCAGGCAGGTATAGGGCCCGTTGCGGTCCTGGCGGCGTTCGGCGGATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                                         <  2:567022/149‑1 (MQ=21)
                                                              tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGAAGCTACCGTGCCCCTGGTAt                                                                     >  1:586922/1‑149 (MQ=255)
                                                               ccACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAtt                                                                    >  1:321023/1‑149 (MQ=255)
                                                                cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAttt                                                                   >  1:619122/1‑149 (MQ=255)
                                                                  cTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGc                                                                 <  1:293606/149‑1 (MQ=255)
                                                                  cTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGc                                                                 <  2:182233/149‑1 (MQ=255)
                                                                    atagTCAGTTCAGGGTGATGCCCATCAGGCACGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                               <  1:452469/146‑1 (MQ=255)
                                                                       gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAACGCCACACGCTACCGTGCCACTGGCATTTGCCCGTg                                                            >  2:527425/1‑149 (MQ=255)
                                                                        tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGc                                                           >  2:629219/1‑149 (MQ=255)
                                                                           gTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCCAGGACTCCGGCGATATACGACACGCCCCTGCCCGGCTAATGACAGAAGCTACCGTGCCACTGGAATTTGCCGGTGCact                                                        >  2:163120/1‑147 (MQ=255)
                                                                              cAGGGGGTTGCTCATAAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGGCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGc                                                     <  2:149919/149‑1 (MQ=25)
                                                                                 ggTGATGCTCATCAGGCAGGTATAGGGCCCTTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCTCTGGTATTTGCCGGTGCAGTAGCGaa                                                  <  1:136179/149‑1 (MQ=255)
                                                                                    gATGCTCATCAGGCAGGTATAGGGCCCATTGGGGTCATGGCGGCGTTCGGCGTATAGCGCTAGGACTCCGGCGGTATCCGGAACGTCCCTGCCGGTGTCATGACAGACGCTACCGTGCCCCCGGTCTTAGCCGGTGCCGTGGCGcaaca                                               >  2:606191/1‑147 (MQ=21)
                                                                                    gATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAgt                                               >  2:578934/1‑148 (MQ=255)
                                                                                     aTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGAt                                              >  1:384722/1‑149 (MQ=255)
                                                                                        cccttcaGGTATGTATTGGGCCTATTGCGGCTCTGGCGGTGTTCGGCGTATATCGCTATGCTGCCGGCGATATCCGGAACGTCCGTGCCGGTGTATTGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                           <  2:216472/145‑1 (MQ=18)
                                                                                         tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCTTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                          <  1:325710/149‑1 (MQ=255)
                                                                                         tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                          <  1:212601/149‑1 (MQ=255)
                                                                                           atcaGGCAGGTATAGGGCCCCTTGCGGTCCTGGCGGGCTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                                        <  2:266299/149‑1 (MQ=255)
                                                                                               ggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                    <  2:122426/149‑1 (MQ=255)
                                                                                               ggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCa                                    <  1:441459/149‑1 (MQ=255)
                                                                                                   ggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGAGTATACAGGTAGGCGTCGGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAt                                <  2:619122/149‑1 (MQ=18)
                                                                                                       taGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTg                            <  2:306739/149‑1 (MQ=255)
                                                                                                            cccATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGATGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGt                       >  1:13143/1‑149 (MQ=255)
                                                                                                                   cGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCa                                                                            <  1:358880/89‑1 (MQ=255)
                                                                                                                   cGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCa                                                                            >  2:358880/1‑89 (MQ=255)
                                                                                                                      tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGGGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCCCTGGCCAGATATGGAGTGCTCAGGATGCTGGCGGTATATCGTCAt             >  1:442528/1‑149 (MQ=17)
                                                                                                                          ggcggcGGTCGGCGTATACGGGTGGGATTCCGGCGGTATCCGGAACGTCCCTGCCGGTGGAATGACTGACGCTACCGGGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGGCTCAGGATGTTGGCGGTATATCGTCATTGcc         <  2:76748/149‑1 (MQ=14)
                                                                                                                                 ttCGGCGTACACCGCTAGGTGTGCGCCGATATCCGGAACGCCCTTGCCGGTGTAATGAGAGTCGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGtt  <  2:586922/149‑1 (MQ=14)
                                                                                                                                             cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTc               >  2:220546/1‑124 (MQ=255)
                                                                                                                                                |                                                                                                                                     
GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTT  >  NC_000913_3_pae_tpiA/263954‑264231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑