Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,616,302 (G)7→6 intergenic (‑42/‑354) purL ← / → yfhD phosphoribosylformylglycinamidine synthase/putative periplasmic binding protein of transport system; b2558_1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,616,2960G.94.7% 259.7 / 8.5 57intergenic (‑36/‑360)purL/yfhDphosphoribosylformylglycinamidine synthase/putative periplasmic binding protein of transport system; b2558_1
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base . (31/23);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

ATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACC  >  REL606/2616221‑2616367
                                                                           |                                                                       
aTGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAg                                                                         <  2:939039/76‑1 (MQ=255)
aTGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAg                                                                         >  1:1367557/1‑76 (MQ=255)
aTGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAg                                                                         <  1:932048/76‑1 (MQ=255)
 tcccGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAgg                                                                        <  1:727272/74‑1 (MQ=255)
 tGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAgg                                                                        >  1:955807/1‑76 (MQ=255)
 tGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAgg                                                                        >  2:357562/1‑76 (MQ=255)
  gCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAggg                                                                       <  1:438024/76‑1 (MQ=255)
     gACCGTGCCGGCGAACCCCGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAgggggg                                                                    <  1:751214/76‑1 (MQ=255)
        aGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGaaa                                                                 <  1:1528478/76‑1 (MQ=255)
         gTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGaaac                                                                >  1:1627434/1‑75 (MQ=255)
         gTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGaaac                                                                <  1:240499/76‑2 (MQ=255)
          tGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTTCCCAAGGGGGGAAACg                                                              <  2:1147057/76‑1 (MQ=255)
          tGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACg                                                              >  1:1166773/1‑76 (MQ=255)
          tGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACg                                                              >  1:755363/1‑76 (MQ=255)
          tGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACg                                                              <  2:1267318/76‑1 (MQ=255)
           gCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACgg                                                             <  1:1525080/76‑1 (MQ=255)
            cAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACggg                                                            >  2:599031/1‑76 (MQ=255)
             aGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGc                                                           <  2:145847/76‑1 (MQ=255)
             aGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGc                                                           >  2:587216/1‑76 (MQ=255)
               gCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGt                                                         >  2:337160/1‑76 (MQ=255)
               gCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGt                                                         <  1:213505/76‑1 (MQ=255)
                 gAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTCCCCAAGGGGGGAAACGGGCGTCa                                                       <  2:682530/76‑1 (MQ=255)
                 gAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCa                                                       <  2:981392/76‑1 (MQ=255)
                 gAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCa                                                       <  2:994179/76‑1 (MQ=255)
                  aaCCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCAt                                                      >  1:1241639/1‑76 (MQ=255)
                     cACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATtat                                                   <  1:497144/76‑1 (MQ=255)
                     cACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATtat                                                   <  1:990117/76‑1 (MQ=255)
                        gCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATaaa                                                <  1:1024749/76‑1 (MQ=255)
                        gCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATaaa                                                >  2:1776358/1‑76 (MQ=255)
                         cAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAg                                               <  2:1690769/76‑1 (MQ=255)
                            aaTTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAAt                                            >  1:667761/1‑76 (MQ=255)
                              tttCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCt                                          <  1:907931/76‑1 (MQ=255)
                               ttCAATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTg                                         >  2:1622661/1‑76 (MQ=255)
                                   atcatAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATgc                                     >  1:826071/1‑76 (MQ=255)
                                     catAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCCTTATAAAGAATCTGATgcgc                                   >  2:416324/1‑76 (MQ=255)
                                      atAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCt                                  >  2:1430265/1‑76 (MQ=255)
                                       tAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTg                                 >  2:1429051/1‑76 (MQ=255)
                                            ctctcGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGaa                            >  2:1087924/1‑76 (MQ=255)
                                            ctctcGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGaa                            <  2:1347823/76‑1 (MQ=255)
                                            ctctcGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGaa                            <  2:344122/76‑1 (MQ=255)
                                              ctcGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAAc                          >  1:198865/1‑76 (MQ=255)
                                              ctcGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAAc                          >  2:1812753/1‑76 (MQ=255)
                                               tcGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAAcc                         <  2:361240/76‑1 (MQ=255)
                                                cGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCg                        >  2:263612/1‑76 (MQ=255)
                                                cGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCg                        >  1:1821524/1‑76 (MQ=255)
                                                 gTCTTCTTAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGt                       >  2:1138890/1‑76 (MQ=255)
                                                 gTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGt                       >  2:1256329/1‑76 (MQ=255)
                                                  tcttctAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGtt                      <  2:762550/76‑1 (MQ=255)
                                                  tcttctAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGtt                      >  1:1444216/1‑76 (MQ=255)
                                                     tctAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGc                   >  1:363225/1‑76 (MQ=255)
                                                      ctAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCg                  >  1:934636/1‑76 (MQ=255)
                                                      ctAAGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCg                  >  1:1214928/1‑76 (MQ=255)
                                                        aaGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTg                <  2:1719814/76‑1 (MQ=255)
                                                        aaGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTg                <  2:307974/76‑1 (MQ=255)
                                                         aGCTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTgg               <  2:355449/76‑1 (MQ=255)
                                                           cTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGaa             <  1:335422/76‑1 (MQ=255)
                                                           cTTTCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGaa             >  1:968863/1‑76 (MQ=255)
                                                            tttCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGaaa            >  1:900819/1‑76 (MQ=255)
                                                            tttCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGaaa            >  1:822967/1‑76 (MQ=255)
                                                             ttCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAAt           >  1:652521/1‑76 (MQ=255)
                                                             ttCGGTGTACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAAt           <  1:1187957/76‑1 (MQ=255)
                                                                  tgtACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAAt      >  2:1297074/1‑76 (MQ=255)
                                                                  tgtACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAAt      <  1:750788/76‑1 (MQ=255)
                                                                    tACCCAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCa    >  2:471053/1‑76 (MQ=255)
                                                                      cccAAGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCAcc  >  1:312352/1‑76 (MQ=255)
                                                                           |                                                                       
ATGCCGACAGTGCAGGCGAACCACGCAGAATTTCCATCATAAGTCTCTCGTCTTCTAAGCTTTCGGTGTACCCAAGGGGGGGAAACGGGCGTCATTATAAAGAATCTGATGCGCTGACGAAACCGTTTGCGTGGAAATAAAATCACC  >  REL606/2616221‑2616367

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