Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,816,986 A→T intergenic (+385/‑172) ygdH → / → sdaC hypothetical protein/predicted serine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,816,9860AT97.5% 306.7 / ‑0.5 79intergenic (+385/‑172)ygdH/sdaChypothetical protein/predicted serine transporter
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base T (34/43);  total (34/45)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ATTCTCTTTTACTTTTAATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGATCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTA  >  REL606/2816912‑2817058
                                                                          |                                                                        
aTTCTCTTTTACTTTTAATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGAt                                                                         <  2:1110886/76‑1 (MQ=255)
   ctctTTTACTTTTAATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCaa                                                                      <  1:727048/76‑1 (MQ=255)
   ctctTTTACTTTTAATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCaa                                                                      >  1:196264/1‑76 (MQ=255)
   ctctTTTACTTTTAATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCaa                                                                      <  1:1695913/76‑1 (MQ=255)
   ctctTTTACTTTTAATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCaa                                                                      <  1:1337376/76‑1 (MQ=255)
     cttttacttttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGa                                                                    >  2:446419/1‑76 (MQ=255)
     cttttacttttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGa                                                                    >  1:1234623/1‑76 (MQ=255)
      ttttacttttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGAt                                                                   <  1:1127674/76‑1 (MQ=255)
      ttttacttttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGAt                                                                   <  2:1707327/76‑1 (MQ=255)
        ttacttttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCa                                                                 >  2:1377513/1‑76 (MQ=255)
        ttacttttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCa                                                                 <  2:1543473/76‑1 (MQ=255)
           cttttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTg                                                              >  1:362368/1‑76 (MQ=255)
            ttttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTga                                                             >  2:1959556/1‑76 (MQ=255)
            ttttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTga                                                             >  2:292221/1‑76 (MQ=255)
             tttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTgat                                                            >  1:634312/1‑76 (MQ=255)
             tttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTgat                                                            <  2:1022685/76‑1 (MQ=255)
             tttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTgat                                                            >  1:2246766/1‑76 (MQ=255)
              ttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTgata                                                           <  2:1436143/76‑1 (MQ=255)
              ttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTgata                                                           <  2:1523393/76‑1 (MQ=255)
              ttaATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTgata                                                           >  2:2233364/1‑76 (MQ=255)
               taATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTgatag                                                          >  1:1722402/1‑76 (MQ=255)
                   ttACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATAc                                                      <  1:116409/76‑1 (MQ=255)
                    tACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACa                                                     >  1:1317777/1‑76 (MQ=255)
                    tACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACa                                                     <  1:666839/76‑1 (MQ=255)
                     aCAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACAt                                                    <  1:325579/76‑1 (MQ=255)
                     aCAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACAt                                                    >  1:416729/1‑76 (MQ=255)
                     aCAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACAt                                                    >  1:517257/1‑76 (MQ=255)
                     aCAGAGGAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACAt                                                    <  1:1540175/76‑1 (MQ=255)
                      cAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATa                                                   <  2:217928/76‑1 (MQ=255)
                       agagTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATaa                                                  <  2:379509/76‑1 (MQ=255)
                           tAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTtg                                              >  1:1576353/1‑76 (MQ=255)
                           tAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTtg                                              >  2:219876/1‑76 (MQ=255)
                             aaGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTtgtg                                            <  2:1115825/76‑1 (MQ=255)
                             aaGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTtgtg                                            <  2:927595/76‑1 (MQ=255)
                              aGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTtgtgt                                           >  1:1798426/1‑76 (MQ=255)
                                  ttGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTAt                                       >  1:678263/1‑76 (MQ=255)
                                  ttGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTAt                                       >  1:1569501/1‑76 (MQ=255)
                                  ttGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTAt                                       >  1:1848032/1‑76 (MQ=255)
                                  ttGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTAt                                       >  1:571420/1‑76 (MQ=255)
                                    gggATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCt                                     <  1:839535/76‑1 (MQ=255)
                                     ggATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCtt                                    <  2:905811/76‑1 (MQ=255)
                                     ggATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCtt                                    <  1:119762/76‑1 (MQ=255)
                                     ggATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCtt                                    <  1:1226928/76‑1 (MQ=255)
                                      gATAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCttt                                   <  2:1699162/76‑1 (MQ=255)
                                       aTAATCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTc                                  <  2:84493/76‑1 (MQ=255)
                                         aaTCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCg                                <  2:193916/76‑1 (MQ=255)
                                         aaTCTATCTTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCg                                <  2:1461211/76‑1 (MQ=255)
                                            ctatctTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGccc                             >  1:2024729/1‑76 (MQ=255)
                                             tatctTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCt                            <  1:999971/76‑1 (MQ=255)
                                              atctTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTc                           >  1:981404/1‑76 (MQ=255)
                                              atctTCCAAGTAGATTATTGTATTTGGGATCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTc                           <  1:2091480/76‑1 (MQ=255)
                                               tctTCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCa                          <  1:108326/76‑1 (MQ=255)
                                                 ttCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCaaa                        <  2:441264/76‑1 (MQ=255)
                                                 ttCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCaaa                        >  2:1403672/1‑76 (MQ=255)
                                                  tCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAat                       <  2:483263/76‑1 (MQ=255)
                                                  tCCAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAat                       >  1:36798/1‑76 (MQ=255)
                                                   ccAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAatt                      <  2:1753996/76‑1 (MQ=255)
                                                   ccAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAatt                      <  2:2105383/76‑1 (MQ=255)
                                                   ccAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAatt                      <  2:1640841/76‑1 (MQ=255)
                                                   ccAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAatt                      <  2:2217810/76‑1 (MQ=255)
                                                   ccAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAatt                      <  2:1068752/76‑1 (MQ=255)
                                                    cAAGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAatta                     >  1:1181066/1‑76 (MQ=255)
                                                     aaGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAattat                    >  2:911796/1‑76 (MQ=255)
                                                     aaGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAattat                    <  1:19117/76‑1 (MQ=255)
                                                      aGTAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAattatt                   <  2:2198249/76‑1 (MQ=255)
                                                        tAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTAc                 >  2:2139104/1‑76 (MQ=255)
                                                        tAGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTAc                 >  2:1050356/1‑76 (MQ=255)
                                                         aGATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTAcg                <  2:510161/76‑1 (MQ=255)
                                                          gATTATTGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTAcgg               <  1:301491/76‑1 (MQ=255)
                                                            ttattGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTAcggcg             >  1:221743/1‑76 (MQ=255)
                                                              attGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGt           >  1:596245/1‑76 (MQ=255)
                                                               ttGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTa          >  1:1351448/1‑76 (MQ=255)
                                                               ttGTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTa          <  2:1502236/76‑1 (MQ=255)
                                                                 gTATTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTaaa        >  1:690669/1‑76 (MQ=255)
                                                                   aTTTGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATg      >  1:1366276/1‑76 (MQ=255)
                                                                    tttGGGTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGa     <  2:1291974/76‑1 (MQ=255)
                                                                       gggTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTa  <  2:46499/76‑1 (MQ=255)
                                                                       gggTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTa  <  2:491384/76‑1 (MQ=255)
                                                                       gggTTCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTa  <  1:1307576/76‑1 (MQ=255)
                                                                          |                                                                        
ATTCTCTTTTACTTTTAATTTACAGAGTAAAGCGTTGGGATAATCTATCTTCCAAGTAGATTATTGTATTTGGGATCAAGATCACTGATAGATACATAACTTGTGTGTATCTTTCCGCCCTCAAATTATTACGGCGGTAAATGATTA  >  REL606/2816912‑2817058

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