Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,100,824 A→T I136N (ATC→AAC)  hslU ← ATP‑dependent protease ATP‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,100,8240AT97.8% 164.5 / ‑4.4 45I136N (ATC→AAC) hslUATP‑dependent protease ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (16/28);  total (16/29)
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

GGAATGCCTGACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGATCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGATAGCCTGGACG  >  REL606/4100753‑4100895
                                                                       |                                                                       
ggAATGCCTGACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAg                                                                     <  1:758354/76‑1 (MQ=255)
ggAATGCCTGACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAg                                                                     <  2:804707/76‑1 (MQ=255)
  aaTGCCTGACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCa                                                                   <  1:1081699/76‑1 (MQ=255)
    tGCCTGACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACg                                                                 >  1:740299/1‑76 (MQ=255)
    tGCCTGACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACg                                                                 >  1:1579446/1‑76 (MQ=255)
      ccTGACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTc                                                               <  2:609778/76‑1 (MQ=255)
        tGACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCga                                                             >  2:64671/1‑76 (MQ=255)
         gACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCgag                                                            >  1:522298/1‑76 (MQ=255)
         gACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCgag                                                            <  2:1804082/76‑1 (MQ=255)
                   ggACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGttct                                                  <  2:1290117/76‑1 (MQ=255)
                   ggACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGttct                                                  <  1:623909/76‑1 (MQ=255)
                          tCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCa                                           >  2:218201/1‑76 (MQ=255)
                          tCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCa                                           >  1:1060075/1‑76 (MQ=255)
                             tgctgctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGtt                                        <  2:1530893/76‑1 (MQ=255)
                             tgctgctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGtt                                        <  2:1515905/76‑1 (MQ=255)
                             tgctgctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGtt                                        <  1:1395932/76‑1 (MQ=255)
                              gctgctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCCGttc                                       >  1:703845/1‑76 (MQ=255)
                               ctgctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGttct                                      >  2:926982/1‑76 (MQ=255)
                               ctgctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGttct                                      <  2:1515781/76‑1 (MQ=255)
                                 gctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGttcttc                                    <  2:630454/76‑1 (MQ=255)
                                 gctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGttcttc                                    <  2:1393230/76‑1 (MQ=255)
                                  ctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCa                                   <  2:989158/76‑1 (MQ=255)
                                  ctgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCa                                   >  2:261641/1‑76 (MQ=255)
                                   tgTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAg                                  <  1:261641/76‑1 (MQ=255)
                                    gTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAgc                                 <  2:241560/76‑1 (MQ=255)
                                        ggTCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGa                             <  1:80358/76‑1 (MQ=255)
                                          tCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATa                           >  2:1117733/1‑76 (MQ=255)
                                          tCTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATa                           <  1:1424827/76‑1 (MQ=255)
                                           cTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATaa                          <  1:1057621/76‑1 (MQ=255)
                                           cTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATaa                          <  2:835783/76‑1 (MQ=255)
                                           cTGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATaa                          >  1:242991/1‑76 (MQ=255)
                                            tGTCCCCAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAAc                         <  2:484569/76‑1 (MQ=255)
                                               ccccAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGt                      >  1:1296172/1‑76 (MQ=255)
                                                 ccAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGttt                    <  2:1089145/76‑1 (MQ=255)
                                                  cAGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGtttt                   >  2:1139442/1‑76 (MQ=255)
                                                   aGTTGTTTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGttttt                  <  2:950618/76‑1 (MQ=255)
                                                       gttTTTAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCg              <  2:366385/76‑1 (MQ=255)
                                                         ttttAGCAGGTGGGATCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGAt            <  1:275964/76‑1 (MQ=255)
                                                            tAGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGATAGc         >  2:1574675/1‑76 (MQ=255)
                                                             aGCAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGATAGcc        >  1:478885/1‑76 (MQ=255)
                                                               cAGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGATAGCCTg      <  2:516687/76‑1 (MQ=255)
                                                                aGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGATAGCCTgg     <  1:1610457/76‑1 (MQ=255)
                                                                aGGTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGATAGCCTgg     <  1:524595/76‑1 (MQ=255)
                                                                 ggTGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGATAGCCTGGa    >  1:1830893/1‑76 (MQ=255)
                                                                   tGGGTTCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGATAGCCTGGAcg  <  1:1420533/76‑1 (MQ=255)
                                                                       |                                                                       
GGAATGCCTGACGAGCAGCGGACGGTTCCTGCTGCTGTTCGGTCTGTCCCCAGTTGTTTTTAGCAGGTGGGATCAGCACGTCGAGAATACGTTCTTCTGCCAGTTCTTCAGCGCGATAACGGTTTTTCTCGATAGCCTGGACG  >  REL606/4100753‑4100895

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