Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,101,027 G→T intergenic (+95/‑250) ycdB → / → phoH hypothetical protein/conserved protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,101,0270GT98.5% 254.1 / ‑5.6 68intergenic (+95/‑250)ycdB/phoHhypothetical protein/conserved protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base T (33/33);  minor base C (1/0);  total (35/33)
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.98e-01

TTTGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACAC  >  REL606/1100954‑1101095
                                                                         |                                                                    
tttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTAt                                                                    >  2:111549/1‑76 (MQ=255)
tttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTAt                                                                    <  1:1167056/76‑1 (MQ=255)
 ttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGtata                                                                   <  1:925718/76‑1 (MQ=255)
 ttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGtata                                                                   >  1:1582629/1‑76 (MQ=255)
 ttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGtata                                                                   <  2:1487842/76‑1 (MQ=255)
   gtgCGAAGATGCCTAATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAg                                                                 <  2:324243/76‑1 (MQ=25)
    tgCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAgg                                                                >  2:841247/1‑76 (MQ=255)
    tgCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAgg                                                                <  1:267874/76‑1 (MQ=255)
    tgCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGAGTGGGTTTGGTAGGCTGTATAAgg                                                                >  2:841243/1‑76 (MQ=21)
     gcgcAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGc                                                               <  2:1230991/76‑1 (MQ=255)
        cAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGtt                                                            <  1:1242773/76‑1 (MQ=255)
        cAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGtt                                                            <  2:510013/76‑1 (MQ=255)
        cAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGtt                                                            <  2:633944/76‑1 (MQ=255)
         aGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTc                                                           >  2:1544932/1‑76 (MQ=255)
         aGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTc                                                           >  2:663421/1‑76 (MQ=255)
          gATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCa                                                          <  1:411851/76‑1 (MQ=255)
          gATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCa                                                          <  2:1367165/76‑1 (MQ=255)
           aTGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAc                                                         <  1:29269/76‑1 (MQ=255)
            tGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcg                                                        >  2:552964/1‑76 (MQ=255)
             gCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgc                                                       <  1:441312/76‑1 (MQ=255)
             gCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgc                                                       <  1:485268/76‑1 (MQ=255)
             gCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgc                                                       <  1:1503357/76‑1 (MQ=255)
             gCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgc                                                       <  1:1409588/76‑1 (MQ=255)
              ccTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgcc                                                      >  2:1446213/1‑76 (MQ=255)
               cTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgccg                                                     >  2:943010/1‑76 (MQ=255)
               cTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgccg                                                     >  1:391392/1‑76 (MQ=255)
               cTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCAcgccg                                                     >  1:49727/1‑76 (MQ=255)
                tGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgccgc                                                    >  2:57718/1‑76 (MQ=255)
                tGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgccgc                                                    >  2:1628073/1‑76 (MQ=255)
                 gATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCa                                                   >  2:1551395/1‑76 (MQ=255)
                    gCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATcc                                                <  2:1394964/76‑1 (MQ=255)
                        cGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCg                                            >  2:13186/1‑76 (MQ=255)
                         gCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGa                                           <  1:1543430/76‑1 (MQ=255)
                          cTTGCGCGTCTTATCATGCCTACAATCGGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGAt                                          >  1:705460/1‑76 (MQ=25)
                               gcgTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGc                                     <  2:1333391/76‑1 (MQ=255)
                                cgTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCa                                    >  1:804929/1‑76 (MQ=255)
                                cgTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGCATAAGGCGTTCACGCCGCATCCGGCGATCGTGcc                                    >  1:817302/1‑75 (MQ=255)
                                  tCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCAct                                  >  1:866996/1‑76 (MQ=255)
                                  tCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCAct                                  >  2:891367/1‑76 (MQ=255)
                                   cTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActg                                 >  2:1614263/1‑76 (MQ=255)
                                   cTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActg                                 >  2:1287886/1‑76 (MQ=255)
                                   cTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActg                                 >  2:1619019/1‑76 (MQ=255)
                                    ttATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActga                                >  1:931031/1‑76 (MQ=255)
                                     tATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActgat                               >  2:1131657/1‑76 (MQ=255)
                                      aTCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActgatg                              >  2:446449/1‑76 (MQ=255)
                                            ccTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActgatgcctgat                        <  2:605604/76‑1 (MQ=255)
                                              tACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActgatgcctgatgc                      <  2:144274/76‑1 (MQ=255)
                                                cAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCaa                    >  1:411584/1‑76 (MQ=255)
                                                cAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCaa                    >  2:431051/1‑76 (MQ=255)
                                                 aaTCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCaaa                   <  2:1106671/76‑1 (MQ=255)
                                                 aaTCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCaaa                   <  2:581703/76‑1 (MQ=255)
                                                  aTCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAAt                  <  2:6267/76‑1 (MQ=255)
                                                  aTCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAAt                  <  2:1484499/76‑1 (MQ=255)
                                                   tCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATc                 <  1:634225/76‑1 (MQ=255)
                                                     aGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCct               >  2:128064/1‑76 (MQ=255)
                                                     aGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCct               >  1:1436209/1‑76 (MQ=255)
                                                       tGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgc             <  1:1612117/76‑1 (MQ=255)
                                                       tGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgc             <  2:1153682/76‑1 (MQ=255)
                                                        gCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgct            <  2:1567498/76‑1 (MQ=255)
                                                         cGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgctg           <  1:311236/76‑1 (MQ=255)
                                                         cGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgctg           <  1:907362/76‑1 (MQ=255)
                                                           ggTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGaa         <  2:828216/76‑1 (MQ=255)
                                                            gTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGaaa        >  2:544718/1‑76 (MQ=255)
                                                            gTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGaaa        >  2:1582348/1‑76 (MQ=255)
                                                             tttGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAg       >  2:1169529/1‑76 (MQ=255)
                                                              ttGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGc      <  1:327158/76‑1 (MQ=255)
                                                              ttGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGc      <  1:307483/76‑1 (MQ=255)
                                                              ttGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGc      <  1:278423/76‑1 (MQ=255)
                                                                ggTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGcac    >  2:852885/1‑76 (MQ=255)
                                                                 gTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGcaca   >  2:311856/1‑76 (MQ=255)
                                                                 gTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGcaca   >  2:1279296/1‑76 (MQ=255)
                                                                  tAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGcacac  >  1:489692/1‑76 (MQ=255)
                                                                         |                                                                    
TTTGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACAC  >  REL606/1100954‑1101095

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