Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,991 (T)5→4 intergenic (‑192/‑131) yifB ← / → ilvL predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,911,9870T.91.1% 175.3 / 14.0 45intergenic (‑188/‑135)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base . (16/25);  total (17/28)
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.96e-01

CGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGT  >  REL606/3911914‑3912056
                                                                          |                                                                     
cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtattttattt                                                                     >  2:1110651/1‑76 (MQ=255)
 ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTattttatttt                                                                    <  1:396151/76‑1 (MQ=255)
  gACAAAGGATGAACTACCAGGAAGGGAACAACATTCATACTGAAATTGAAATTTTTTCACTCACTatttaattnta                                                                   >  2:1143756/1‑75 (MQ=255)
      aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAaaaa                                                               <  2:217326/76‑1 (MQ=255)
        ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaaaaac                                                             >  2:1042878/1‑75 (MQ=255)
           tGANCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTNTTTTCACTCACTATTTTATTTTAAAAaacaac                                                         <  2:1124806/76‑1 (MQ=255)
             aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAAcaacaa                                                       >  2:754846/1‑76 (MQ=255)
               cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAAtt                                                     <  1:658957/76‑1 (MQ=255)
                 aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTta                                                   <  1:1189202/76‑1 (MQ=255)
                 aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTta                                                   <  2:136871/76‑1 (MQ=255)
                 aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTta                                                   >  1:994655/1‑76 (MQ=255)
                  cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTtat                                                  <  1:229790/76‑1 (MQ=255)
                    aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTatat                                                >  2:1068798/1‑76 (MQ=255)
                        aGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGaa                                            <  2:86272/76‑1 (MQ=255)
                            aacaacATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAatt                                         <  2:309763/76‑1 (MQ=255)
                            aacaacATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAatta                                        <  2:267112/76‑1 (MQ=255)
                            aacaacATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAatta                                        <  2:18024/76‑1 (MQ=255)
                             acaacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAattat                                       <  1:262247/76‑1 (MQ=255)
                              caacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAattatt                                      >  1:436117/1‑76 (MQ=255)
                                acaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaa                                    >  1:705676/1‑76 (MQ=255)
                                 caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaaa                                   <  1:415802/76‑1 (MQ=255)
                                     cATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca                               >  2:148404/1‑76 (MQ=255)
                                     cATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca                               <  2:325406/76‑1 (MQ=255)
                                      aTACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat                              <  1:214607/76‑1 (MQ=255)
                                       tACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc                             >  1:506170/1‑76 (MQ=255)
                                        aCTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca                            >  1:488443/1‑76 (MQ=255)
                                         gtcaaaTTGAATTNTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat                           <  1:156759/73‑1 (MQ=255)
                                          tGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                          <  1:559936/76‑1 (MQ=255)
                                          tGAAATTGAATTTTTGTCACTCTCTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                          <  2:36188/76‑1 (MQ=255)
                                            aaaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGTAAGTATTAAACGCATCATaaa                        >  1:1133304/1‑76 (MQ=255)
                                            aaaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                        <  2:270070/76‑1 (MQ=255)
                                            aaaTTGAANTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                        <  1:1064312/76‑1 (MQ=255)
                                             aaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                       <  2:962073/76‑1 (MQ=255)
                                             aaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                       <  2:725032/76‑1 (MQ=255)
                                             aaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                       <  2:607879/76‑1 (MQ=255)
                                             aaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                       <  1:544640/76‑1 (MQ=255)
                                              aTTGAATTNTATTCACTCANTATTTTATTTTTAAAAAACAACNATTTATATTGAAATTATTAAACGCATCATaaaa                       <  2:1002840/76‑1 (MQ=255)
                                               ttGAATTTTTTTCACTCACTATTTTATTTTAAGAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt                     <  2:130120/76‑1 (MQ=255)
                                                tGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATc                    >  1:250765/1‑76 (MQ=255)
                                                        ttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa             <  2:253780/76‑1 (MQ=255)
                                                         tttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaa            >  1:845388/1‑76 (MQ=255)
                                                           tcactcacTATTTTATTTTATAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaaa          >  1:2983/1‑76 (MQ=255)
                                                           tcactcacTATTTTATTTTAAAAATCAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaaa          >  2:328798/1‑76 (MQ=255)
                                                              ctcactATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAata        <  1:240795/76‑1 (MQ=255)
                                                               tcactATTTTATTTTAAAAAGCCACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATc      >  1:643/1‑76 (MQ=255)
                                                               tcactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATc      >  1:94637/1‑76 (MQ=255)
                                                                 actATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt    >  2:678658/1‑76 (MQ=255)
                                                                  ctATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTg   <  2:529789/76‑1 (MQ=255)
                                                                   tattttatttTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTg   <  1:968767/76‑1 (MQ=255)
                                                                   tattttatttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt  >  2:1206592/1‑76 (MQ=255)
                                                                          |                                                                     
CGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGT  >  REL606/3911914‑3912056

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: