Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,997 (A)6→5 intergenic (‑198/‑125) 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,9920A.97.1% 156.5 / ‑1.2 34intergenic (‑193/‑130)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (18/15);  minor base T (1/0);  total (19/15)
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.40e-01

ACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTA  >  REL606/3911917‑3912064
                                                                           |                                                                        
aCAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTa                                                                          <  1:814099/76‑1 (MQ=255)
  aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaa                                                                        >  2:870352/1‑76 (MQ=255)
    aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaa                                                                      <  2:106626/76‑1 (MQ=255)
       aTGAACTACGAGGAAGGGAGCAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaacaa                                                                   >  1:1030326/1‑73 (MQ=255)
         gAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAacaaca                                                                <  1:889122/76‑1 (MQ=255)
          aaCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAcaacaa                                                               <  2:1161790/76‑1 (MQ=255)
           aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAAt                                                              <  1:44872/76‑1 (MQ=255)
            cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAAtt                                                             <  1:385533/76‑1 (MQ=255)
             tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAAttt                                                            <  1:919698/76‑1 (MQ=255)
              aCGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTta                                                           >  1:1208125/1‑76 (MQ=255)
                gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTtata                                                         >  1:1027666/1‑76 (MQ=255)
                gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTtata                                                         <  2:255134/76‑1 (MQ=255)
                 aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAATTTTAAAAACAACTNTTTatat                                                        <  2:983498/76‑1 (MQ=255)
                        gAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAatt                                                 >  1:415814/1‑76 (MQ=255)
                         aacaacATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAatta                                                >  2:598509/1‑76 (MQ=255)
                               aTTCAACCTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAAc                                          >  1:1114566/1‑76 (MQ=255)
                                ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACg                                         <  2:642606/76‑1 (MQ=255)
                                   aTACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGcat                                      <  1:358327/76‑1 (MQ=255)
                                    tACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGcatc                                     >  2:867145/1‑76 (MQ=255)
                                     aCTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGcatca                                    <  2:631062/76‑1 (MQ=255)
                                      cTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGcatcat                                   <  1:470054/76‑1 (MQ=255)
                                       tGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                                  >  1:1026090/1‑76 (MQ=255)
                                        gAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa                                 >  1:631204/1‑76 (MQ=255)
                                        gAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa                                 <  1:698860/76‑1 (MQ=255)
                                        gAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa                                 <  1:829244/76‑1 (MQ=255)
                                          aaTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                               <  2:747887/76‑1 (MQ=255)
                                             tGAATTTTTTTCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATc                            >  1:424374/1‑76 (MQ=255)
                                                  ttttttCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa                       >  1:1097531/1‑76 (MQ=255)
                                                   tttttCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaa                      >  2:1041988/1‑76 (MQ=255)
                                                    ttttCACTCACTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa                     >  1:422968/1‑76 (MQ=255)
                                                       tcactcacTATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaaa                  >  2:70035/1‑76 (MQ=255)
                                                        cactcactATTTTATTTTTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAat                 >  2:1018552/1‑76 (MQ=255)
                                                                 ttttattttTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAc        >  1:950226/1‑76 (MQ=255)
                                                                 ttttattttTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTAc        <  1:748884/76‑1 (MQ=255)
                                                                  tttattttTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACt       >  2:593181/1‑76 (MQ=255)
                                                                  tttattttTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACt       >  1:581346/1‑76 (MQ=255)
                                                                   ttattttTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTa      >  1:664014/1‑76 (MQ=255)
                                                                     tattttTAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAt     >  2:1070789/3‑76 (MQ=255)
                                                                        ttttAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTa  <  1:711332/76‑1 (MQ=255)
                                                                        ttttAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTa  <  1:554869/76‑1 (MQ=255)
                                                                        ttttAAAAACAACAATTTATATAGAAANTANTGNACGCATCNTACAAACCGGCCAAAAAATATCTTGTACTATTTa  <  1:1052711/76‑1 (MQ=255)
                                                                           |                                                                        
ACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTA  >  REL606/3911917‑3912064

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