Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,997 (A)6→7 intergenic (‑198/‑125) 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,9911.A80.0% 92.0 / 17.5 30intergenic (‑192/‑131)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base A (10/14);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATT  >  REL606/3911919‑3912062
                                                                            |                                                                         
aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaa                                                                        >  2:270083/1‑76 (MQ=255)
 aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaa                                                                       >  2:2645/1‑76 (MQ=255)
   ngATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaaa                                                                   <  1:1004409/75‑1 (MQ=255)
   ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaaa                                                                   <  1:635835/76‑1 (MQ=255)
    gATGAACTACGAGGTAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaaaa                                                                  <  2:747464/76‑2 (MQ=255)
           tACGAGGGAGGGAGCAACATTCATACTGATATTGCATTTTTTTCACTCACTATTTTATTTTTAAAAATGCAACAAt                                                            <  2:1032078/76‑1 (MQ=255)
           tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTATTTTCACTCACTATTTTATTTTTAAAAAAACAACAAt                                                            >  1:840060/1‑76 (MQ=255)
             cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTCTTAAAAAAacaaca                                                              <  1:31175/72‑1 (MQ=255)
             cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAacaaca                                                              >  2:31175/1‑72 (MQ=255)
             cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAAttt                                                          <  1:1131028/76‑1 (MQ=255)
                  aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAAAAAAAATTTTTAtt                                                     >  2:714258/1‑76 (MQ=255)
                     ggAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGaa                                                  <  1:200711/76‑1 (MQ=255)
                      gAACAACACTCATACTGAAATTGAATTTTTTTCATTCACTATTTTATTCTTAAAAAAACAACGATTTATATTGaaa                                                 >  2:1030757/1‑76 (MQ=255)
                          aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAattat                                             <  2:507663/76‑1 (MQ=255)
                            caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAAttatta                                           >  1:1128498/1‑76 (MQ=255)
                             aTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTaa                                          <  2:1148780/76‑1 (MQ=255)
                             aTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTaa                                          <  1:437857/76‑1 (MQ=255)
                              ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTaaa                                         >  1:617932/1‑76 (MQ=255)
                                  tACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGca                                     >  2:844429/1‑76 (MQ=255)
                                    cTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGcatc                                   <  2:518153/76‑1 (MQ=255)
                                     tGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGcatca                                  <  1:710758/76‑1 (MQ=255)
                                     tGAAATTGAATTTTTTTCACCGACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGcatca                                  <  2:973930/76‑1 (MQ=255)
                                       aaaTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                                <  2:332400/76‑1 (MQ=255)
                                         aTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                              <  2:668277/76‑1 (MQ=255)
                                            gAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt                           >  1:684227/1‑76 (MQ=255)
                                                  ttttCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                     <  1:218631/76‑1 (MQ=255)
                                                  tattCACTCACTATTTTATTTATAAAAAAACATCAATTTATATTGAAATTATTAATCGCATCATAAAAATCGGCCa                     <  1:140512/74‑1 (MQ=255)
                                                     tcactcacTATTTTATCTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaa                  <  2:486802/76‑1 (MQ=255)
                                                      ccctcactATTTTATTTTTAAAAAAACAACAATTTATATTGCAATTATTAAAAGCATCATAAAAATCGGCCaaaaa                 >  2:403353/3‑76 (MQ=255)
                                                      cactcactATTTTATTTTTAAAAAAACCACAATTTATATTGAAATTATTAAACGCATCATAAAAATCCTCCaaaaa                 <  2:1036171/76‑1 (MQ=255)
                                                      cactcactATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaa                 >  2:66131/1‑76 (MQ=255)
                                                           actATTTTATTTTTAAAACAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATc            >  2:1048875/1‑76 (MQ=255)
                                                             tattttatttTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt          >  1:634018/1‑76 (MQ=255)
                                                              attttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTg         >  2:360167/1‑76 (MQ=255)
                                                               ttttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt        <  2:1147380/76‑1 (MQ=255)
                                                               ttttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt        <  2:730109/76‑1 (MQ=255)
                                                                  atttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa    >  1:435277/3‑76 (MQ=255)
                                                                  atttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa    >  1:824685/3‑76 (MQ=255)
                                                                   nttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAt   <  1:1002455/75‑1 (MQ=255)
                                                                   nttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAt   <  1:1002450/75‑1 (MQ=255)
                                                                    ttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAtt  >  1:537519/1‑76 (MQ=255)
                                                                            |                                                                         
AAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATT  >  REL606/3911919‑3912062

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