Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,972 (T)7→8 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,9651.T98.1% 202.4 / ‑2.7 52intergenic (‑166/‑157)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (28/23);  total (29/23)
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.73e-01

AATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA  >  REL606/3911892‑3912034
                                                                          |                                                                      
aaTGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtt                                                                      <  1:1545921/76‑1 (MQ=255)
 atgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttt                                                                     >  1:1280002/1‑76 (MQ=255)
 atgatgCTCACAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttt                                                                     <  1:854612/76‑1 (MQ=255)
  tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                   >  1:218224/1‑76 (MQ=255)
   gatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttt                                                                  <  1:990370/76‑1 (MQ=255)
   gatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttt                                                                  <  2:1201231/76‑1 (MQ=255)
    atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttt                                                                 >  1:483995/1‑76 (MQ=255)
      gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttttttt                                                               <  2:402614/76‑2 (MQ=255)
       cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttttc                                                              <  2:1440809/76‑3 (MQ=255)
        tCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttttca                                                             <  1:467644/76‑4 (MQ=255)
           cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactc                                                           <  2:426267/76‑1 (MQ=255)
            aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtcactca                                                          >  2:133164/1‑76 (MQ=255)
                aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtat                                                      >  1:216367/1‑76 (MQ=255)
                aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtat                                                      >  1:270611/1‑76 (MQ=255)
                 tGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatt                                                     <  1:65783/76‑1 (MQ=255)
                   cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtatttt                                                   <  2:1402705/76‑1 (MQ=255)
                      cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttatt                                                >  1:212464/1‑76 (MQ=255)
                      cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttatt                                                >  1:422690/1‑76 (MQ=255)
                      cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttatt                                                >  2:681693/1‑76 (MQ=255)
                       ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACtattttattt                                               <  1:344628/76‑1 (MQ=255)
                        gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTattttatttt                                              >  2:221383/1‑76 (MQ=255)
                         aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTAttttt                                             <  1:363228/76‑1 (MQ=255)
                          cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTa                                            >  1:473499/1‑76 (MQ=255)
                              ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTaaaaa                                        <  1:246402/76‑1 (MQ=255)
                                  gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAaacaa                                    <  1:813074/76‑1 (MQ=255)
                                     cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                 >  2:1752318/1‑76 (MQ=255)
                                      tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                >  1:948585/1‑76 (MQ=255)
                                      tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAt                                >  1:1454830/1‑76 (MQ=255)
                                          aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            >  1:883792/1‑76 (MQ=255)
                                          aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtat                            <  2:484423/76‑1 (MQ=255)
                                             aaGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATAtt                         >  1:1167211/1‑76 (MQ=255)
                                              aGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        <  2:1156786/76‑1 (MQ=255)
                                              aGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        >  2:666763/1‑76 (MQ=255)
                                              aGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        <  2:10901/76‑1 (MQ=255)
                                              aGGGAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        <  1:268364/76‑1 (MQ=255)
                                                ggAACAACATTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                      <  1:1349750/76‑1 (MQ=255)
                                                   acaacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                   <  1:983107/76‑1 (MQ=255)
                                                    caacaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                  >  1:278867/1‑76 (MQ=255)
                                                      acaTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                >  1:1019121/1‑76 (MQ=255)
                                                       caTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta               <  1:529457/76‑1 (MQ=255)
                                                        aTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa              >  1:922301/1‑76 (MQ=255)
                                                        aTTCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa              >  1:160268/1‑76 (MQ=255)
                                                         ttCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa             >  2:721774/1‑76 (MQ=255)
                                                         ttCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa             >  1:342233/1‑76 (MQ=255)
                                                         ttCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa             <  1:199672/76‑1 (MQ=255)
                                                         ttCATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaaa             <  1:1429276/76‑1 (MQ=255)
                                                           cATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAGAAAACAACAATTTATATTGAAATTATTAAACg           >  1:316747/1‑76 (MQ=255)
                                                           cATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg           >  1:1441424/1‑76 (MQ=255)
                                                           cATACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACg           >  1:1115179/1‑76 (MQ=255)
                                                             tACTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca         >  2:528356/1‑76 (MQ=255)
                                                              aCTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat        >  1:1550740/1‑76 (MQ=255)
                                                              aCTGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTNATATTGAAATTATTAAACGcat        <  2:244044/76‑1 (MQ=255)
                                                                tGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca      <  2:843755/76‑1 (MQ=255)
                                                                tGAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca      <  2:1108150/76‑1 (MQ=255)
                                                                 gAAATTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat     >  2:264186/1‑76 (MQ=255)
                                                                  aaaTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa    <  2:473499/76‑1 (MQ=255)
                                                                  aaaTTGAATTTCTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa    >  2:390780/1‑76 (MQ=255)
                                                                   aaTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa   <  2:617699/76‑1 (MQ=255)
                                                                   aaTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa   <  1:636352/76‑1 (MQ=255)
                                                                    gttGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  >  1:477513/2‑76 (MQ=255)
                                                                    aTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  <  1:527373/76‑1 (MQ=255)
                                                                    aTTGAATTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  >  1:477499/1‑76 (MQ=255)
                                                                          |                                                                      
AATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA  >  REL606/3911892‑3912034

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