Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,219,546 T→G I97S (ATT→AGT)  yjbA → phosphate‑starvation‑inducible protein PsiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,219,5460TG94.8% 225.9 / 1.1 58I97S (ATT→AGT) yjbAphosphate‑starvation‑inducible protein PsiE
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base G (35/20);  total (37/21)
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.83e-01

CTGATTGTGAAGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCATTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGCTGCGATCCTGTTACT  >  REL606/4219473‑4219618
                                                                         |                                                                        
cTGATTGTGAAGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGtg                                                                        <  1:297542/76‑1 (MQ=255)
 tGATTGTGAAGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCATtgt                                                                       >  2:1759502/1‑76 (MQ=255)
 tGATTGTGAAGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCATtgt                                                                       >  1:1719820/1‑76 (MQ=255)
  gATTGTGAAGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGtgtg                                                                      >  2:1148875/1‑76 (MQ=255)
  gATTGTGAAGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGtgtg                                                                      <  1:1630059/76‑1 (MQ=255)
    ttGTGAATTTCTTTCAGTCCGTTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTgcg                                                                    <  1:400820/76‑1 (MQ=255)
     tgtgAAGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTgcgc                                                                   >  2:164812/1‑76 (MQ=255)
       tgAAGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCtt                                                                 >  1:1052482/1‑76 (MQ=255)
       tgAAGTACTTTCAGTCCGGTTTTCACGTCCCGCTCCGTTGCTTTGTCTATTTTGGGATCACCGCCAGTGTGCGCtt                                                                 <  2:668019/76‑1 (MQ=255)
         aaGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGa                                                               >  2:1535593/1‑76 (MQ=255)
          aGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGat                                                              <  1:588428/76‑1 (MQ=255)
            tACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGatca                                                            >  2:309689/1‑76 (MQ=255)
              cTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGatcatc                                                          >  2:397663/1‑76 (MQ=255)
              atttCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGatcatc                                                          >  1:959283/2‑76 (MQ=255)
                ttCAGTCCTGTTTTCCCTTCCCGCTCGGTTACTTTGTCTATATTGGGATCGCCGCCATTGTGCGCTTGATCAtcgt                                                        <  1:131686/76‑1 (MQ=255)
                ttCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATGGGGATCACCGCCAGTGTGCGCTTGATCAtcgt                                                        <  2:1315814/76‑1 (MQ=255)
                  cAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCAtcgtcg                                                      >  1:1405485/1‑76 (MQ=255)
                  cAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCAtcgtcg                                                      >  1:258294/1‑76 (MQ=255)
                   aGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGa                                                     >  1:492209/1‑76 (MQ=255)
                     tCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATc                                                   >  2:1229621/1‑76 (MQ=255)
                        ggTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATa                                                <  2:1263172/76‑1 (MQ=255)
                         gTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATaa                                               <  2:666408/76‑1 (MQ=255)
                         gTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATaa                                               <  2:1577770/76‑1 (MQ=255)
                          ttttCACTTCCCTCTACGTTACTTTGTCTATATTGGGATCTCCGCCAGTGTGCGCTTGATCATCGTCGATCATaaa                                              <  2:1628783/76‑1 (MQ=255)
                          ttttCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATaaa                                              >  2:367041/1‑76 (MQ=255)
                           tttCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAAt                                             <  2:1704398/76‑1 (MQ=255)
                           tttCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAAt                                             <  1:1633155/76‑1 (MQ=255)
                             tCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGGCGATCATAAATCg                                           >  1:919170/1‑76 (MQ=255)
                                cTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCa                                        >  2:1170819/1‑76 (MQ=255)
                                 ttCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCAc                                       >  2:266907/1‑76 (MQ=255)
                                  tCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACt                                      >  2:1764932/1‑76 (MQ=255)
                                  cccccCTACGTTACTTTGTCTATATTGGAATCCCCGCCAGTGTGCGCTGGATCTTCGTCGGTCATAAATCGCCACt                                      <  2:914061/75‑1 (MQ=255)
                                    ccGCTGCGTTTCTTTGTATATATGGGGATCCCCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTgg                                    <  1:828611/76‑1 (MQ=255)
                                     cGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGa                                   >  1:58101/1‑76 (MQ=255)
                                     cGCTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGa                                   >  2:1334396/1‑76 (MQ=255)
                                       ctgcGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGAtg                                 >  1:1340384/4‑76 (MQ=255)
                                       cTACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGAtg                                 <  2:1088046/76‑1 (MQ=255)
                                        tACGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGAtgt                                <  1:1444825/76‑1 (MQ=255)
                                         aCGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGAtgtg                               >  1:768545/1‑76 (MQ=255)
                                         aCGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGAtgtg                               >  1:1328372/1‑76 (MQ=255)
                                         aCGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGAtgtg                               >  1:1294326/1‑76 (MQ=255)
                                         aCGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGAtgtg                               >  1:969634/1‑76 (MQ=255)
                                          cGTTACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATgtgt                              >  2:1171323/1‑76 (MQ=255)
                                            ttACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTg                            >  2:1507137/1‑76 (MQ=255)
                                             tACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGa                           >  1:1657011/1‑76 (MQ=255)
                                             tACTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGa                           >  1:373638/1‑76 (MQ=255)
                                              aCTTTGTTTATGTTGGAATCGCCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGAt                          <  1:480342/76‑1 (MQ=255)
                                              aCTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGAt                          >  1:566332/1‑76 (MQ=255)
                                              aCTTTGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGAt                          <  2:211186/76‑1 (MQ=255)
                                                 ttGTCTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTa                       <  1:1375466/76‑1 (MQ=255)
                                                     cTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCg                   >  2:1566703/1‑76 (MQ=255)
                                                     cTATATTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCg                   >  1:783860/1‑76 (MQ=255)
                                                      tataTTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCgg                  <  1:701519/76‑1 (MQ=255)
                                                      tataTTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCgg                  >  1:28953/1‑76 (MQ=255)
                                                       atatTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGc                 >  1:713242/1‑76 (MQ=255)
                                                       atatTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGc                 >  2:248926/1‑76 (MQ=255)
                                                        tatTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGCt                >  2:1304530/1‑76 (MQ=255)
                                                        tatTGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGCt                >  1:1412930/1‑76 (MQ=255)
                                                           tGGGATCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGCTGCg             >  1:884076/1‑76 (MQ=255)
                                                               aTCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGCTGCGATcc         >  2:553730/1‑76 (MQ=255)
                                                                tCACCGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGCTGCGATCCt        >  2:1099761/1‑76 (MQ=255)
                                                                   ccGCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCCACTCGGCTGCGATCCTGtt     <  2:1469719/76‑1 (MQ=255)
                                                                     gCCAGTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGCTGCGATCCTGTTAc   <  1:1002577/76‑1 (MQ=255)
                                                                      ccAGTGTGCGCTGGAGCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCGACTCGGCTGCGATCCTGTTACt  <  2:319895/76‑1 (MQ=255)
                                                                         |                                                                        
CTGATTGTGAAGTACTTTCAGTCCGGTTTTCACTTCCCGCTACGTTACTTTGTCTATATTGGGATCACCGCCATTGTGCGCTTGATCATCGTCGATCATAAATCGCCACTGGATGTGTTGATCTACTCGGCTGCGATCCTGTTACT  >  REL606/4219473‑4219618

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