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.91.5% 191.4 / 14.5 47intergenic (‑188/‑135)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (23/20);  total (25/22)
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.96e-01

AGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTA  >  REL606/3911912‑3912057
                                                                           |                                                                      
aGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtattttat                                                                        >  2:1407735/1‑76 (MQ=255)
 ncGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCNCTCACtattttatt                                                                       >  2:220239/2‑76 (MQ=255)
  cGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtattttattt                                                                      >  2:838290/1‑76 (MQ=255)
   ggACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTattttatttt                                                                     <  2:506175/76‑1 (MQ=255)
       aaaGGATGAACTANGAGGAAGGGAACAACAATCATACTGAAATTGTATTTTTTTCACTCACTATTTTATTTTaaaa                                                                 <  2:1403582/76‑1 (MQ=255)
       aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTaaaa                                                                 >  1:777579/1‑76 (MQ=255)
        aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAaaaa                                                                <  2:1418856/76‑1 (MQ=255)
         aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTCCACTCACTATNTTATTTTaaagag                                                               >  1:206334/1‑73 (MQ=255)
            aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaa                                                           >  2:841827/1‑76 (MQ=255)
             tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaac                                                          <  1:1024757/76‑1 (MQ=255)
             tGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAaacaac                                                          <  2:366981/76‑1 (MQ=255)
               aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAAcaacaa                                                        <  1:993952/76‑1 (MQ=255)
                   aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTta                                                    <  1:1478594/76‑1 (MQ=255)
                    cGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTtat                                                   >  2:767539/1‑76 (MQ=255)
                     gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTtata                                                  >  2:702468/1‑76 (MQ=255)
                       ggAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATAtt                                                <  1:678985/76‑1 (MQ=255)
                         aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGa                                              >  2:415145/1‑76 (MQ=255)
                         aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGa                                              <  1:992251/76‑1 (MQ=255)
                         aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGa                                              <  1:291297/76‑1 (MQ=255)
                          aGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGaa                                             <  2:508031/76‑1 (MQ=255)
                          aGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGaa                                             >  2:1136424/1‑76 (MQ=255)
                            ggAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAat                                           <  2:470363/76‑1 (MQ=255)
                             gAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAatt                                          <  1:1407736/76‑1 (MQ=255)
                              aacaacATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAatta                                         >  2:80104/1‑76 (MQ=255)
                               acaacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAattat                                        <  1:860439/76‑1 (MQ=255)
                               acaacaTTCATACTGAAATTGAATTTTTTTCACNCACTATTTTATTTTAAAAAACANCAATTTATATNGAAattat                                        >  2:1289183/1‑76 (MQ=255)
                                 aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAAttatta                                      >  2:968453/1‑76 (MQ=255)
                                  acaTTCATACTGAGATTGAATTTCTTTCACTCACTATTTTATTTTAAAAAACAAGAATTTATATTGAAATTATTaa                                     <  1:1318560/76‑1 (MQ=255)
                                  acaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaa                                     <  1:335570/76‑1 (MQ=255)
                                  acaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaa                                     >  2:664759/1‑76 (MQ=255)
                                   caTTCATACTGAAATTGAATTTTTTTCNCTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaaa                                    >  1:1522186/1‑76 (MQ=255)
                                   caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaaa                                    >  1:224122/1‑76 (MQ=255)
                                   caTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaaa                                    <  1:398426/76‑1 (MQ=255)
                                   caTTCATACTGAAATGGAATTTGTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTaaa                                    <  2:53267/76‑1 (MQ=255)
                                     ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACg                                  <  2:1420063/76‑1 (MQ=255)
                                     ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACg                                  <  2:701038/76‑1 (MQ=255)
                                       cATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGca                                >  1:1470355/1‑76 (MQ=255)
                                              aaaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                         >  2:1210531/1‑76 (MQ=255)
                                              aaaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                         >  1:1113332/1‑76 (MQ=255)
                                              aaaTTGAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                         <  1:810332/76‑1 (MQ=255)
                                                   gAATTTTTTTCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCg                    <  1:1136518/76‑1 (MQ=255)
                                                        tttttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaa               >  2:24598/1‑76 (MQ=255)
                                                          tttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAGTTATTAAACGCATCATAAAAATCGGCCaaaa             >  1:614591/1‑76 (MQ=255)
                                                          tttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaa             >  1:313909/1‑76 (MQ=255)
                                                           ttCACTCACTATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaa            >  2:827831/1‑76 (MQ=255)
                                                             cactcactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAat          >  1:1013814/1‑76 (MQ=255)
                                                             cactcactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAat          >  1:669899/1‑76 (MQ=255)
                                                               ctcactATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAANGCATCATAAAAATCGGCCAAAAAatat        <  1:1533442/76‑1 (MQ=255)
                                                                tcactATTTNATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATc       >  1:33159/1‑76 (MQ=255)
                                                                 cactATTTTATTTTAAACAACAGCAATTTATATTGAAATTATTAAACGCATCATAAAAATCTGCCAAAAAATATCt      >  2:233652/1‑76 (MQ=255)
                                                                  actATTTTATTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt     >  1:275670/1‑76 (MQ=255)
                                                                    tattttatttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt   <  2:654275/76‑1 (MQ=255)
                                                                    tattttatttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt   <  1:747607/76‑1 (MQ=255)
                                                                    tattttatttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt   <  1:386306/76‑1 (MQ=255)
                                                                    tattttatttTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt   >  1:240879/1‑76 (MQ=255)
                                                                     attttattttAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTa  >  1:737319/1‑76 (MQ=255)
                                                                           |                                                                      
AGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTA  >  REL606/3911912‑3912057

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: