Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,972 (T)7→6 intergenic (‑173/‑150) 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,9660T.98.3% 277.9 / ‑1.4 58intergenic (‑167/‑156)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (27/30);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAT  >  REL606/3911891‑3912037
                                                                           |                                                                       
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                         <  2:990094/76‑1 (MQ=255)
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                         >  2:1644763/1‑76 (MQ=255)
gAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAt                                                                         <  1:765684/76‑1 (MQ=255)
 aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt                                                                        >  1:1116119/1‑76 (MQ=255)
     atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                    <  1:1152227/76‑1 (MQ=255)
       gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttca                                                                  <  1:766233/76‑3 (MQ=255)
        cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttcac                                                                 >  2:1081296/1‑73 (MQ=255)
           gCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactca                                                             <  2:1071791/76‑1 (MQ=255)
              aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACta                                                          >  1:610061/1‑76 (MQ=255)
              aaaaTGCAGCAGACAAAGGATGAACTACGAGGAAAGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACta                                                          <  1:986939/76‑1 (MQ=255)
              aaaaTGCAGCAGACAAAGGATGAACTACGAGGAAAGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACta                                                          <  1:986931/76‑1 (MQ=255)
               aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtat                                                         >  2:259405/1‑76 (MQ=255)
               aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtat                                                         >  2:1897587/1‑76 (MQ=255)
                  tGCAGCGGACAGAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttt                                                      <  1:1360337/76‑1 (MQ=255)
                    cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttat                                                    <  2:5715/76‑1 (MQ=255)
                     aGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttatt                                                   >  1:1562594/1‑76 (MQ=255)
                       cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTattttatttt                                                 <  2:95315/76‑1 (MQ=255)
                       cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTattttatttt                                                 <  2:686739/76‑1 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTattttatttta                                                >  2:1134768/1‑75 (MQ=255)
                        ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTAttttt                                                <  2:1433320/76‑1 (MQ=255)
                         gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTa                                               <  1:654061/76‑1 (MQ=255)
                          aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaa                                              <  2:1635034/76‑1 (MQ=255)
                           cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                             <  1:519654/76‑1 (MQ=255)
                            aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaa                                            >  1:1331667/1‑76 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaa                                           >  1:1439957/1‑76 (MQ=255)
                             aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaa                                           >  1:1439974/1‑76 (MQ=255)
                               ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaaa                                          <  1:388113/76‑1 (MQ=255)
                               ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaac                                         <  1:793324/76‑1 (MQ=255)
                               ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaac                                         <  2:1539153/76‑1 (MQ=255)
                               ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaac                                         <  1:1032325/76‑1 (MQ=255)
                                 aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                       >  2:91916/1‑76 (MQ=255)
                                   gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                     <  1:172356/76‑1 (MQ=255)
                                    aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                    <  1:567541/76‑1 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  >  2:1875719/1‑76 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  <  2:997112/76‑1 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  >  1:743713/1‑76 (MQ=255)
                                      cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAtt                                  >  1:411855/1‑76 (MQ=255)
                                       tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                                 <  1:62909/76‑1 (MQ=255)
                                        aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                                <  1:988691/76‑1 (MQ=255)
                                        aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                                >  2:1825810/1‑76 (MQ=255)
                                        aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                                <  1:1810482/76‑1 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                             >  2:297180/1‑76 (MQ=255)
                                           aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                             >  1:1311584/1‑76 (MQ=255)
                                            ggAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATAtt                            <  2:476076/76‑1 (MQ=255)
                                              aaGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                          >  2:770198/1‑76 (MQ=255)
                                              aaGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                          >  2:1556599/1‑76 (MQ=255)
                                              aaGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                          >  2:1556580/1‑76 (MQ=255)
                                                gggAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaaa                        <  2:1493611/76‑1 (MQ=255)
                                                gggAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaaa                        <  2:1493633/76‑1 (MQ=255)
                                                   aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                     >  1:44194/1‑76 (MQ=255)
                                                   aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                     >  2:1828169/1‑76 (MQ=255)
                                                   aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                     >  2:215779/1‑76 (MQ=255)
                                                     caacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                   >  2:1685676/1‑76 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                  <  1:274911/76‑1 (MQ=255)
                                                      aacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                  >  2:268263/1‑76 (MQ=255)
                                                       acaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa                 <  1:72466/76‑1 (MQ=255)
                                                        caTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa                <  1:827647/76‑1 (MQ=255)
                                                         aTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc               >  1:552006/1‑76 (MQ=255)
                                                           tCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGc             <  1:796552/76‑1 (MQ=255)
                                                             aTACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat           <  2:1810855/76‑1 (MQ=255)
                                                              tACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc          >  1:890947/1‑76 (MQ=255)
                                                              tACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc          <  2:1302133/76‑1 (MQ=255)
                                                                cTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat        <  1:1443080/76‑1 (MQ=255)
                                                                     aTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa   >  1:1518050/1‑76 (MQ=255)
                                                                      ttGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt  >  2:1851984/1‑76 (MQ=255)
                                                                           |                                                                       
GAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAT  >  REL606/3911891‑3912037

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: