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.0% 236.9 / ‑1.4 50intergenic (‑167/‑156)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (22/27);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA  >  REL606/3911893‑3912034
                                                                         |                                                                    
atgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttt                                                                    >  2:385774/1‑76 (MQ=255)
 tgatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                   <  2:480968/76‑1 (MQ=255)
 agatgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAtttt                                                                   >  1:1172186/2‑76 (MQ=255)
   atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACTAGGAAGGGAACCACAGNCATACTGAAATTGAAtttttt                                                                 >  2:279661/1‑76 (MQ=255)
   atgCTCGCAAAAATGCAGCGGACAAAGGATGAACTACTAGGAAGGGAACCACAGNCATAATGAAATTGACtttttt                                                                 >  2:279680/1‑76 (MQ=255)
      cTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAAttttttcac                                                              <  1:214347/76‑4 (MQ=255)
         gCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTtcactca                                                          >  2:18095/1‑76 (MQ=255)
            aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACta                                                       <  1:1033998/76‑1 (MQ=255)
             aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtat                                                      <  1:992248/76‑1 (MQ=255)
              aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatt                                                     >  1:1621702/1‑76 (MQ=255)
               aTGCAGCGGACAAAGGATGAANCACGCGGAAGGGAACATCATTCATACTGAAATTGAATTTTTTCACTCACtattt                                                    <  2:269460/76‑1 (MQ=255)
               aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattt                                                    <  1:1657498/76‑1 (MQ=255)
               aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattt                                                    <  1:28627/76‑1 (MQ=255)
                 gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttta                                                  <  1:275895/76‑1 (MQ=255)
                 gCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtatttta                                                  <  1:1066573/76‑1 (MQ=255)
                  cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACtattttat                                                 <  2:352898/76‑1 (MQ=255)
                      ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTAttttt                                             >  2:873784/1‑76 (MQ=255)
                         cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                          <  2:953324/76‑1 (MQ=255)
                         cAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaa                                          <  2:586227/76‑1 (MQ=255)
                          aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaa                                         >  2:633077/1‑76 (MQ=255)
                           aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTaaaaa                                        >  1:679414/1‑76 (MQ=255)
                             ggATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaac                                      <  1:455879/76‑1 (MQ=255)
                              gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAaaca                                     <  1:170615/76‑1 (MQ=255)
                                 gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                  <  1:1252735/76‑1 (MQ=255)
                                 gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAacaaca                                  <  2:1173713/76‑1 (MQ=255)
                                  aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                 >  2:1432086/1‑76 (MQ=255)
                                  aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAAcaacaa                                 >  2:1282727/1‑76 (MQ=255)
                                     tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAAttt                              >  1:949676/1‑76 (MQ=255)
                                      aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTta                             <  2:1324192/76‑1 (MQ=255)
                                        gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata                           >  2:1369186/1‑76 (MQ=255)
                                        gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata                           >  2:146041/1‑76 (MQ=255)
                                        gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata                           <  1:1371135/76‑1 (MQ=255)
                                        gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTtata                           <  2:221075/76‑1 (MQ=255)
                                         aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                          <  1:1267173/76‑1 (MQ=255)
                                         aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTatat                          >  2:203903/1‑76 (MQ=255)
                                           gAAGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                        <  1:1051201/76‑1 (MQ=255)
                                            aaGGGAACAACATTCATACTGAAATTGGATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                       >  2:191794/1‑76 (MQ=255)
                                            aaGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGa                       >  1:75225/1‑76 (MQ=255)
                                             aGGGAACAACATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaa                      <  1:1526946/76‑1 (MQ=255)
                                                 aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                  <  1:1035936/76‑1 (MQ=255)
                                                 aacaacATTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatta                  >  2:461073/1‑76 (MQ=255)
                                                  acaacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                 >  1:838724/1‑76 (MQ=255)
                                                  acaacaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                 >  2:1573867/1‑76 (MQ=255)
                                                     acaTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa              <  1:1407300/76‑1 (MQ=255)
                                                      caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTaa              >  2:541341/1‑76 (MQ=255)
                                                       aTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc            <  2:1471653/76‑1 (MQ=255)
                                                       aTTCATACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAAc            <  2:595932/76‑1 (MQ=255)
                                                           aTACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat        <  1:385774/76‑1 (MQ=255)
                                                           aTACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat        >  2:607382/1‑76 (MQ=255)
                                                           aTACTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat        >  1:993487/1‑76 (MQ=255)
                                                             aCTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatca      >  1:227902/1‑76 (MQ=255)
                                                              cTGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat     <  2:804847/76‑1 (MQ=255)
                                                               tGAAATTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa    >  2:898504/1‑76 (MQ=255)
                                                                 aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  >  2:1438307/1‑76 (MQ=255)
                                                                 aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  >  1:1173052/1‑76 (MQ=255)
                                                                 aaaTTGAATTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa  <  2:1334479/76‑1 (MQ=255)
                                                                         |                                                                    
ATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAA  >  REL606/3911893‑3912034

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