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.88.5% 225.0 / 29.4 61intergenic (‑188/‑135)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base . (26/28);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAC  >  REL606/3911916‑3912058
                                                                       |                                                                       
gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAttttt                                                                     >  1:1486095/1‑76 (MQ=255)
  cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaa                                                                   >  2:1169646/1‑76 (MQ=255)
   aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaaa                                                                  <  1:1055833/76‑1 (MQ=255)
   aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaaa                                                                  >  2:583846/1‑76 (MQ=255)
    aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAaaaa                                                                 <  1:1168567/76‑1 (MQ=255)
     aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaaaaa                                                                >  1:126993/1‑76 (MQ=255)
     aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaaaaa                                                                >  1:1407511/1‑76 (MQ=255)
       gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaaca                                                             <  1:818548/76‑1 (MQ=255)
         tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaac                                                           <  1:188373/76‑1 (MQ=255)
          gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAacaaca                                                          <  1:1482764/76‑1 (MQ=255)
           aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAAcaacaa                                                         >  1:428616/1‑76 (MQ=255)
           aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAAcaacaa                                                         <  1:187150/76‑1 (MQ=255)
               aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTta                                                     >  2:34091/1‑76 (MQ=255)
               aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTta                                                     >  2:148280/1‑76 (MQ=255)
               aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTta                                                     >  2:382704/1‑76 (MQ=255)
                cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTAACTATTTTATTTTAAAAAACAACAATTtat                                                    >  2:229146/1‑76 (MQ=255)
                 gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTtata                                                   <  1:86780/76‑1 (MQ=255)
                   ggAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATAtt                                                 >  2:792246/1‑76 (MQ=255)
                    gAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTg                                                <  2:1459231/76‑1 (MQ=255)
                    gAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTg                                                <  2:543470/76‑1 (MQ=255)
                    gAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTg                                                <  1:1079303/76‑1 (MQ=255)
                     aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGa                                               <  1:760717/76‑1 (MQ=255)
                     aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGa                                               <  1:24973/76‑1 (MQ=255)
                       gggAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                                              <  1:212333/76‑1 (MQ=255)
                            caacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAattatt                                        <  1:261640/76‑1 (MQ=255)
                            caacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAattatt                                        <  2:139001/76‑1 (MQ=255)
                            caacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAattatt                                        <  1:579845/76‑1 (MQ=255)
                              acaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaa                                      >  2:660219/1‑76 (MQ=255)
                              acaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaa                                      >  2:809574/1‑76 (MQ=255)
                                 ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACg                                   <  2:347232/76‑1 (MQ=255)
                                 ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACg                                   <  1:732612/76‑1 (MQ=255)
                                 ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACg                                   <  2:138989/76‑1 (MQ=255)
                                   cATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca                                 >  1:953284/1‑76 (MQ=255)
                                    aTACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat                                >  2:1484404/1‑76 (MQ=255)
                                    aTACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat                                >  2:1177022/1‑76 (MQ=255)
                                    aTACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat                                >  2:1177012/1‑76 (MQ=255)
                                     tACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc                               >  1:223031/1‑76 (MQ=255)
                                        tGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                            >  1:786995/1‑76 (MQ=255)
                                         taaaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa                           >  2:148402/2‑76 (MQ=255)
                                           aaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                         >  1:417063/1‑76 (MQ=255)
                                            aTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa                        <  1:898945/76‑1 (MQ=255)
                                            aTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa                        >  2:572534/1‑76 (MQ=255)
                                            aTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa                        >  2:572529/1‑76 (MQ=255)
                                              tGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATc                      >  2:614248/1‑76 (MQ=255)
                                               gAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCg                     >  1:504906/1‑76 (MQ=255)
                                                  tttttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGcc                  <  1:1468262/76‑1 (MQ=255)
                                                   ttttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                 <  1:1210200/76‑1 (MQ=255)
                                                   ttttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                 <  1:287385/76‑1 (MQ=255)
                                                   ttttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                 >  1:1318682/1‑76 (MQ=255)
                                                   ttttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                 <  1:287401/76‑1 (MQ=255)
                                                    tttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaa                <  1:1343962/76‑1 (MQ=255)
                                                    tttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaa                <  1:1567998/76‑1 (MQ=255)
                                                     ttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa               <  1:981895/76‑1 (MQ=255)
                                                      tttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaa              <  1:1169184/76‑1 (MQ=255)
                                                         cactcactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAat           <  1:1030989/76‑1 (MQ=255)
                                                         cactcactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAat           >  2:1027482/1‑76 (MQ=255)
                                                          actcactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAata          >  1:1578131/1‑76 (MQ=255)
                                                             cactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCt       >  2:40673/1‑76 (MQ=255)
                                                              actATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt      >  2:233018/1‑76 (MQ=255)
                                                              actATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt      <  2:129778/76‑1 (MQ=255)
                                                              actATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt      <  2:746092/76‑1 (MQ=255)
                                                               ctATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTg     >  2:1558334/1‑76 (MQ=255)
                                                                tattttatttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGACCAAAAAATATCTTGt    >  1:1635699/1‑76 (MQ=255)
                                                                tattttatttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGACCAAAAAATATCTTGt    >  1:732685/1‑76 (MQ=255)
                                                                 attttattttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt    <  2:1417917/76‑1 (MQ=255)
                                                                 attttattttAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTa   >  1:947176/1‑76 (MQ=255)
                                                                  ttttattttaAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAc  <  1:1574004/76‑1 (MQ=255)
                                                                  ttttattttaAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAc  >  1:1580854/1‑76 (MQ=255)
                                                                  ttttattttaAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAc  >  1:307567/1‑76 (MQ=255)
                                                                       |                                                                       
GACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAC  >  REL606/3911916‑3912058

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