Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,911,997 (A)6→7 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,9911.A97.7% 174.1 / ‑2.7 44intergenic (‑192/‑131)yifB/ilvLpredicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (20/23);  total (20/24)
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.50e-01

GACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATT  >  REL606/3911916‑3912062
                                                                            |                                                                       
gACAAAGGATGAACTACGAGGAAGGGAACAACCTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTAttttt                                                                          >  1:1076001/1‑76 (MQ=255)
   aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaa                                                                      >  2:78308/1‑76 (MQ=255)
   aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaa                                                                      >  2:1738686/1‑76 (MQ=255)
   aaaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaa                                                                      <  2:247174/76‑1 (MQ=255)
     aGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaa                                                                    <  2:1138617/76‑1 (MQ=255)
       gATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaaaa                                                                  <  2:1535404/76‑2 (MQ=255)
        aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTaaaaaaac                                                                 <  2:1125691/76‑3 (MQ=255)
            aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAacaaca                                                              <  1:1333651/76‑1 (MQ=255)
            aCTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAacaaca                                                              >  2:442919/1‑76 (MQ=255)
                 gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTta                                                         >  1:355878/1‑76 (MQ=255)
                  aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTtat                                                        <  1:820366/76‑1 (MQ=255)
                  aGGAAGGGAACAACATTCAGACTGAAATTTAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTgat                                                        <  2:487352/76‑1 (MQ=255)
                     aaGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATAtt                                                     <  1:340692/76‑1 (MQ=255)
                           acaacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAatt                                               >  2:238319/1‑76 (MQ=255)
                           acaacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAatt                                               <  2:1814805/76‑1 (MQ=255)
                            caacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAatta                                              <  2:76653/76‑1 (MQ=255)
                             aacaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAattat                                             <  2:1391311/76‑1 (MQ=255)
                              acaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAattatt                                            >  2:1172100/1‑76 (MQ=255)
                              acaTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAattatt                                            <  2:1588914/76‑1 (MQ=255)
                                aTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTaa                                          >  1:1730398/1‑76 (MQ=255)
                                aTTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTaa                                          >  1:1405010/1‑76 (MQ=255)
                                 ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTaaa                                         >  2:740408/1‑76 (MQ=255)
                                 ttCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTaaa                                         <  2:918810/76‑1 (MQ=255)
                                  tCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAAc                                        >  1:831259/1‑76 (MQ=255)
                                   cATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACg                                       >  1:1189169/1‑76 (MQ=255)
                                    aTACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGc                                      >  2:1414835/1‑76 (MQ=255)
                                      aCTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGcat                                    >  2:892557/1‑76 (MQ=255)
                                      aCTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGcat                                    >  2:1742329/1‑76 (MQ=255)
                                        tGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGcatca                                  >  2:1106339/1‑76 (MQ=255)
                                        tGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGcatca                                  <  2:1685528/76‑1 (MQ=255)
                                         gAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat                                 <  1:1106476/76‑1 (MQ=255)
                                          aaaTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                                <  1:50455/76‑1 (MQ=255)
                                          aaaTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATa                                <  2:231299/76‑1 (MQ=255)
                                            aTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                              <  2:994668/76‑1 (MQ=255)
                                            aTTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa                              >  1:657872/1‑76 (MQ=255)
                                             ttGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                             <  1:1578208/76‑1 (MQ=255)
                                             ttGAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaa                             <  1:644705/76‑1 (MQ=255)
                                              tGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt                           <  1:1659183/76‑1 (MQ=255)
                                               gAATTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAAt                           >  1:1850033/1‑76 (MQ=255)
                                                aaTTTTTTTCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATc                          <  1:941947/76‑1 (MQ=255)
                                                   ttttttCACTCACTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGc                       >  1:24110/1‑76 (MQ=255)
                                                        tcactcacTATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaa                  <  1:841620/76‑1 (MQ=255)
                                                         cactcactATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaa                 <  2:719614/76‑1 (MQ=255)
                                                          tctcactATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaaaaa                >  2:1913475/2‑76 (MQ=255)
                                                           ctcactATTCTATTTTTAACAAACAACACTTTATATTGAAATTATTAAACGCATCCTAACAATCAGCCACCAAata              >  2:914301/1‑76 (MQ=255)
                                                               ctATTTTATTTTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCt           >  1:948976/1‑76 (MQ=255)
                                                                tattttatttTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt          <  2:1773267/76‑1 (MQ=255)
                                                                tattttatttTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt          >  2:1640022/1‑76 (MQ=255)
                                                                tattttatttTTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt          <  1:1464794/76‑1 (MQ=255)
                                                                 attttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTg         <  1:152443/76‑1 (MQ=255)
                                                                  ttttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGt        <  1:1592461/76‑1 (MQ=255)
                                                                   tttattttTAAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTa       >  1:1267653/1‑76 (MQ=255)
                                                                       ttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAtt  <  1:292058/76‑1 (MQ=255)
                                                                       ttttaAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTAtt  <  2:905047/76‑1 (MQ=255)
                                                                            |                                                                       
GACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATT  >  REL606/3911916‑3912062

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